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ABSTRACT 

This  study  is  an  evaluation  of  a pressurized-water  reactor  (PWR)  1 
accident  as  defined  by  WASH  1400  for  the  proposed  nuclear  reactor  site 
at  Cementon,  Using  an  extension  of  the  Environmental  Pro- 

tection Agency's  AIREM  computer  code,  the  following  were  analyzed  for 

1 

up  to  50  miles  for  16  compass  directions:  (^)  whole  body  doses  due  to 

'J. 

cloud  submersion,  inhalation  and  groimd  deposition;  and  (^)  thyroid, 
lung,  and  gastrointestinal  tract  doses  due  to  inhalation.  In  addition, 
a model  was  developed  to  simulate  population  evacuation  up  to  both  10 
and  50  miles  and  again  determine  doses.  Subsequently,  the  nuclide  con- 
centration in  drinking  water,  milk,  and  fish  in  the  area  of  interest 
is  examined,  as  well  as  doses  due  to  swimming  in  a nearby  contaminated 
lake.  T 

Iodine  131  is  the  isotope  of  concern,  being  the  major  contributor 
to  the  inhalation  dose  which  itself  contributes  60%  of  the  total  dose 


received  in  this  instance.^  Comparison  of  this  study  with  WASH-1400 
shows  agreement  with  respect  to  total  man-rem  and  individual  doses 
received,  and  resulting  acute  lethalities. 

Finally , "^evacuation , a significant  mode  of  dose  reduction's  ex- 

- 

amined#  and  an  alternate  plan  of  seeking  shelter  in  the  event  of  an 
accident  is  suggested  because  at  close  proximities  to  the  site,  the 
average  winds  would  be  expected  to  carry  the  radioactive  cloud  down- 
wind much  too  quickly  to  permit  successful  evacuation. 


PART  I 


INTRODUCTION 


With  the  controversy  and  public  concern  over  the  safety  of  nuclear 
power  reactors,  and  the  proposed  construction  of  a pressurized  water  ^ 
reactor  at  Cementon,  New  York,  it  is  relevant  to  examine  a severe  acci- 
dent in  which  a significant  portion  of  the  isotope  inventory  is  re- 
leased to  the  atmosphere.  Rather  than  examine  the  design  base  accident 
in  which  there  is  a large  loss  of  coolant,  the  engineered  scifeguards 
fimction  properly,  and  therefore,  the  core  does  not  melt,  but  leakage 
to  the  atmosphere  takes  place  at  a design  leak  rate,  the  less  probable 
PWR  1 accident  is  assumed  to  occur.  In  Ref.  2,  this  event  assumes  a 
failure  of  all  the  engineered  safeguards  resulting  in  a core  meltdown, 
a steam  explosion  upon  contact  of  molten  fuel  and  water  in  the  pressure 
vessel  causing  the  vessel  head  to  penetrate  containment  with  a subse- 
quently large  portion  of  the  radioactive  inventory  being  released  out- 
side containment.  No  mechanism  is  postulated  here  for  how  this  woxild 
take  place  since  the  assumption  is  made  that  the  event  does  occur. 

In  WASH- 1400,  with  the  use  of  fa\ilt-tree  analysis,  the  worst 
possible  events  are  postulated  and  their  probabilities  of  occurrence 
determined.  Since  neither  an  event  this  serious  has  ever  occurred  in 
the  nuclear  ind\istry,  these  statistical  studies  must  be  used  as  the 
next  best  alternative;  nor  is  there  taken  into  account  the  fact  that 
the  probability  of  a PWR  1 type  accident  for  the  average  case  is 

_7 

7x10  per  reactor-year,  and  for  the  worst  meteorological,  and  popula- 
tion conditions  which  combined  results  in  the  highest  number  of 
lethalities  is  7x10  or  one  accident  of  this  type  in  one  billion 
reactor  years  or  with  100  reactors  running,  one  in  10  million  years. 


9*. 


Although  this  probeJjility  is  much  smaller  than  other  man-made  or 
naturally-occurring  disasters,  the  psychology  of  the  American  public 
and  indeed  of  any  general  populace,  but  the  United  States  in  particular, 
first  of  all  has  no  concept  of  veiy  large  or  very  small  numbers  and 
second,  has  come  to  equate  the  word  "nuclear”  with  great  destructive 
power  ever  since  1945.  Since  the  industry  has  grown  in  secrecy,  only 
a small  fraction  of  the  population  understands  the  concepts,  benefits 
and  risks  of  nuclear  power.  With  the  emotional  vulnerability  of  the 
vast  majority  of  people  who  only  wish  a light  to  appear  when  they 
turn  it  on.  they  are  more  susceptible  to  emotional  argximents  than  to 
a logic  which  they  either  cannot  understand,  or  has  not  been  presented 
clearly.  Most  of  the  more  important  decisions  made  by  people  in 
general  are  ultimately  made  emotionally,  not  logically,  and  it  is 
this  concept  which  must  be  evident  to  both  sides  of  the  nuclear 
question . 

This  study,  however,  will  only  deal  further  with  the  psychology 
of  people  in  the  portion  dealing  with  population  evacuation  in  the 
case  of  an  accident.  The  central  idea  is  to  assess  th.,  short-term 
results  to  the  surroxmding  countryside,  population,  bodies  of  water, 

1 

and  milk  production,  and  to  determine  the  doses  to  individuals  at  i 

various  locations  from  the  reactor  sites  from  three  mechanisms:  in- 
halation, cloud  submersion,  and  groimd  deposition.  Since  acute  fatal- 
ities are  of  most  concern  in  this  study,  a 30-day  integrated  dose  due 
to  a single  short-term  exposure  is  used  for  inhalation  doses  for  the 
whole  body  and  lung,  and  essentially  infinite  (50  years)  doses  to  the 

— 'i 

thyroid  and  GI  tract.  Also,  the  population  is  assumed  to  be  a year  i 

1990  projection  in  a 50-mile  radius  from  Cementon,  and  the  prevailing  j 

I 

1 


2 


J 


winds  vised  are  those  on-site  readings  taken  while  writing  the  Cementon 
Construction  Stage  Environmental  Report.  The  release  rate  and  meteor- 
ology is  assumed  constant  over  both  space  and  the  period  of  release 
(several  hours)  and  the  study  concerns  itself  primcirily  with  dry 
weather  conditions.  It  also  did  not  assume  increased  scrubbing  out  of 
contamination  due  to  land  profile  (mountains,  trees,  buildings). 

The  area  under  study  is  shown  in  Figs.  1 and  2. 
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PART  II 


HISTORICAL  REVIEW 

The  basis  for  the  dose  calcxilation  is  the  application  of  diffu- 
sion theory  to  the  atmospheric  case.  In  this  regard,  the  state  of  the 
art  will  be  briefly  sketched. 

Atmospheric  diffusion  theories  first  developed  by  Pick  were 
attempts  to  use  macroscopic  models  of  the  molecular  processes  of 
heat  and  mementum  transfer  in  which  the  transfer  was  described  in 
terms  of  its  gradient  and  a coefficient  of  diffusivity.  These  early 
models  then  were  attempts  to  solve  the  differential  equations  and 
determine  the  coefficients  which  would  accurately  represent  the  physi- 
cal process.  The  Tnajor  difficulty  encountered  was  that  the  coefficients 
were  found  to  vary  with  the  size  of  the  process,  the  height  above  the 
ground,  and  the  physical  state  of  the  atmosphere. 

By  the  1920 's,  G.  I.  Taylor  had  introduced  a statistical-motion 
approach  to  diffusion  which  was  based  on  continuoxas  movement  of  particles. 
This  theory  did  not  consider  dynamic  and  thermodynamic  causes  of  tur- 
bxilence,  but  it  began  to  techniques  employed  in  many  models  used  today. 
Sutton's  equation  of  diffusion  are  based  on  this  approach.  Because  of 
the  complex  nature  of  the  atmosphere,  diffusion  theory  has  always  been 
linked  with  knowledge  of  the  processes  of  turbulence. 

Both  Sutton  and,  more  recently,  Pasquill,  have  devised  methods 
involving  Gaussian  interpolation  in  a formula  that  relates  downwind 
concentrations  to  the  vertical  and  horizontal  concentrations  in  a plume 
as  a function  of  time  and  distance.  Sutton's  equations  relate  concen- 
trations to  the  vertical  gradient  of  the  wind  speed,  atmospheric  tur- 
bulence and  surface  roughness  while  Pasquill 's  relate  concentrations 

i 


6 


to  wind  direction  fluctuations  and  wind  speed  only.  This  study  uses 
a form  of  Sutton’s  eqiiation,  the  derivation  of  which  will  be  developed 
later.  There  is  currently  work  being  done  in  the  area  of  the  earlier 
neglected  dynamic  and  thermodynamic  mechanisms  of  turbulence  which 
will  hopefully  provje  to  yield  the  next  improvements  in  diffusion 
techniques 


THEORY 


A.  Development  of  the  Diffusion  Equation 

Becaiise  this  study  is  dependent  on  and  limited  by  the  diffusion 
equation  used,  it  is  appropriate  to  show  how  this  equation  has  been 
used.  As  mentioned  previoiisly,  diffusion  according  to  the  gradient 
transport  theory  developed  by  Fick  is  proportional  to  the  local  con- 
centration gradient.  Pick's  theory  considers  fluid  motion  from  a 
fixed  co-ordinate  system  while  statistical  theories  generally  con- 
sider motion  following  fluid  particles.  The  fundamental  equation  of 
Fickian  diffusion  theory  is  for  isotropic  diffusion, 

= (4tr  Ktr^/2  g-r^/4Kt 

where  K is  a constant  diffusivity  coefficient,  Q is  the  source 
strength,  the  total  amount  released,  and  q is  the  mean  value  of  a 
conservative  air  property  per  unit  mass  of  air. 

For  the  non-isotropic  case 

q(x.y,z,t)  _ (4irtr®''^(K  K K J’^/^exp  f zl 

Q ■ 4t^K  K K^ 

L X y 2 J 

Pick’s  diffusion  equation  is  solved  statistically  with  a Gaussian 
distribution  which  has  been  described  as  a continuous  souirce  diffusion 
model  by  Sutton,  and  does  not  apply  for  short  diffusion  times  or  in- 
bomogeneovis  or  non- stationary  conditions  although  it  will  show  trends 
if  these  conditions  are  altered.  Average  plume  diffusion  formulas 
usually  assume  an  instantaneous  point  source  of  material  diffusing  in 
three  dimensions.  The  Gaussian  formula  for  an  instantaneous  point 
soxirce  is: 


r ’I 


x(x,y.z)  . 2^ 


•3/2 


r o'  r 

X *7  z 


exp 


+ JL_  + £_> 
2^2  2^2  2 
X y z 


Eqn.  A 


where  X is  the  average  concentration, r is  the  standard  deviation  of  the 
distribution  of  material  in  a plume  and  Ji.  is  the  average  wind  component. 

There  are  two  restriction,  however,  in  using  the  above  equation. 
First,  it  is  assumed  that  the  equation  is  separable  in  x,  y and  z,  that 
is,  there  are  no  cross-wind  components.  The  second  is  based  on  the 
concept  of  relative  diffxision  which  states  that  the  equation  cannot  be 
used  to  describe  the  spread  of  a single  puff  of  material.  This  can  be 


seen  by  the  following.  The  maximum  rate  of  diffusion  from  a fixed  axis 

2 


is  proportional  to  time  t since  the  cloud  can  diffuse  only  in  two  di- 
rections (downwind  diffusion  is  negligible  due  to  the  equal  concentra- 
tions in  that  direction).  Average  diffusion  from  the  center  of  mass  of 

3 

a puff,  however,  is  proportional  to  t since  there  are  three  possible 
directions  in  which  diffusion  can  take  place. 

A continuous  point-so\arce  diffusion  formula  is  obtained  by  assximing 
that  the  plume  is  composed  of  an  infinite  number  of  overlapping  average 
puffs  carried  ed.ong  the  x-axis  by  an  average  wind  u.  Again,  assuming 
there  is  no  diffusion  in  the  x direction,  equation  A is  integrated  over 
t from  0 to  «•  as. 


X(x,  y,  z)  _ (2  n*  O'  O'  u)  ^ exp 

= y Z 

Q* 


where  Q'  is  the  source  strength-time  rate  of  material  emission  from  a 
continuous  point  source.  The  point  source  is  asstjmed  to  be  located  at 
or  near  the  earth's  surface  which  presents  a physical  barrier  to  diffu- 
wion  and  as  in  heat  conduction  theory,  assimiing  an  image  source  located 
symmetrically  with  respect  to  the  ground  plane,  and  assximing  the 


I' 
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1 


receptor  is  at  ground  level. 


2^ 

Q' 


^ O'z'" 


e^qp 


2 «r 


2 


Eqn.  B 


where  h is  the  effective  elevation  of  the  source  above  the  ground  plane. 

2 


Statistically  the  above  equation  is  completely  determined  by  o' 
and  it  agrees  reasonably  well  with  much,  but  not  all  atmospheric  diffu- 
sion data;  for  example,  the  equation  breaks  down  in  the  case  of  severe 
wind  shear . 

Integrating  Equation  B with  respect  to  y from  ao  yields  what 


is  known  as  the  crosswind  integrated  concentration,' 


CWI 


exp 


-h‘ 


2 r 


z J 


Over  some  time  t,  the  mean  wind  u will  generally  shift.  With  an  average 
concentration  estimate  over  time  which  is  much  greater  than  over  which 


the  mean  wind  is  computed,  X can  be  multiplied  by  the  frequency  with 

X 


which  the  wind  blows  toward  a particular  location  during  the  period. 


11 


Dividing  this  result  by  the  distance  from  the  source  at  each  location 


for  each  sector  (where  each  sector  width  is 


2trr 


n 


■;  r = distance , n = 


number  of  sectors),  adding  the  effect  of  simple  radiological  decay,  and 


summing  the  wind  frequencies  at  each  sector  yields, 

% 

f exp  [ — — exp  (- 


|r  =(|-)  ^ f exp  exp  (-  Xty'(5 


Eqn.  C 


which  is  the  basic  equation  used  in  this  study.  The  Gaussian  depletion 
model  assumes  that  the  shape  of  the  concentration  profile  in  the  vertical 
direction  is  not  altered  by  deposition.  If  on  the  other  hand  there  is 
vertical  mixing  and  the  cloud  is  depleted  at  the  bottom,  the  point  of 
maximum  concentration  moves  up  the  cloud  as  t increases,  carrying 


10 


more  material  over  a larger  area,  reducing  the  ground  level  concentrations 
with  subsequent  reductions  in  doses  and  dose  levels.  This  assumption 
then  yields  conservative  results. 

The  diffusion  equation  spreads  the  activity  uniformly  across  a 
given  22%  degree  compass  sector.  As  long  as  the  release  occurs  over 
several  hours,  this  averaging  may  be  quite  realistic  since  the  wind, 
even  if  relatively  steady,  would  probably  show  minor  shifts  dxiring  the 
period.  If  this  is  the  case,  the  Gaussian  single  plume  model  will 
yield  high  results.  In  this  study,  the  time  is  assumed  long  enough  to 
allow  the  assiomption  of  a long  plume  to  be  formed  but  short  enough  to 
minimize  wind  shifts  out  of  the  respective  wind  sector  being  studied. 

The  computer  code  is  not  applicable  nor  is  there  a study  of  strong  wind 
shifts  or  wind  shears  that  the  plxjme  may  pass  through. 


The  MODAIREM  Prog 


The  major  portion  of  the  main  AIPEM  Program  used  is  described  in 
Reference  10.  The  major  additions  to  the  main  AIREM  Program  are  the 
calculations  of  doses  due  to  ground  deposition  of  various  nuclides,  and 
the  effect  of  evacuation  on  total  doses  to  the  individuals  and  popula- 
tions concerned.  The  groxind  deposition  dose  calculations  are  considered 
in  Section  Illi. 

The  evacuation  model  employed  is  as  described  in  Reference  2.  The 
MODAIREM  Code  reads  population  data  and  computes  individual  and  popula- 
tion doses  for  the  imevacuated  case.  In  this  code,  two  additional 
alternatives  are  possible:  evacuation  outside  the  outer  boimdary  of  the 
outermost  radius  under  consideration  (50  miles),  and  evacuation  beyond 


t 
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the  outer  boundary  of  the  middle  radii;is  (10  miles)  in  a direction  away 
from  the  impending  radioactive  cloud.  Therefore,  those  who  have  success- 
fully evacxiated  before  cloud  passage  do  not  contribute  to  the  population 
dose.  The  original  population  is  modified  by  the  factor  where, 

Fj  = A + (1-A)  e “ 

where  A = the  fraction  of  the  population  that  for  some  reason  does  not 


evacuate,  X=  the  evacuation  rate  in  days 


.693/1%  where  1%  is  the 


time  necessary  to  evacuate  half  of  the  remaining  population,  t^,  = the 
time  for  the  radioactive  cloud  to  move  to  mesh  point  i,  T = the  time 
between  the  awareness  of  the  core  melt  and  leakage,  and  TL=  the  time  lag 
necessary  to  interpret  the  data  and  issue  the  warning  to  evaciiate.  T 
and  TL  may  be  fully  or  partially  concurrent. 

With  the  use  of  the  evacuation  model,  the  code  then  recomputes  the 


dose  to  the  individual,  and  the  total  doses  from  the  pathways  desired. 

For  ground  deposition  doses  during  evacuation,  the  passage  of  the  cloud 
begins  the  delivery  of  dose  to  a receptor  and  continues  until  departxire 
from  the  contaminated  area.  In  the  inhalation  and  cloud  submersion 
case,  the  man-rem  dose  depends  on  the  population  remaining  at  the  time 
of  cloud  passage  while  the  ground  dose  is  a fvinction  of  the  people 
present  at  cloud  passage  and  those  present  at  the  end  of  a given  period 
of  interest.  Since  the  man-rem  dosage  due  to  ground  deposition  is  so 
small,  a linear  relationship  which  is  conservative  is  used  to  approximate 


■f  ' ) 


the  total  dose. 


C.  Subroutine  RIVCON 


This  sxibroutine  computes  the  concentrations  in  PCi/1  of  the 
radionuclides  listed  in  the  main  program  for  a river  flowing  in  the 


i 


vicinity  of  a nuclear  power  plant  which  has  undergone  an  accidental 

release  of  fission  products.  The  river  is  divided  into  reaches  which 

/ 

correspond  to  mesh  points  (MM  II),  and  the  concentration  of  a certain 
nuclide  in  a reach  is  the  sum  of  the  direct  deposition  on  the  reach, 
the  decayed  concentration  of  nuclide  that  has  travelled  from  the  previous 
reach,  and  the  concentration  discharged  into  the  river  by  the  released 
cooling  water. 

Two  major  assumptions  in  this  subroutine  are  that  the  deposited 
nuclides  mix  well  with  the  river  water  to  yield  a uniform  concentration, 
and  that  the  particle  size  is  small  enough  that  transportation  of  the 
deposited  nuclide  not  in  solution  is  well  approximated  by  solution 
transport.  Considering  the  lengths  of  the  reaches  involved  and  the 
dispersal  of  the  nuclides  over  large  amounts  of  river  surface  area,  the 
first  approximation  is  reasonable.  The  second  will  result  in  conserva- 
tively high  calculations  for  untreated  water.  The  effect  becomes  more 
noticeable  for  treated  water  since  water  treated  by  municipal  treatment 
facilities  has  all  or  nearly  all  of  its  suspended  nuclides  removed. 

For  this  reason,  the  solubility  of  the  nuclide  is  taken  into  account 
for  treated  water. 

The  program  takes  into  account  the  varying  times  and  velocities  of 
hi^  and  low  tide,  allows  up  to  19  reaches  and  up  to  20  pijfiles  (cross 
sections)  per  reach.  Up  to  11  depths  per  profile  are  used  in  a Simpson's 
Rule  integration  to  determine  the  river  cross  sectional  area  at  each 
location  from  which  flow  volumes  in  each  reach  are  determined. 


The  solution  transport  equation  is. 


where  C ^ = the  concentration  (pCi/1)  at  reach  x at  time  t,  Ax  = the 
distance  of  the  previous  reach,  At  = the  time  it  takes  to  travel  the 

/ 

previous  reach  and  arrive  in  the  current  reach,  = the  direct  depo- 

sition (pCi/sec)  on  the  river  and/or  the  discharge  of  the  plant  through 

cooling  water,  Q = the  streamflow  (1/sec)  in  reach  x at  time  t,  t = 

x,^ 

time  (days)  concentration  is  measured  after  contamination. 

Foiar  types  of  water  treatment  facilities  are  considered.  A Class  1 
treatment  plant  utilizes  only  water  chlorination  and  does  not  reduce 
dissolved  or  suspended  nuclide  amounts.  A Class  2 facility  utilizes 
filtration  and  chlorination  to  remove  all  suspended  but  no  dissolved 


radionuclides.  A Class  3 facility  involves  chemical  processing  that 
removes  variable  amoimts  of  dissolved  nuclides  and  nearly  all  suspended 
nuclides  and  a Class  4 plant  filters  and  chemically  processes  water  to 
remove  all  suspended  and  variable  amounts  of  dissolved  nuclides. 

Ihe  treated  water  concentration  for  nuclide  1 in  reach  i,  C., 

il  can 


be  represented  as: 


C.,  s F'W'.,  + F"  W"., 
il  il  il 


Where  W'.,  + the  dissolved  concentration  of  nuclide  1 in  reach  i, 
il 

W'\  = the  suspended  concentration  of  nuclide  1 in  reach  i,  F'  = the 

fraction  of  dissolved  nuclide  that  passes  through  the  treatment  plant, 

F"  = the  fraction  of  suspended  nuclide  that  passes  through  the  treatment 


plant. 


In  this  study,  F'  = 1 nuclides  for  which  no  experimental  data 
has  been  developed  to  the  contraiy.  For  plant  types  2,  3 and  4,  F"=  0 
and  for  types  1 and  2^ F'  = 1.  The  subroutine  also  sums  the  nuclide 
concentrations  by  reaches  before  and  after  water  treatment.  If  RIVCON 


is  not  desired,  one  blank  card  can  be  inserted  instead  of  the  data  cards 


i 

i 


When  computing  water  concentrations  in  river  reaches , the  time 
it  takes  a particle  of  water  to  travel  through  the  reach  must  be  deter- 
mined. To  do  so,  the  particle  must  be  moved  with  the  tide  up-  and 
down-stream  according  to  the  up-  and  down-stream  velocities  and  the 
time  of  ebb  and  flow  tides . RIVCON  follows  the  particle  from  the  first 
reach  at  the  outer  limits  of  interest  beginning  with  the  tide  going 
either  in  or  out  as  specified  to  the  last  reach  while  retaining  the 
times  of  travel  in  each  reach. 


In  each  reach,  the  vertical  cross  sectional  area  is  computed 
using  Simpson's  Rule  or. 


CROSS 


(yo  + 4yj  + 2y2  + 4yg  + 2yj^  + ...  + 


■'2m-l  ■’^2m 


where  yj^  = f(Xj^),  y^^  corresponding  to  profile  depths  and  to  each 

eqtiidistant  width.  The  interval  a to  b is  broken  into  2m  equal  parts 

of  width  A X = so  that  x = a,  x.  = a + x,  x-  + 2 x ...  x.  = b. 

2m  o 1 2 2m 


The  constant  a is  set  equal  to  zero  at  one  and  of  the  shore  of  the  river 

at  each  profile  and  b corresponds  to  the  river  width.  In  RIVCON,  from 

3 to  11  depths  can  be  \ised  to  chart  the  river  bottom  with  up  to  20 

profiles  per  reach  to  provide  as  much  detail  as  possible  to  find  the 

2 2 

average  cross  sectional  area.  Several  constants  used  are:  .0929  m /ft  , 
9.99972x10^  and  4.470193x10  ^ m/sec  per  mi/  hr. 

RIVCON  overlays  a river  on  a ground  deposition  pattern.  Sedimenta- 
tion, resuspension  and  the  effects  of  dams  are  not  included  in  this 
subprogram. 
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D.  Subroutine  LAKEDS 


This  subroutine  computes  the  treated  and  untreated  water  concentra- 
tions in  pCi/  1 of  the  radionuclides  listed  in  the  main  program  for  a 
lake  located  within  the  area  of  interest  sxirrounding  a nuclear  power 
plant  which  has  undergone  an  accident.  The  concentrations  are  used 
to  determine  individual  doses  due  to  swimming  and/or  boating  on  the 
contaminated  lake . 

The  lake  is  assxiraed  to  fit  primarily  in  one  mesh  point  (im,  II)  ’ 

with  mixing  of  the  deposited  nuclides  in  the  water  to  obtain  an  average 

I 

concentration  throughout  the  lake  water.  Small  particle  size  is  ] 

assumed  in  using  solution  transport  as  the  water  is  removed  to  a i 

treatment  facility  before  consumption.  Tne  subroutine  allows  for  decay 
of  the  nuclides  in  the  lake,  but  once  the  water  has  been  treated,  no 
further  decay  is  assumed  since  the  water  may  be  used  at  any  time. 

The  user  has  several  options  using  LAKEDS.  With  the  xise  of  one 
blank  card,  the  calculations  will  be  deleted  and  no  other  data  cards 
are  then  necessary.  When  using  LAKEDS,  the  doses  can  be  foimd  due  to 
swimming  and  boating  for  whatever  time  is  spent  in  or  near  the  water. 

Either  of  the  two  dose  calculations  can  be  deleted  if  desired.  For  each 
activity,  skin  and/or  whole  body  doses  can  be  calculated  and  either  may 
be  deleted  if  desired.  In  the  water  treatment  section,  entering  a 
class  1-4  will  result  in  concentrations  of  treated  water  being  displayed 
while  a 0 entered  will  assume  no  water  treatment  and  terminate  the 
subroutine . 

— In  LAKEDS,  the  whole  body  and  skin  dose  of  factors  for  swimming 
for  each  isotope  are  entered,  and  the  boating  dose  factors  are  computed 
as  one  half  of  the  emersion  dose  factors  due  to  the  geometry.  The  lake 

L - 


The  whole  body  dose  due  to  swimming,  the  whole  body  dose  due  to  boating. 


the  skin  dose  due  to  swimming,  and  the  skin  dose  due  to  boating  are 
found  with  equations  of  the  type: 

DOSE  (mrem)  = CON(pCi/l  )(Dose  factor  (time  spent  (hr)  in  activity) 

(pCi/1  ) 

As  in  RIVCON,  water  treatment  facility  usage  is  included. 

E.  Subroutines  FISHR  and  FISHL 

These  subroutines  compute  the  amount  of  various  nuclides  due  to 
their  concentration  in  fish.  FISHR  computes  the  concentrations  in 
fish  in  each  reach  of  a river  in  the  area  of  interest.  FISHL  does  the 
same  in  a lake. 

The  equation  used  for  each  reach  and  the  lake  is, 

C = CF  ( k) 

^LL  LL  ^ ^MM,II,LL 

where  Cj^  = the  concentration  of  nuclide  LL  in  fish  in  the  lake  or  a 

reach  in  the  river  (pCi/kg);  CF  = fresh  water  concentration  factor 

LL 

. . ll 

(pCi/kg  per  PCi/1  )(see  Table);  A = concentration  (pCi/1  ) of 

MM, II 

nuclide  LL  in  the  river  or  lake  located  in  mesh  point  MM, II. 

The  above  concentrations  assume  that  the  fish  do  not  move  signifi- 
cantly away  from  the  original  mesh  point  location  where  they  obtained 
the  concentration.  The  fish  are  assxmied  to  rapidly  assimilate  the 
nuclides  from  the  surroundings. 

Sedimentation  and  resuspension  are  implicitly  considered  in  the 
fish  ingestion  dose  calculations  through  the  use  of  the  reconcentration 
factors  for  fresh  water  fish. 
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F.  Subroutine  MILK 


This  subroutine  computes  the  concentration  of  radionuclides  in 
milk  due  to  the  ingestion  of  nuclides  by  the  dairy  cow  which  depends 
upon  its  mesh  point  location.  The  first  calculations  are  for  a given 
number  of  days  after  the  accident  specified  in  the  data.  MILK  allows 
either  the  evacuation  of  the  dairy  cow  after  one  day  of  exposure 
(equivalent  to  a change  to  uncontaminated  feed)  or  its  remaining  in 
place  during  the  period  of  interest  and  eating  the  contaminated 
vegetation.  Three  drinking  water  options  are  available:  river,  lake, 
or  underground  well.  Each  mesh. point  is  assumed  to  use  mderground 
wells  unless  entered  otherwise  as  data.  If  a mesh  point  area  uses 
river  water  and  is  not  located  on  the  river,  siibroutine  REACH  finds  the 
closest  reach  and  uses  water  from  that  location. 

If  the  evacuation  mode  is  used,  decay  is  allowed  according  to  the 


effective  half-life,  i.e.,  X 


eff  ~ ^radiological 


Tijbiological 


If  evacuation  is  not  used,  the  concentration  for  the  initial  period  is 
determined.  Subsequent  concentraticns  are  added  to  earlier  concentrations 
which  have  been  decayed  by  the  effective  half  life.  Concentrations  are 
printed  for  specified  increments  of  days  for  a specified  length  of  time 
unless  all  doses  are  less  than  10  pCi/1  at  which  time  the  subroutine 
will  terminate. 


MILK  assumes  that  there  is  no  decay  in  transmission  because  there 
is  rapid  transmission  of  water  from  its  origin  to  its  place  of  consump- 
tion, This  subroutine  considers  only  forage  crop  contamination  since 
it  is  2issumed  that  stored  feed  and  grain  would  be  covered  and  probably 
not  sold  or  used  if  contaminated.  Therefore,  the  use  of  100%  forage 
crop  is  conservative. 
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Ilia  aquation  for  nuclide  concentration  in  milk  is,  (see  Table  11) 

but  is  modified  due  to  accident  rather  than  long  term  analysis  and 
redefinition  of  some  terms  as 


k [fj  '“f’  * ^ "■f]  " v j« 


ingestion. 


Terns  similar  to  the  ones  enclosed  in  brackets  for  the  forage  crop  are 

found  for  grain  and  stored  feed  but  are  neglected  as  stated.  Likewise, 

any  concentration  due  to  inhalation  of  radionuclide  particles  is  very 

difficult  to  determine  but  is  much  smaller  than  concentrations  from 

31 

other  sources.  The  second  term  is  due  to  water  ingestion. 


If  A=  Q, 


'■f] 

1:  It 


^ pCi/1 
d pCi/day 

I and  B = L^W,  then  for  the  case  of  no  evacuation. 


(l-e“  '^®ei 


)6-  t Be- 


where  Xe  is 'the  effective  decay  in  the  cow  for  nuclide  i or. 


^ 

c.  Th 


.693  .693 


1 1 


where  and  T^  are  the  biological  and  radiological  half-lives 

respectively  for  nuclide  i,  and  '^e  is  the  effective  nuclide  decay 

®i 

on  the  ground  for  nuclide  i,  or: 

A _ .693  .693 

= 1% = 


'H 


and  t is  the  time  after  nuclide  deposition.  = forage  deposition 

2 2 
(pCi/ra  ) and  F-  = soil  deposition  (pCi/m  ). 


19 


_ Xe  ^ 

n»e  e term  describes  the  rate  of  loss  of  nuclide  i from 

vegetation  and  radioactive  decay  while  the  remainder  of  the  first  ex- 
ponential term  describes  the  i^^  nuclide  concentration  due  to  forage 
crop  intake,  and  the  last  term  describes  the  i^^  nuclide  concentration 
due  to  ingested  drinking  water.! 

For  the  case  of  evacuation  after  one  day  of  ingestion  of  contaminated 
forage  or  switching  to  uncontaminated  feed,  the  equation  is  identical 
to  the  previous  one  for  up  to  one  day,  then  it  becomes, 

C,  = S.I.e  “ ^®c.  ^ 
d d d 1 

Comparison  of  Equation  A with  a Federal  Radiation  Council  study 
89  90 

of  Sr  and  Sr  shows  that  it  is  initially  conservative  and  approaches 

22  131 

FRC  calculations  closely  by  ten  days.  Comparison  with  a study  of  I 
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in  milk  also  shows  this  equation  to  be  conservative,  both  maximum 
concentrations  peaking  approximately  five  days  after  contamination. 


G.  Inhalation  Dcses 


Factors  for  computing  the  inhalation  doses  to  the  whole  body,  Ixing, 
thyroid,  and  GI  tract  are  listed  in  Table  1.  These  factors  have  been 

3 

taken  from  Reference  2 and  converted  to  mrem/sec  per  Ci/m  using  a 

-4  3 

breathing  rate  of  2.2  x 10  m /sec.  The  MODAIPJIM  Code  multiplies 
these  factors  by  the  cloud  concentration  of  each  nuclide  yielding  doses 
froiveach  mesh  point  and  these  values  multiplied  by  the  population  per 
sector  yield  their  respective  man-rem  doses. 

Assumptions  made  are  as  in  Reference  2 with  respect  to  nuclide 


distribution  and  retention.  Vfhole  body  inhalation  doses  are  integrated 
over  30  days  and  are  conservative  since  using  an  LD_  of  400  R is  based 
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on  short  teiro  exposure.  The  total  dose  during  any  period  of  time  from  the 


assumed  rapid  intake  at  time  t=o  to  a period  t~f rom  isotope  j to  organ  i is , 

^i  \ R. .(t 

jL 
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5IE. .1. 


D..  = 
i: 


:)exp(-  Agfft)  dt 


M. 

1 


where  D. . = rem  to  organ  i from  isotope  j;  E,,  = f .F. .E.(QF). .(MF). . ; 

= intake  in  >>.Ci;  = organ  burden  in  = mass  of  organ  i in  grams; 

Ej  = energy  of  disintegration  mode  j in  Mev;.fj  = probability  per  disintegration 
of  mode  j's  production;  F^^  = fraction  of  energy  emitted  by  the  mode  absorbed 
in  organ  i;  QF  = quality  factor;  and  MF  = modifying  factor. 

The  compartmental  model  of  internal  dose  with  feedback  is  approximated 
with  an  exponential  reduction  of  organ  burden  by  both  the  radiological  and 
biological  half-life  resulting  in  an  effective  decay 


H.  Whole  Body  Cloud  Submersion  Dose 

2 14 

The  cloud  submersion  doses  of  WASH  1400  and  WASH  1258  were  compared 
and  the  WASH  1258  values  were  used  as  being  slightly  more  current.  Both  sets 
of  values  are  seen  to  be  quite  similar,  being  derived  from  a semi-infinte 
cloud  model.  The  2‘Tr  geometry  was  used  to  determine  dose  factors,  with  gamma 
radiation  equal  from  all  directions  and  beta  with  short  penetrating  power  irra- 
diating only  from  one  side.  The  dose  factors  have  been  corrected  for  the 

fractional  penetration  of  beta  and  gamma  radiations  to  the  appropriate  tissue 

-3 

depth  of  5 cm  for  whole  body  and  7x10  cm  for  skin. 
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The  WASH  1400  report  accepted  the  year  2000  study  with  respect  to 

\ 

dose  factors  since  it  was  the  best  data  available,  but  the  WASH  1258 

report  recalculated  the  dose  factors  using  the  latest  available  decay 

schemes.  With  respect  to  both  air  and  water  submersion  doses,  an 

external  dose  depends  on  the  penetrating  power  of  the  emissions  of 

the  nuclide  under  study.  Beta  and  gamma  radiation  which  could  not 
2 

penetrate  7mg/cra  of  tissue  was  not  considered  in  the  skin  dose.  Both 
sets  of  dose  factors  are  listed  in  Table  26. 


I.  Ground  Deposition 


To  determine  the  effect  of  radionuclide  deposition  resulting  in 

groxmd  contamination,  dose  rates  in  terms  of  activity  per  unit  area 

must  be  found.  The  dose  rates  were  calculated  in  the  manner  shown  in 

Reference  32  for  a height  of  one  meter.  At  this  height,  ji  radiation 

is  not  significant,  and  since  a short  term  accident  has  been  postulated, 

all  contamination  is  assumed  to  be  on  the  surface  of  the  ground.  The 

ground  deposition  dose  factors  were  found  using  the  equation, 

n 

DF  = .869  V A,  DR.  P. 

111 

i=l 

where  A^  is  the  fractional  abundance  of  photon  i in  the  nuclide  under 

consideration,  DR^  is  the  exposxire  rate  at  one  meter  above  an  infinitely 

2 

smooth  plane  contaminated  to  one  pCi/m  , and  is  the  fraction  of 

surface  dose  which  pene'^rates  to  whole  body  depth  (5cm)  and  where  .869  = 

rads  in  air  1 rem  . -i-w  - o..-  t 

TT-: : X 3—  for  Y . An  example  is  shown  below  in  Section  J. 

R in  air  rad 

ITie  above  equation  is  modified  by  . 5 in  Reference  31  to  account  for 
ground  roughness  and/or  heavy  clothing.  The  dose  factors  listed 
therefore  are  half  the  values  shown  in  Reference  2.  It  is  the  unmodi- 
fied value  which  is  shown  in  the  following  calculation.  The  WASH  1400 
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values  m\ist  be  modified  by  1/3  to  account  for  the  shielding  effect  of 

houses,  office  buildings,  cars,  etc.,  and  even  this  value  is  conserva- 

/ 

tive  since  a more  realistic  value  would  be  less  than  this  since  more 

people  are  in  protected  or  partially  protected  environments  at  a 

given  time.  Using  1/3  accounts  for  the  fact  that  many  people  will 

leave  buildings  and  enter  cars  during  an  evacuation,  and  is  more 

appropriate  in  this  case  than  ^ or  Skin  dose  calculations  were  done 

-3 

in  the  same  manner  as  whole  body  using  a depth  of  7 x 10  cm  since 


is  affected  by  e 


(depth) 


With  a decay  energy  of  .66  MeV,  93%  of  the  time,  and  assuming  an 

a'* 

infinite,  thin  area  source,  -fov  C» 


. - V r , , -12  Cl  / Im  1 3.7 

- — — f""  ^ = 1 >■  1°  — ( Tors;)  V 

cm  sec  m ^ ' ' 


Ci-sec 


= 3.7  X 10 


-6  Y 


2 

cm  sec 


Hie  dose  is  calculated  at  a height  of  1 meter  ( ''•  3 ft)  in  air: 


jkS  ( /o  )(height)  = .077—  » .00129S\  x 

^ 1°  ’ era  ^3 

<•66  ' Ej(.099)  = 5. 


100  cm  = .099 


3.7  X lO"®  — I MS)  = 9.25  X lo"®  — j 

cm  sec  cm  : 


• *.66  ^.33  = 

cm  sec 


2 

cm  sec 


Dose  Rate  = 6.58  x lo'®  . I = 6.58  x 10~®»  ( .031)x  7.33  x 10~® 

^Jair  o 


= 1.52  X lO"^^  R/hr  = 1.52  x lo”®  MR/hr 


Dose  Factor  = .869  S A.  DR.P.  P,  = = e = e = 

^ 1 1 t i I 


~ 93%  (only  1 mode  in  decay  energy) 
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8.2  Rem/hr 


Ci/rn^ 

This  calculation  is  within  2%  of  the  value  obtained  in  Column  4 of 

Table  21.  Also  shown,  there  are  the  average  gamma  energies,  the  skin 

dose  factors  (for  reference  although  no  such  calculations  are  done  in 

this  study),  and  the  nuclide  solubility  in  cold  water.  As  stated 

previously,  the  modified  whole  body  dose  factors  are  used  and  are  one 

third  the  value  of  the  original  factors. 

To  obtain  the  dose  rate  DR  due  to  ground  deposition  in  rem/day, 

DR  = (Activity  density  in  i^)(^^^^10^^pCi)(Dose  Factor  in 

m Ci/m 

X (24  hr/day). 

To  find  the  total  dose  TD  over  an  integrated  time  from  the  beginning 


of  exposure  at  time  t=o  to  some  time  t then, 

ft  . DR  ^ 

DR  e‘  ''Mt  = (l-e~  ) 


uhci-PA  DR  is  dnse  T'ate  al-  in  T«eTn/dav  ealcnl  atad  above  and  X is 


PART  IV 

THE  MODAIREM  COMPUTER  CODE 


Figures  3 and  4 through  4.7  show  the  basic  flow  of  control 
between  subroutines  and  a simplified  flow  diagram  for  each  subroutine. 


t 


I 


PRINT 

MAN-REM 

TOTALS 


I WHOLE  BODY 

COMPUTE  MAN-REM  TOTAL 
FOR  EACH  NUCLIDE 


EVACUATION  v YES 

Vdesiredx 


PRINT  CLOUD  TIMES, 
POPULATION  FRACTIONS 
REMAINING  AT  CLOUD 
ARRIVAL  AND  AFTER 
"TIME"  DAYS 


MODIFY  POPULATION 
DUE  TO  EVACUATION 


SUBROUTINE  NOT  DESIRED 


NO,  OR  EVACUATION  SATISFIED 


FIGURE  4.2 


COMPUTE  WATER  TRAVEL  TIME 


COMPUTE  WATER  TRAVEL  TIME 


PRINT  WATER  TRAVEL  TIME 


COMPUTE  RIVER  CROSS  SECTIONAL  AREA 
USING  SIMPSON’S  RULE 


A 


y PROFILE  y 
X COUNTER  > 
S^ATISFIE^ 

YES 


PRINT  WATER 

COMPUTE  AVERAGE  RIVER  FLOW,  AVERAGE 

CONCENTRATIONS 

CONCENTRATION  OF  NUCLIDE  IH  WATER 

NO  PLANT 


/TYPE  X 

LANT  / WATER  X TYPE  4 

Xtreatment/ 

TYPE  1 Xp'amtX  3 

[ N/  1 

[ T type  2 f 

COMPUTE  NUCLIDE  CONCENTRATION 
AFTER  WATER  TREATMENT  FOR  THE 
TYPE  TREATMENT  PLANT 


y/NUMBER\ 
'^F  REACHES^ 
^COMPLETED' 


'^NUCLIDES  X NO 
SATISFIED/” 


1 


NO  PLANT 


TYPE  1 


FIGURE  4,4 

G) 


Xtype  \ 

WATER  \ TYPE  4 
TREATMENT  / 

V PLANT  X type  3 


f f TYPE  2 Y f 

COMPUTE  NUCLIDE  CONCENTRATION 
AFTER  WATER  TREATMENT  FOR  THE 
TYPE  TREATMENT  PLANT 


IF  APPLICABLE 

COMPUTE  TOTAL  CONCENTRATION  OF 
NUCLIDES  IN  TREATED  WATER 


I PRINT  J 

concentrations] 


RETURN 


SUBROUTINE  FISHL 


READ  CONCENTRATION  FACTORS 
COMPUTE  CONCENTRATION  BY  NUCLIDE 
IN  FISH  IN  LAKE 

COMPUTE  CONCENTRATION  TOTAL  IN  FISH 


RETURN 


, A.  Definitions 

I 

j 

■ There  are  many  new  variables  used  in  addition  to  the  original 
AIREM  program.  For  purposes  of  clarity,  definitions  of  the  more  im- 
portant terms  are  given  in  Tables  7 through  11. 

B.  Initial  Data 

The  initial  data  used  was  entered  according  to  Reference  10  and 
the  definitions  aur'e  as  listed. 


Number  of  sectors 

Number  of  stability  classes 

Number  of  radii 

Number  of  isotopes 

Effective  release  height 

SI6MAX 

In-plant  holdvq)  time 
Rainfall  fraction 
Washout  factor 

Wind  frequencies  by  stability  class  (%) 
Wind  speeds  by  stability  class 
Annular  ring  population  (yr  1990) 
Midpoint,  lower,  and  v^per  radii 
Cloud  dose  model 
Source  strengths 


16 

6 

12 

19+19+7  =45 
31  m 


13 


742  m 
0 

0 or  1 
.0002 
Table  2 
Table  3 
Table  4 
Table  5 
semi -infinite 
Table  6 


13 


10 
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C,  Additional  Data  Required 


I 


To  obtain  whole  body  or  skin  doses  due  to  deposition  of  radionuclide 

% contamination  on  the  ground,  and  to  use  one  of  the  two  possible  modes 

of  evacuation,  41  data  cards  are  necessary  for  20  isotopes.  Ground 

deposition  dose  factors  for  whole  body  or  skin  for  each  nuclide  are 

2 -1 

required  in  units  of  rem/hr  per  Ci/m  , as  are  decay  factors  in  days 

1 

In  addition,  the  evacuation  parameters  must  be  given.  Further  detail 
is  provided  in  Figure  5. 

Data  required  for  subroutine  RIVCON  will  vary  considerably 
according  to  the  detailed  data  available  for  the  river  under  study. 

Needed  are  ebb  and  flow  velocities  and  durations , reach  data  which 
profiles  the  river  at  various  locations , nuclide  solubilities  in  cold 
water  and  dissolved  nuclide  fractions  which  can  pass  through  the  various 
type  water  treatment  facilities  available . Information  as  to  the  river 
segments  (reaches)  and  their  relative  location  to  the  16  compass 
sector,  12  radii  map  is  also  needed.  Specifics  are  provided  in  Figxire  6. 

Following  RIVCON  is  subroutine  FISHR  which  requires  one  data 
card  per  isotope  to  pinavide  nuclide  concentration  factors  for  fresh 
water  fish.  Subroutine  LAKEDS  requires  24  data  cards  for  20  isotopes 
and  the  entering  of  swimming  and/or  boating  dose  factors,  lake 

surface  areas  and  volumes , its  mesh  point  location  and  water  treatment 

• ! 

facility  type  used.  Details  on  data  card  production  and  their  order 

8U?e  given  in  Figure  7.  Subroutine  FISHL  uses  the  data  entered  in 
FISHR  for  computation  so  no  additional  cards  are  necessary.  Siibroutine 
MILK  requires  local  or  average  data  on  crop  yields,  drinking  water 
sources  if  available,  forage  crop  concentration  factors  in  pCi/kg  of 
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fresh  plant  per  pCi/kg  of  dry  soil  for  each  nuclide,  and  the  coefficient 

of  transfer  of  each  nuclide  from  the  diet  to  milk  in  pCi/1  per  pCi/day. 

( 

Both  a biological  half-life  on  the  crop  and  in  the  cow  mxist  be  given. 
Figxire  8 lists  the  pertinent  data  required.  Subroutine  REACH  requires 
no  additional  data. 


ti 

.'f! 

t.| 


li 
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ill 
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FIGURE  5 

DATA  REQUIRED-ADDITION  TO  MAIN  PROGRAM 


1 •• 


After  five  blank  cards  in  AIREM  (semi-infinite  cloud  only  used): 


Card  Sequence  Columns 
L Cards^ 


1 Card 


Title 


Format 


1-10 

Ground  deposition  dose  factors 
( Rem/hr/ Ci/m 2)  for  WB  or  skin  for  the 

nuclide . ^ 

F10.4 

11-25 

Decay  factor  in  days  for  the  nuclide. 

E15.3 

26-30 

Integrated  time  of  exposure  to  contaminatd 

ground  (days). 

F5.1 

31-35 

Type  Dose,  0=  whole  body,  1 = skin 

15 

1-10 

Fraction  of  population  unaffected  by 

evacuation  ^ 

F10.2 

11-20 

Evacuation  rate  (days  ■'■) 

F10.2 

21-30 

Time  lag  between  interpretation  of  data 

and  order  to  evacuate 

F10.2 

31-40 

Time  between  awareness  of  core  melt  and 

leakage . 

F10.2 

41-45 

Evacuation  mode 

. -1  = Evacuation  outside  furthest 
sector. 

0 = no  evacuation. 

1 = Evacuation  outside  midway  radius 

being  considered. 

15 

If  evacuation  = 1 or  -1  repeat  L cards. 


FIGURE  6 


Card  Sequence 
1 Card 

1 Card 


L Cards 


IIRECH  Cards 


IIRECH  X 
NPROF  cards 

IIRECH  cards 


DATA  REQUIRED  SUBROUTINE  RIVCON 


Columns 

Title 

Format 

1-10 

Delay  time  (days)  before  computing  water 
concentrations.  IF  = 0.,  then  subroutine 

RIVDOS  is  not  used. 

FIO.O 

1-5 

Time  (hrs)  the  tide  flows  out. 

F5.1 

6-10 

Time  (hrs)  the  tide  flows  in. 

F5.1 

11-15 

Water  velocity  (mi /hr)  with  the  tide 

flowing  out. 

F5.1 

16-20 

Water  velocity  (mi/hr)  with  the  tide 

flowing  in. 

F5.1 

21-25 

Tide  direction  at  the  time  of  accident. 

0 = out 

1 = in 

15 

26-30 

Nimiber  of  reaches  (19  Max).  - 

15 

31-35 

Number  of  depths/cross  sectional 
profile  at  equidistant  intervals  (3,  5,  7, 

9 or  11  only). 

15 

36-40 

Number  of  profiles/reach  (20  Max). 

15 

1-15 

Plant  discharge  (pCi/sec) 

E15.6 

16-25 

Nuclide  solubility  in  cold  water. 

0.  = insoli±)le 
.5  = slightly  soluble  or 
decomposes 

1.0  = soliiile  or  unknown 

F10.2 

26  -35 

Dissolved  nuclide  fraction  passing 

water  treatment  facilities  of  Type  3 or  4 

F10.2 

1-10 

Length  of  reach  ^miles) 

F10.3 

11-25 

Surface  area  (ft  ) of  river  in  the  reach. 

E15.6 

26-30 

Water  treatment  facility  status 

0 = no  facility 

1 = class  1 facility 

2 = class  2 facility 

3 = class  3 facility 

4 = class  4 facility 

15 

1-6 

River  Width  (ft). 

F6.1 

7-12 

Equidistant  river  depths  (ft) 

6 columns  each,  3,  5,  7,  9 or  11  depths. 

F6.1 

DO  ALL 

PROFILES  FOR  1 REACH,  THEN  NEXT  REACH. 

1-5 

Reach  number 

6-10 

MM  mesh  point  location 

1-15 

II  of  the  reach 

15 

40 


FIGURE  7 


DATA  REQUIRED-SUBROUTINE  LAKEDS 


Ceu?d 

Sequence 

Columns 

Title 

Format 

1 

Card 

1-5 

Setup  indicator  number  of  LAKES 

1 = run  LAKEDS  x times 

0 - skip  LAKEDS 

15 

(IF  0 - no  other  cards  required) 

FlO.l 

1 

Card 

1-15 

2 

Lake  surface  area  (ft  ) 

E15.7 

16-25 

Time  spent  swimming  (hrs) 

FlO.l 

26-35 

Time  spent  boating  (hrs) 

FlO.l 

36-50 

Lake  volume  (gallons) 

E15.7 

1 

Card 

1-5 

Desired  doses  0 = whole  body  only 

1 = skin  dose  only 

2 = WB  and  skin  doses 

15 

6-13 

Name  of  lake 

2A4 

L Cards 

1-15 

Whole  body  dose  factor  due  to  swimming 

in  the  lake. 

E15.7 

16-30 

Skin  dose  factor  due  to  swimming  in  the 

lake. 

E15.7 

1 

Card 

1-5 

Radius  of  majority  of  the  lake 

15 

6-10 

Sector  of  majority  of  the  lake 

15 

11-15 

Water  treatment  facility  used 

(same  as  in  RIVCON). 

15 

IF  MORE  THAN  ONE  LAKE  IS  USED  ENTER  1 LAKEDS  DATA  SET,  THAN  CORRESPONDING  FISHL 
DATA,  IHEN  NEXT  LAKEDS  DATA  SET. 


FIGURE  8 

DATA  REQUIRED-SUBROUTINE  MILK 


Card  Sequence  Columns  Title  Format 


1 1-5  Number  of  mesh  points  not  using  wells 

-1  = skip  subroutine  15 

1 1-10  Forage  ingestion  by  cow  (Kg/day)  F10.2 

11-20  Deposition  retension  factor  F10.2 

21-30  Forage  crop  yield  (Kg/m^)  F10.2 

1 1-10  Ingestion  of  drinking  water  by  dairy 

cow  (L/day)  - F10.2 

11-20  lypical  drinking  trough  size  (ft  ) F10.2 

21-30  Typical  trough  capacity  (gal)  F10.2 

NUMBER  1-5  Drinking  water  type 

1 = river 

2 = lake  15 

6-10  Radius  of  drinking  water  source  15 

11-15  Sector  of  drinking  water  source  15 

16-20  Water  treatment  facility  type  15 

1 1-5  Maximum  number  of  days  for  which  milk 

calculations  are  to  be  calculated  15 

6-10  Increments  (days)  between  successive 

concentration  calculations  15 

1 1-5  Evacuation  mode/use  of  uncontaminated 

* 4? 

fes!<?.:  0 = no  evacuation 

1 = evacuation  after  1 day  15 

6-15  Biological  half-life  (days)  of  nuclide 

on  forage  crop  F10.2 

L 1-10  Forage  crop  concentration  factor  (B) 

pCi/Kg  fresh  plant 

pCi/Kg  dry  soil  E10.2 

11-15  Coefficient  of  transfer  S , from  diet  to 

milk  (pCi/1  per  pCi/day)  E15.1 

1 1-10  Nuclide  half  life  in  the  dairy  cow  F10.2 
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Subroutine  RIVCON 


A partial  listing  of  the  initial  data  used  in  this  subroutine  is: 

Delay  time  1 day 

Time  of  ebb  flow  (hrs)  6.0 

Time  of  flow  tide  (hrs)  5.9 

Average  water  velocity  at  ebb  tide  (mi/hr)  2.1 

Average  water  velocity  at  flow  tide  (mi/hr)  1.5 

Tide  direction  at  accident  time  FLOW 

Number  of  reaches  19 

Depths /pro file  5 

Profiles /reach  10 

Plant  discharge  (pCi/sec)  0. 

Nuclide  solubility  in  cold  water  Table  21 

Nuclide  fraction  passing  water  treatment  facilities  Table  22 

Reach  length,  surface  area,  water  treatment  facility 
type,  reach  number,  and  mesh  point  locations 
of  the  reaches  Table  23 

Stibroutine  LAKEDS 

Thirteen  lakes  are  examined  in  this  study  and  whole  body  doses 
only  cire  calculated.  These  lakes  with  their  associated  drinking  water 
facility  data  is  given  in  Table  24.  Water  submersion  dose  factors 
used  are  listed  in  Table  25. 

Subroutine  FISHR  and  FISHL 

The  values  used  in  subroutines  FISHR  and  FISHL  for  fresh  water 


fish  concentration  factors  ar*e  listed  in  Table  31. 


Subroutine  MILK 


In  this  study,  18  mesh  points  using  other  than  underground  well 

water  were  examined.  The  cow  is  assumed  to  ingest  65  kg  of  forage  per 

day,  and  the  fraction  of  deposition  retained  by  the  forage  crop  is 

2 

taken  as  .25.  Also  the  forage  crop  yield  is  a^jsumed  1.8  kg/m  . Milk 
concentrations  are  computed  generally  for  60  days  in  increments  of  10 
days  and  generally  for  the  case  of  no  evacuation.  The  biological  half- 
life  of  the  forage  crop  on  the  ground  is  14  days  and  in  the  cow,  it  is 
assumed  to  be  two  days  (an  overly  conservative  estimate  for  the  cov; 
of  .5  days  yields  concentrations  of  10-20%  higher  than  for  two  days). 
Table  32  lists  the  concentration  factors  (B)  of  various  isotopes  on 
the  forage  crop  and  the  coefficients  of  transfer  of  the  various  nuclides 
from  the  cow's  diet  to  the  milk. 


PART  V 


RESULTS 


Typical  values  of  radioactivity  are  1.873x10  Curies  deposited  in 

490  square  miles  (50  mile  radius),  specifically  in  one  wind  sector  of 

I 

22^5  degrees  and  31  square  miles  which  represents  approximately  20%  of 
the  total  Curies  released  in  the  accident. 


A.  Water  Concentrations  and  Resultant  Doses 


Subroutines  RIVCON  and  LAKEDS  give  nuclide  concentrations  for  the 
various  nuclides  in  pCi/1  before  and  after  watei*  treatment.  The  calcu- 


lations below  are  for  treated  water  in  Alcove  Reservoir  after  northerly 

winds,  however,  the  same  methods  can  be  used  in  any  of  the  affected 

locations.  The  computed  nuclide  concentrations  are  listed  in  Column  2 

of  Table  27.  These  concentrations  are  converted  to  approximate  doses 

26 

with  the  use  of  the  equation. 


„ AfETe  r 

D = ^ (1  - expC 

m 


1. 26x10*^ 1.6x10  ® X 3.2x10^^)  rem 
T .693x10^  '' 


where  D = the  total  dose;  f = the  fractional  uptake  to  the  organ  of 

interest  by  ingestion;  E=  the  effective  decay  energy;  T^  = the  effective 

half-time  of  the  nuclide  in  the  organ  of  interest;  m = the  mass  of  the 

^6  15 

organ  of  interest  in  grams;  1.6x10  = erg/Mev;  3.2x10  = disintegrations/ 

2 

day  per  Ci;  .693x10  = erg/gm/rad;  A = 1 to  calculate  rem  per  Curie 


ingested.  For  example,  using  the  values  of  E,  T^,  m,  and  f from  ICRP-II 

the  dose  to  the  thyroid  due  to  the  ingestion  of  one  liter  of  contaminated 
131 

water  due  to  I is, 

„ _fl  rem.>  , oo^/''T  ,--1. 26x10^^, 1.6x10  ®x3.2xl0^^\ 

° 131  (.3)(.23)(7.6)(l-exp(  ) 

THY 

= 1.94x10®  rem/Ci 


r 


si 

‘i 


k 


!• 

(• 

i 


I 

\ 


1 


! 
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This  is  done  for  each  nuclide  for  each  organ  of  interest,  and  multiplying 
these  dose  factors  times  the  concentrations  yields  the  values  in 
Columns  6 and  7 of  Table  27.  Although  the  results  are  displayed  only 
for  Alcove  Reservoir,  the  same  results  can  be  found  from  the  data  for 
Basic  Reservoir,  lakes  near  Kingston,  Saugerties,  Cooper,  Staatsburg, 

Hollister,  Churchtown,  Coxsackie,  Potuck  Lake,  Onteora  Pond,  Lake  Forest 

eUid  Ashokan  Reservoir  when  the  winds  blow  in  the  appropriate  direction.  i 

B.  Radionuclide  Concentrations  in  Fish  1 

] 

j 

The  individual  and  total  radionuclide  inventory  ingested  by  an 

i 

average  fish,  both  in  the  various  river  locations  and  in  the  major 

lakes  was  calculated  using  the  concentration  factors  as  previously  ^ 

listed.  Total  concentrations  as  high  as  3.1  Ci/Kg  were  recorded.  It  ' 

is  probable  that  this  level  of  radioactive  material  will  not  actually  j 

i 

be  ingested  for  two  reasons . First , there  could  be  rapid  death  of 

the  fish  depending  on  the  nuclides  ingested  and  the  resultant  total 

dose  received,  and  second,  the  assumption  of  one  fixed  concentration 

factor  will  not  be  valid  when  extrapolated  to  these  large  values.  The  \ 

uptake  of  these  nuclides  will  vary  in  response  to  environmental  changes 

and  in  response  to  previously  ingested  levels. 

C.  Milk  Concentrations 

The  major  assumption  in  producing  the  results  shown  in  Tables  28 
and  29  is  that  the  cows  live  and  produce  the  milk  from  which  the  concen- 
trations are  determined.  Because  the  milk  concentration  equations  are  ^ 

semi-empirical,  they  provide  only  an  estimate  of  concentration  at  each 

1 
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period  of  time,  although  the  values  for  I,  Sr,  and  Sr  fit  other 
specific  equation  reasonably  well.  The  initial  concentrations  for  one 
day  are  overly  conservative  but  agree  closely  with  observed  data  in 
Reference  22  for  the  above  three  nuclides  by  11  days.  Comparison  of 
the  two  tables  shows  the  value  of  evacuation  of  the  dairy  cows,  or 
more  practically,  changing  to  xincontaminated  feed.  Converting  from 
concentrations  to  doses  due  to  ingestion  of  one  liter  of  milk  are 
accomplished  as  for  water. 

D.  Lying  Inhalation  Dose 

Lung  doses  to  an  individual  due  to  inhalation  range  from  22,000  rem 
at  mile  to  39  rem  at  miles  in  the  south-southwesterly  direction 
from  25,200  rem  to  47  rem  to  the  west,  and  17,100  to  36  rem  to  the 
north.  Relative  doses  between  sectors  are  as  with  others  reported  in 
this  study  and  can  be  seen  in  Table  16.  The  man-rem  totals  to  the  lung 
are  shown  in  Table  17.  The  nuclides  which  contribute  most  to  the 
man-rem  dosage  are: 


106^ 

Ru 

18% 

^°3ru 

15% 

132^e 

14% 

131i 

13% 

10% 

99 

Mo 

8% 

133j 

6% 

These  nuclides  con^irise  84%  of  the  total  dose  integrated  over  30  days. 
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E.  Gastrointestinal  Tract  Dose  Due  To  Inhalation 


Gastrointestinal  doses  to  an  individual  due  to  inhalation  vary 
from  48,600  rem  at  one-half  mile  in  the  south-southwesterly  direction 
to  32  rem  at  47^5  miles.  To  the  west,  the  doses  run  from  48,600  rem 
at  one-half  mile  to  36  rem  at  47^5  miles.  To  the  north  15,900  rem  to 
31  rem  are  calculated,  and  in  the  direction  of  least  consequence , 
southwest,  from  14,600  to  29  rem  doses «re  found  for  each  person. 
Specifics  can  be  found  in  Table  14.  The  gastrointestinal  man-rem 
totals  are  displayed  for  each  direction  in  Table  15.  Five  nuclides 
contribute  76%  of  the  total  GI  dose  with: 


132^e 

30% 

^^Mo 

22% 

103^ 

Ru 

10% 

Ru 

9% 

140„ 

Ba 

5% 

These  doses  are  integrated  over  a 50  year  period. 

F.  Thyroid  Inhalation  Dose 

Thyroid  doses  to  an  individual  due  to  inhalation  of  the  nuclides 

in  this  study  vary  from  300,000  rem  at  h mile  south  southwesterly  to 

4100  rem  at  15  miles  to  660  rem  at  47^5  miles.  In  the  westerly  direction 

where  the  highest  individual  doses  occur,  the  mile  dose  is  710,000  rem, 

4000  rem  at  22%  miles  and  1000  rem  at  47%  miles.  The  smallest  doses 

are  in  the  south-southeasterly  direction  with  the  doses  in  the 

northerly  direction  towards  Albany  slightly  larger  as  can  be  seen  in 

Table  12.  The  thyroid  man-rem  doses  are  listed  in  Tcdjle  13.  In  this 

131 

case,  iodine  contributes  97%  of  the  dose  to  the  thyroid  with  I 


f 


T 


133  135 

contributing  60%,  I,  30%,  and  I contributing  7%  of  the  dose.  The 
integration  time  used  is  assumed  essentially  infinite  since  a 50  year 
integration  time  is  used. 

G.  Whole  Body  Inhalation  Dose 

Individual  whole  body  doses  due  to  inhalation  of  45  nuclides  in 
this  study  vary  from  1286  rem  in  the  south-southwesterly  direction 
where  the  maximum  rem  dose  occurs  at  one-half  mile  from  the  site, 
dropping  sharply  to  100  rem  at  MSg  miles  and  to  16  rem  at  15  miles.  The 
highest  dosage  to  an  individual  occurs  in  the  westerly  direction  from 
a maximum  of  3123  rem  at  one-half  mile  to  90  rem  at  miles  to  13  rem 
at  22^  miles.  In  the  northerly  direction  the  dose  is  990  rem  at  \ mile, 
13  rem  at  15  miles,  and  approximately  2 rem  at  Albany.  (See  Figures  9 
and  10.)  Approximately  three-quarters  of  the  dose  can  be  attributed  to 
the  following  nuclides : 


131j. 

25% 

132^^ 

15% 

133j 

11% 

^"°Ba 

10% 

®®Mo 

7% 

5% 

H.  Whole  Body  Cloud  Submersion  Dose 

Whole  body  doses  to  an  individual  due  to  radiation  from  a passing 
nuclide  cloud  vary  from  760  rem  at  h mile  to  .61  rem  at  473$  miles  in 
the  south-southwesterly  direction,  from  1730  to  .57  rem  to  the  west 
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MAN -REM  INTEGRATED  DOSE -WHOLE  BODY  DUE  TO 
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FIGURE  9 

jPOPULATION  DOSE  vs.  COMPASS  SECTOR 
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FIGURE  10 
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and  580  to  .62  retn  to  the  north.  Doses  in  selected  directions  can 

be  seen  in  Table  18.  Man-rem  totals  can  be  compared  in  Table  19.  The 

135 

iodines  contribute  75%  of  the  man-rem  dose  in  this  mode  with  I 
contributing  30%,  at  19%,  at  12%, at  9%,  and  at 


I.  Whole  Body  Ground  Deposition  Dose 


Whole  body  doses  due  to  standing  on  contaminated  ground  for  one 
day  immediately  following  the  reactor  accident  are  calculated  and 
displayed  in  Figure  11.  The  man-rem  dosage  calculated  shows  the  same 
general  shape  as  Figure  1 for  whole  body  inhalation.  The  variation 
in  individual  dose  can  be  seen  in  Table  20  which  shows  the  factor  of 
1000  decreases  in  dose  at  approximately  50  miles  from  initial  ^ mile 
dose  values,  and  the  man-rem  dosage  by  direction  are  listed  in 
Table  20A. 

Five  nuclides  contribute  88%  of  the  total  ground  deposition  dose: 

^°^Ru  40% 

^°®Ru  15% 

17% 

^^“^Cs  12% 

4% 


J.  The  Effect  of  Meandering  Winds  on  Dose 


The  iodine-131  whole  body  dose  is  shown  in  Figure  12.  The  second 
column  shows  the  dose  due  to  winds  blowing  into  the  north.  The  third 
column  shows  the  doses  due  to  the  winds  blowing  into  the  two  directions 
shown.  The  doses  may  be  higher  or  lower  than  the  average  based  on  the 
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Figure  12 


Iodine  131  Dose  (Rem)  vs.  Wind  Variation 
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weighted  average  of  the  fractional  length  of  time  the  wind  blows  in 
each  sector.  The  fourth  column  is  for  winds  blowing  from  north  to 
east  with  all  wind  sectors  having  equal  probabilities  of  occurrence 
and  all  stability  classes  having  the  same  probability  of  occurrence. 
This  figure  shows  the  decreasing  dose  with  increased  meandering  of 
the  wind  since  slightly  more  radioactivity  is  spread  over  much  larger 
areas.  This  does  not  mean,  however,  that  there  will  necessarily  be 
a man-rem  reduction  based  on  the  population  distribution  encountered. 

K.  Swimming  Doses 

As  with  the  lake  doses,  the  swimming  doses  are  only  discussed 
for  Alcove  Reservoir  but  apply  to  eill  lakes  in  the  area  of  interest. 

For  a person  who  happens  to  swim  in  the  contaminated  lake,  for 
periods  shortly  after  the  deposition,  the  total  dose  is  .47  rem/hr  and 
will  diminish  as  the  nuclides  decay.  Although  the  percentage  contri- 
bution of  the  nuclides  is  not  the  same  in  general  the  same  nuclides 
contribute  to  the  total  dose  as  to  the  water  ingestion  dose.  Actually, 
it  should  be  noted  that  swimming  and  boating  are  not  permitted  in 
Alcove  Reservoir. 


L.  Evacuation  Model 

I 

I 

Using  the  evacuation  parameters  of  evacuation  rate  A=  8.30  days”^ 
(1%  = two  hours),  the  fraction  not  evacuated  A = .10,  the  time  between 
the  aweu?eness  of  core  meltage  and  actual  leakage  T = .02  days  (%  hour), 
the  time  needed  to  interpret  the  data  and  issue  the  warning  TL  = .02 
days,  and  the  average  wind  speeds,  the  population  remaining  after  passage 
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of  the  cloud  range  from  96-99%  in  the  nearest  radius  to  from  19-35% 
in  the  outermost  sector  when  evacuation  is  carried  out  to  a full  50 
miles  downwind  from  the  accident  site.  It  is  interesting  to  note 
that  the  direction  of  highest  dose  per  individual,  west,  is  one  of  the 
directions  allowing  the  greatest  time  for  evacuation  since  the  sub- 
mersion and  inhalation  doses  at  each  location  depend  on  the  stay  time 
of  the  cloud.  The  westerly  direction  allows  from  4%  to  88%  of  the 
population  (according  to  their  distance  from  the  accident  site)  to 
evacuate  while  the  direction  of  smallest  dose  per  individual  south- 
sovrtheast  allows  1%  at  ^ mile  to  65%  at  50  miles.  In  the  Troy-Albany 
area,  approximately  60%  of  the  people  could  be  evacuated  under  average 
conditions.  The  sector  of  highest  man-rem  dosage  south-southwest 
allows  2-73%  evacuation. 

With  the  above  evacuation,  man-rem  dosage  can  be  significantly 
reduced.  To  the  north,  the  reduction  in  man-rem  dosage  is  found  to  be 
56%;  to  the  south-southeast,  a 33%  reduction;  to  the  south-southwest, 
40%;  and  to  the  west,  39%;  and  the  overall  reduction  in  dose  is 
found  to  be  approximately  40%. 

Evacuation  to  10  miles  in  each  accident  case  reduces  the  man-rem 
dosage  only  approximately  8%.  Because  of  the  average  cloud  speeds, 
evacuation  to  10  miles  would  not  remove  enough  people  to  satis- 
factorily reduce  the  total  dosage . The  highest  winds  recorded  in 
the  area  were  71  miles  per  hour  from  the  northwest.  Using  this 
representative  wind  speed,  the  individual  doses  are  reduced  4 to  6 
times  as  are  the  amounts  of  deposition.  When  rainfall  is  considered 
to  occur  during  the  accident  period,  no  significant  variations  are 
seen  in  the  dose  to  an  individual.  Since  more  of  the  nuclides  are 
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washed  out,  the  ground  deposition  dose  reaches  two  to  three  times  its 
reference  value  with  dry  weather,  but  this  is  still  only  a small  per- 
centage of  the  total  dose  received. 

If  the  evacuation  rate  is  taken  as  4.15  days  ^ corresponding  to 
an  evacuation  "half-life"  of  four  hours,  or  one-half  the  reference 
rate,  evacuation  in  the  path  of  westerly  winds  allows  only  2%  to  escape 
at  *5  mile  while  78%  can  evacuate  from  50  miles  away  (10%  less  than  with 
a two  hour  evacuation  "half-life").  To  the  south-southeast,  from  1% 
to  49%  can  evacuate  and  in  the  Troy-Albany  area,  approximately  one- 
half  worild  be  able  to  leave  in  time,  resulting  in  a 45%  reduction  in 
dose  (compared  with  a 68%  reduction  with  the  quicker  evacuation  rate). 
Overall,  there  is  only  a 25%  reduction  in  man-rem  dosage  due  to  evacua- 
tion. With  a four  hour  evacuation  "half-life"  plus  a one  hour  delay 
to  evaluate  the  data  and  issue  the  warning,  there  will  be  no  evacuation 
within  two  miles  and  very  little  out  to  five  miles.  In  Albany,  the 
population  could  be  reduced  by  41%  (or  about  10%  less  than  before)  which 
results  in  a 31%  reduction  in  dose.  Overall  there  is  a 20%  reduction 
in  dose  in  this  evacuation  case . These  results  are  summarized  in 
Figure  13. 
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Figure  13 
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PART  VI 


k 


DISCUSSION  AJJD  CONCLUSIONS 

It  seems  obvious  that  a reactor  accident  will  require  a number  of 
ejq)erts  to  outline  contaminated  areas  and  insure  area  and  water  moni- 
toring. The  greater  the  severity  of  the  release,  the  more  experts  will 
become  involved.  With  northerly  winds,  the  water  in  Alcove  Reservoir 
will  become  contaminated  requiring  alternate  water  supplies.  The  values 
listed  in  Table  27  do  not  allow  for  decay  in  transmission  but  are  high 
enough  that  contamination  will  remain  for  some  time  without  special 
treatment . 

Nuclide  Concentration  In  Milk 

Estimation  of  radionuclide  concentrations  in  milk  is  a difficult 
task,  and  data  taken  from  environmental  samples  is  not  easily  repre- 
sented by  a single  curve  even  for  each  isotope  as  evidenced  by  results 
in  the  literature  (1,  7,  14,  15,  23).  As  a result,  the  concentrations 
of  the  various  nuclides  although  reasonable  figures , are  by  no  means 
exact.  They  do  show,  however,  that  the  concentrations  will  be  extremely 
high  so  that  if  the  cows  do  live  to  produce  milk,  it  will  be  contaminated 

for  months  after  the  event,  with  high  concentrations  of  the  long-lived 
go  137 

isotopes  of  Sr  eind  Cs. 

It  can  be  seen  that  several  protective  actions  can  be  used  to  re- 
duce dose  even  in  this  worst  case.  Replacing  forage  crops  with  uncon- 
taminated feed  is  extremely  important  in  reducing  concentrations  and 

137 

doses.  For  example,  the  tenth  life  of  Cs  can  be  reduced  from  67 
da^s  to  8 days  by  placing  a cow  on  uncontaminated  feed.  As  stated, 
removal  of  the  milk  must  be  accomplished  initially  and  extensive 
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monitoring  over  a wide  area  will  be  in  order.  Taking  protective  actions 
includes  the  weighing  of  risks  of  taking  no  action  against  the  impli- 
cations of  the  action,  but  in  this  case,  there  will  be  no  reason  for 
hesitant  action  by  this  point  in  time.  Assumptions  following  contamina- 
tion are  that: 

I 

1.  The  contaminating  event  is  a single  acute  occurrence  within 
a few  hours. 

2.  Previous  accumulation  of  nuclides  is  negligible. 

3.  The  forage  crop  is  grass  which  supplies  all  of  the 

22 

roughage  consumed,  and  it  is  growing  at  its  normal  rate. 

Total  Dose  Received 

Sxjmmation  of  the  whole  body  doses  due  to  a single  radioactive 
cloud  submersion,  a one  day’s  dose  due  to  ground  contamination,  and 
a 30-day  integrated  dose  due  to  inhalation  of  various  radionuclides 
diiring  the  cloud  passage  yield  representative  totals  useful  in  comparing 
the  relative  doses.  The  highest  dose  per  individual  occurs  in  the 
western  sector  with  4846  rera  at  one-half  mile  to  3.66  rem  at  47^5  miles 
and  1604  rem  to  2.85  rem  to  the  north.  Specific  values  for  three 
directions  are  shown  in  Figure  14.  The  relative  contributions  of  the 
three  sources  of  whole  body  dose  vary  with  distance  as  shown  in  Figure  15 
with  the  inhalation  contribution  running  from  approximately  60%  at  one- 
half  mile  to  over  77%  at  4735  miles  while  the  cloud  submersion  dose  de- 
creases in  the  same  interval  from  36  to  20%.  A possible  explanation 
for  this  effect  is  that  the  cloud  being  carried  downwind  diffuses 
vertically  and  horizontally  with  the  resultant  vertical  midpoint  of  the 
cloud  mass  z rising  with  distance.  In  the  MODAIREM  Code,  equal  cloud 
deposition  from  all  parts  of  the  cloud  is  assumed.  The  receptor  lies 


Figure  14 
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FI6URE  15  RELATIVE  DOSES  vs  DISTANCE  FROM  ACCIDENT 


Figure  16 

RELATIVE  NUCLIDE  CONCENTRATIONS  TO  TOTAL  DOSE 


in  the  center  of  this  semi-infinite  unifonn  medium  and  the  cloud 


radiation  dose  at  any  distance  can  be  found  using  reciprocity  rela- 
tionships, and  it  is  essentially  a three  dimensional  effect  which 
depends  directly  on  the  passage  of  the  cloud.  The  inhalation  dose 
is  essentieilly  a point  effect  depending  on  the  concentration  of  the 
cloud  at  that  location  dxiring  passage.  That  is,  close  to  the  acci- 
dent site  where  the  cloud  is  small,  it  passes  over  the  receptor, 
and  most  of  the  cloud  contributes  to  the  submersion  (radiation)  dose 
and  likewise  for  the  inhalation  dose . Far  from  the  accident , part  of 
the  cloud  which  has  expanded  does  not  contribute  proportionally  as 
much  as  close  to  the  receptor  since  the  potential  target  (receptor) 
relative  area  is  smaller.  The  inhalation  dose,  however,  although 
smaller  due  to  expansion,  and  assuming  homogeneity,  all  contributes 
to  the  total  dose.  Continuing  this  assun^tion,  and  averaging  the 
relative  doses  yields  the  average  contributions  to  the  total  dose  which 
are  69.7%  due  to  inhalation,  28%  due  to  cloud  radiation  and  2%  due  to 
cloud  deposition.  With  these  figures,  the  relative  contributions  of 
each  isotope  to  the  total  dose  can  be  calculated,  and  are  shown  in 
Figure  16. 

Evacuation  vs.  Dose 

One  can  see  from  Figure  14  that  doses  to  the  individual  drop 

sharply  as  the  distance  from  the  site  increases.  An  acute  exposure 

2 

received  of  less  than  150R  will  result  in  no  lethalities.  It  has 
been  seen  already  that  evacuation  of  10  miles  reduced  the  man-rem 
only  slightly,  however,  its  effect  on  the  individual's  dose  is  very 
important  in  developing  an  emergency  plan.  Evacuation  to  over  7^5  miles 
from  the  site  reduces  the  dose  equivalent  to  under  130  rem  in  the  worst 
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case.  In  fact,' in  the  most  probable  wind  direction,  north,  the  total 
is  reduced  to  57  rem.  The  manner  of  addition  of  doses  described  pre- 
viously is  again  emphasized.  Evacuation  to  about  miles  would  effect 
the  towns  of  Catskill,  Saugerties  and  Tivoli. 

A word  must  be  said  concerning  any  proposed  evacuation.  In  a 
28 

study  done  it  was  found  that  successful  evacuation  can  be  accom- 
plished given  sufficient  warning  and  adequate  reaction  time  according 
to  an  inverse  relationship  between  the  population  density  and  the 
evacuation  time.  The  1990  population  of  the  area  near  Cementon  shows 

a population  density  of  between  approximately  200-600  people  per  square 

28 

mile  which  corresponds  to  four  to  six  hours  evacuation  time  necessary. 
This  agrees  with  the  evacuation  half-time  of  two  hours  used  in  the 
model.  This  appears  to  be  an  accurate  representation  and  is  extended 
by  the  model  which  shows  that  the  radioactive  cloud  takes  approximately 
one  hour  to  travel  to  ^ miles  which  results  in  few  individuals  being 
able  to  leave  the  area  in  time.  Then  the  time  to  evacuate  from  a 
distance  closer  to  the  site  is  even  less  and  successful  evacuation  is 
less  probable.  It  has  been  shown  that  evacuation  is  sensitive  to  the 
time  necessary  for  evaluation  of  the  accident  and  issuance  of  a warning. 
Due  to  the  encompassing  consequences  of  taking  such  action,  the  person 
responsible  would  make  stare  he  was  making  a correct  decision,  and  if 
the  decision  is  ultimately  positive,  every  minute  in  decision  enlarges 
the  consequences. 

Evacuation  vs . Shelter 

A special  aspect 'of  evacuation  during  a nuclear  event  is  first, 
who  shotild  evacuate  and  second,  where  should  they  go?  In  answering 
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these  questions,  one  m\jst  ask  first,  in  what  direction  are  the  winds 

which  are  carrying  the  radioactivity  travelling,  and  second,  in  what 

direction  will  they  blow  in  the  future.  Given  evacuation  routes, 

even  if  known,  may  cause  many  travellers  to  move  into  the  path  of  the 

cloud  which  would  be  otherwise  avoided  if  the  wind  should  suddenly 

shift  direction.  Most  evacuation  can  be  carried  out  in  an  orderly 

28 

fashion  without  panic,  according  to  an  EPA  study  but  there  is  poten- 
tial for  trouble  due  to  the  psychological  response  to  an  event  of  this 
type  and  the  very  short  times  required  to  leave  the  area. 

The  worst  case  results  if  a person  does  not  leave  the  area  since 
he  would  receive  the  maximxim  dose.  The  best  situation  is  an  evacuation 
away  from  the  cloud  resulting  in  no  dose.  The  majority  of  the  popula-- 
tion,  as  it  has  been  shown,  will  fall  somewhere  in  between.  An  unsuc- 
cessful attempt  at  evacuation  will  result  in  a needlessly  high  dose 
although  the  farther  away  from  the  site,  the  better  the  chances  for 
evacuation.  Remaining  in  place  and  taking  shelter  is  an  alternative 
that  should  be  considered.  For  example,  in  Cementon  there  eu:^  locations 

designated  as  fallout  shelters  which  provide  room  for  2,224  individuals 

30 

with  a protection  factor  of  at  least  40  for  gamma  radiation.  If  in 
addition,  two  diy  bath  towels  or  similar  material  are  used  as  a make- 
shift filter  in  the  absence  of  civilian  protective  masks,  it  has  been 
estimated  that  85%  of  the  inhalation  dose  can  be  avoided.  Since  inha- 
lation dose  primarily  occurs  during  cloud  passeige,  the  estimated  dose 
can  further  be  reduced.  It  takes  approximately  two  hours  to  exchange 
air  volume  in  a building,  therefore,  a dilution  can  take  place.  The 
better  sealed  a building  is,  the  less  the  inhalation  dose.  Of  course 
absolute  filters,  especially  charcoal  activated  filters  drastically 
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reduce  the  total  dose  through  a reduction  in  the  thyroid  dose.  If  the 
effect  on  inhalation  dose  is  to  lower  it  by  an  additional  25%,  then 
even  at  one  mile  no  protected  person  will  receive  over  150  rem.  Clearly, 
this  is  not  entirely  satisfactory,  but  this  is  a highly  xonlikely  event 
and  is  superior  to  an  attempt  at  evacuation  for  short  distances  from 
the  reactor  where  the  cloud's  travel  time  to  one  mile  is  approximately 
ten  minutes.  Remaining  in  place  resulting  in  a dose  from  the  nuclides 
deposited  on  the  ground  does  not  appreciably  add  to  the  dose.  The  use 
of  immediately  available  shelters  and  either  stockpiled  protective 
masks  or  makeshift  filters  is  reasonably  effective.  In  rural  areas, 
bcisements  cp.uld  be  used  with  similar  success.  As  a public  relations 
gesture  as  well  as  a protective  one,  the  construction  of  the  plant  witK 
a large  hall  equipped  with  a filtration  system  at  an  easily  accessible 
location  in  or  near  Cementon  which  could  provide  a type  of  civic  center 
and  recreational  hall.  This  scheme  has  two  additional  benefits.  First, 
it  takes  advantage  of  the  inertia  of  the  population  who  tend  to  continue 
to  take  actions  they  have  in  the  past  rather  than  to  undertake  new 
patteras.  That  is,  a movement  into  the  cement  plants,  into  a basement 
or  into  a recreational  hall  is  more  likely  to  occur  than  an  evacuation 
of  families,  even  if  time  were  available.  Second,  the  sheltering  of 
individuals  would  allow  various  disaster  mechanisms  and  health  physics 
personnel  to  arrive,  assess  the  situation  and  proceed  on  the  best  course 
of  action  with  a minimum  of  confusion.  The  warning  of  individuals  of 
the  need  for  moving  into  shelters  would  be  no  more  difficult  than 
warning  them  of  the  need  to  evacuate,  in  fact,  it  could  be  much  easier. 
If  the  local  population  is  educated  to  the  facts  and  dangers  involved 
and  simple  protective  actions  to  take,  one  day  in  a shelter  could  be 


tolerable.  A type  of  siren  or  system  of  sirens  which  could  be  heard 
for  at  least  10  miles  could  be  used  as  a warning  device  for  everyone 
in  the  area  but  woxild  require  that  the  population  know  what  it  meant. 

It  is  unfortunate,  however,  that  this  type  of  action,  clearly  the 
quickest  and  therefore  the  best  would  not  be  used  because  of  the 
accentuation  of  alarmist  tendencies  by  the  opponents  of  nuclear  power. 
Again  they  would  question  the  need  for  these  actions  if  nuclear  power 
were  so  safe,  so  in  preparing  for  an  extremely  unlikely  event,  fears 
are  raised  that  an  accident  will  take  place  because  of  preparation 
for  all  eventualities. 

Comparison  With  WASH-1400 

Using  values  for  lethal  doses  of  a five  day  to  the  gastroin- 
testinal tract  of  4,000-6,000  rad,  an  estimated  of  3,000  rad  to 

the  lung,  and  whole  body  LD^  of  225R,  of  400R,  and  of  5-600R, 

the  worst  case  to  the  west  yields  and  estimated  326  acute  deaths  and 
the  most  probable  case  to  the  north  where  the  wind  blows  one  quarter 
of  the  time,  an  estimated  12  acute  deaths  from  this  accident  assuming 
no  protective  actions  whatever.  The  wind  blows  to  the  west  less  than 
2%  of  the  time.  The  number  of  thyroid  modules  is  extremely  high  since 
there  are  over  3,000  rem  per  person  to  the  thyroid  out  to  25  miles  to 
the  west,  and  out  to  15  miles  to  the  north.  There  appears  from 
Table  12  that  nodules  will  be  fomed  in  children  for  downrange  throughout 
the  range  of  study  if  no  protective  actions  are  taken.  Since  the  ma- 
jority of  the  thyroid  dose  comes  from  inhalation,  this  is  another  reason 
to  reduce  the  inhalation  dose.  Comparison  of  these  results  with 
^pendix  VI  of  WASH- 1400  show: 


Highest  man-rem  in  this  study  (SSW)  = 2.79x10® 

Average  man-rem,  WASH  1400  = 2.8  x 10® 

Peak  acute  fatalities , this  study  = 326 

Peak  acute  fatalities,  WASH  1400  = 1,100 

Fatalities  (northerly  winds)  = 12 

Average  acute  fatalities,  WASH  1400  = 34 


Both  studies  are  to  an  extent  averages  over  populations,  areas  and 
meteorological  conditions. 

Thus , the  consequences  of  a veiy  serious  reactor  accident  have  been 
shown.  Rather  than  concentrate  on  one  single  area,  a broad  perspective 
was  taken,  and  it  is  expected  that  several  of  these  areas  for  the  Cementon 
plant  will  be  expanded  upon  and  refined, in  the  future.  Once  again,  the 
probability  of  this  accident  must  be  emphasized  as  being  very  small. 

If  the  population  close  to  the  plant  were  educated  in  basic  protective 
action  and  an  emergency  plan,  there  is  a possibility  of  preventing  acute 
fatalities.  There  would,  however,  be  many  long-term  effects  that  this 
study  has  not  considered.  Putting  the  benefits  of  the  plant  in  perspective 
with  the  risks  involved,  this  study  must  support  the  location  at  Cementon 
as  being  adequate  in  this  respect.  It  is  qualitatively  safer  than  the 
case  studied  in  WASH  1400,  and  the  worst  effects  can  be  considerably 
mitigated. 
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Table  1 


DOSE  INHALATION  FACTORS  (mrem/sec/Ci/m^) 
(Rem/CiXBreathing  rate  of  2.2xlo“**  m^/sec)(1000  mrem/rem) 


Isotope  VfB  (30  day)  Lung  (30  day)  Thyroid  (50  yr)  GI  tract  (50  yr) 


KR  85 

0 . 

0 

0 

0 

KR  85* 

0 

0 

0 

0 

KR  87 

0 

0 

0 

0 

KR  88 

0 

0 

0 

0 

SR  89 

858 

17160 

0 

9900 

SR  90 

352 

41360 

0 

7920 

SR  91 

77 

1210 

0 

5280 

Y 90 

132 

4620 

0 

15620 

Y 91 

572 

18920 

0 

9460 

Zr  95 

1078 

16940 

0 

4840 

Zr  97 

86 

2200 

0 

15620 

Nb  95 

440 

7480 

0 

2420 

Mo  99 

145 

2640 

0 

6600 

Tc-  99* 

2 

17.8 

0 

99 

Ru  103 

170 

8140 

0 

4840 

Ru  105 

14 

352 

0 

2420 

Ru  106 

627 

S1480 

0 

23540 

Rh  105 

22 

572 

0 

2420 

Te  129 

4 

75 

8.8 

462 

Te  129* 

968 

23540 

1474 

12100 

Te  131* 

161 

2420 

220 

7920 

Te  132 

431 

6820 

550 

12100 

I 131 

528 

4620 

323400 

242 

01  131 

528 

4620 

323400 

242 

I 132 

29 

198 

11660 

528 

I 133 

125 

1122 

87120 

484 

I 134 

9 

79 

5500 

484 

1-135 

64 

440 

27060 

462 

01-135 

64 

440 

27060 

462 

Xe  133 

0 

0 

0 

0 

Xe  135 

0 

0 

0 

0 

Cs  134 

3300 

11000 

0 

462- 

Cs  136 

1210 

3300 

0 

242 

Cs  137 

1804 

7700 

0 

242 

Ba  140 

La  140 

1584 

26840 

0 

12100 

184 

3740 

0 

12100 

Ce  141 

178 

5060 

0 

2420 

Ce  143 

74 

2200 

0 

6600 

Ce  144 

1452 

47520 

0 

24200 

Pr-143 

249 

6380 

0 

5280 

Hd-147 

255 

5280 

0 

4840 

Pm  147 

79 

2640 

0 

1210 

Pu  238 

66880 

2200000 

0 

9460 

Pu  239 

62260 

2112000 

0 • 

9460 
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Table  2 


WIND  FREQUENCIES  BY  STABILITY  CLASS  FOR 
WINDS  BLOWING  TO  1 SECTOR  (hourly  observations/%) 


A 

B 

C 

N 

89/4.93 

81/4.49 

73/4.04 

NNE 

31/3.73 

15/. 81 

36/4.33 

NE 

3/1.31 

4/. 69 

7/3.06 

ENE 

1/.81 

1/4.19 

1/.81 

E 

1/.69 

1/4.45 

3/2.07 

ESE 

4/1.29 

13/ 

14/4.52 

SE 

13/2.75 

21/4.42 

26/5.51 

SSE 

23/3.63 

28/ 

33/5.21 

S 

44/5.14 

23/2.69 

32/3.74 

ssw 

26/2.69 

23/2.38 

36/2.73 

sw 

50/11.71 

29/6.79 

29/6.79 

wsw 

10/6.33 

14/8.86 

16/10.13 

w 

5/5.0 

3/3.0 

8/8.0 

WNW 

6/7.14 

5/5.95 

9/10.71 

NW 

8/6.56 

5/4.1 

9/7.38 

NNW 

20/5.08 

18/4.57 

23/5.84 
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Table  2 continued 


D 

E 

F 

N 

613/33.96 

851/47.15 

88/4.88 

KNE 

240/28.88 

396/47.65 

100/12.03 

HE 

62/27.07 

111/48.47 

28/12.23 

ENE 

34/27.42 

55/44.35 

24/19.35 

E 

67/46.21 

48/33.1 

20/13.79 

ESE 

151/48.71 

102/32.9 

19/6.13 

SE 

188/39.83 

175/37.08 

47/9.96 

SSE 

231/37.44 

269/42.43 

43/6.78 

S 

336/39.25 

372/43.46 

- 46/5.37 

SSW 

448/46.42 

383/39.69 

45/4.66 

SH 

158/37.0 

146/34.19 

.14/3.28 

WSW 

56/35.44 

55/34.81 

7/4.43 

W 

39/39.0 

39/39.0 

6/6.0 

WNW 

28/33.33 

32/38.1 

4/4.76 

NW 

45/36.89 

40/32.79 

13/10.66 

NMW 

123/31.22 

182/46.19 

26/6.6 

Table  3 
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WIND  SPEEDS  BY  STABILITY  CLASS  (mph/  — ) 


A 

B 

C 

N 

12,71/5.68 

12.24/5.47 

12.61/5.64 

NNE 

11.37/5.08 

11.8/5.28 

10.81/4.83 

NE 

7.5/3.35 

12.88/5.76 

9.43/4.22 

ENE 

5.5/2.46 

15/6.71 

15.0/6.71 

E 

5.5/2.46 

15/6.71 

11.83/5.29 

ESE 

15/6.7 

12.58/5.62 

13.21/5.91 

SE 

16.42/6.71 

14.24/6.37 

12.56/5.61 

SSE 

18.09/8.09 

15.66/7.0 

14.17/6.33 

S 

17.97/8.03 

14.26/6.37 

14.44/6.46 

SSW 

12.46/5.57 

13.11/5.86 

11.25/5.03 

SW 

10.23/4.57 

10.4/4.65 

8.64/3.86 

wsw 

6.05/2.7 

5.07/2.27 

5.41/2.42 

n 

6.6/2.95 

2/.  89 

5.06/2.26 

WNW 

4.33/1.94 

2/.  89 

4.33/1.94 

NW 

4.19/1.87 

5/2.24 

2.78/1.24 

NNW 

5.75/2.57 

6.42/2.87 

4.59/2.05 
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Table  3 continued 


D 

E 

F 

N 

12.87/5.75 

11.26/5.03 

10.95/4.9 

NNE 

10.22/4.57 

9.53/4.26 

9.14/4.09 

NE 

9. 4/4. 2 

6.77/3.03 

6.41/2.87 

ENE 

7.82/3.5 

6.02/2.69 

5.67/6.12 

E 

9. 4/4. 2 

6.01/2.69 

4.78/5.23 

ESE 

12.03/5.38 

8.8/3.93 

6.21/2.78 

SE 

11.3/5.05 

9.81/4.39 

6.71/3.0 

SSE 

14.37/6.42 

11.95/5.34 

6.84/3.06 

S 

12.74/5.7 

9.76/4.36 

5.89/2.63 

SSH 

11.5/5.14 

8.0/3.58 

5.3/2.37 

SW 

7.53/3.37 

5.55/2.48 

4.64/2.07 

WSW 

5.39/2.41 

4.13/1.85 

3.0/1.34 

H 

4.21/1.88 

3.87/1.73 

3.17/1.42 

WNW 

4.41/1.97 

3.86/1.73 

3.75/1.68 

m 

5.08/2.27 

3.8/1. 7 

6.27/2.8 

NNW 

7.29/3.26 

6.69/2.99 

7,42/3.32 

1 
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Table 

4 

1 

YEAR 

1990  EXPECTED  POPULATION 

13 

Radius  (mi) 

N/S 

NNE/SSN 

NE/SW 

1 

(0-1) 

3 

0 

0 

1 

0 

3 

100 

1 

I (1-2) 

3 

0 

17 

13 

175 

272 

(2-3) 

23 

' 44 

'55 

i 

32 

215 

122 

i 

(3-4) 

48 

52 

139 

73 

543 

- 348 

(4-5) 

180 

1039 

107 

104 

505 

204 

(5-10) 

1337 

5109 

10,915 

872 

2580 

4446 

1 

(10-20) 

3877 

7869 

10,639 

1 

19,868 

58,206 

8790 

1 

(20-25) 

6007 

11,400 

6006 

i 

14,000 

14,951 

3765 

I 

(25-30) 

6007 

11,400 

6006 

I 

14,000 

14,951 

3765 

(30-40) 

114,379 

133,932 

11,436 

182,797 

18,037  ■ 

17,213 

(40-45) 

103,014 

86,034 

5719 

' 

86,843 

38,066 

9389 

: 

(45-50) 

103,014 

86,034 

. 5719 

! 

86,843 

38,066 

9389 

Radius  (mi) 

ENE/WSW 

E/W 

ESE/WNW 

SE/NW 

WWE/NNW 

(0-1) 

0 

3 

3 

0 

0 

41 

11 

• 3 

0 

3 

(1-2) 

193 

127 

213 

308 

245 

0 

11 

129 

0 

8 

(2-3) 

143 

73 

155 

143 

167 

237 

111 

15 

92 

82 

(3-4) 

130 

110 

123 

92 

178 

153 

66 

90 

59 

82 

(4-5) 

130 

167 

76 

110 

216 

579 

343 

211 

262 

85 

(5-10) 

1320 

923 

765 

1168 

595 

15,094 

2758 

2213 

1257 

1266 

(10-20) 

6395 

2769 

1936 

3720 

6499 

6998 

989 

1794 

2550 

3796 

(10-15) 

1791 

4900 

2758 

8067 

6500 

1565 

1117 

• 920 

1167 

1715 

(25-30) 

1791 

4900 

2758 

8067 

6500 

1565 

1117 

920 

1167 

1715 

(30-40) 

47,218 

5085 

9232 

10,032 

28,398 

1737 

4008 

3371 

2628 

6491 

(40-45) 

17,973 

1957 

28,816 

12,517 

24,012 

5759 

1271 

2897 

6298 

6178 

(45-50) 

17,973 

1957 

28,816 

12,517 

24,012 

5759 

1271 

2897 

6298 

6178 
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Table  5 


MESH  POINT  DISTANCES  IN  METERS^^ 


Radl\is 

Midpoint 

Lower  Radius 

Upper  Radius 

1 

805 

0 

1609 

2 

2414 

1609 

3219 

3 

4023 

3219 

4828 

4 

‘ 5633 

4828 

6437 

5 

7242 

6437 

8047 

6 

12,070 

8047 

16,093 

7 

24,140 

16,093 

32,187 

8 

36,210 

32,187 

40,234 

9 

44,257 

40,234 

48,280 

10 

56,327 

48,280 

64,374 

11 

68,397 

64,374 

72,420 

12 

76,444 

72,420 

80,467 
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Tetble  6 

NUCLIDE  SOURCE  DATA^ 
For  1200  MW(e)  After  550  Days  of  Full 


Isotope  Half-Life(d) 


KR85 

3.9E03 

KR85* 

1.8E-01 

KR87 

5.3E-02 

KR88 

1.16E-01 

Sr89 

5.06E01 

Sr90 

1.05E04 

Sr91 

4.0E-01 

Y90 

2.7E00 

Y91 

5.9E01 

Zr95 

6.55E01 

Zr97 

7.0E-01 

Nb95 

3.5E01 

Mo99 

2.8E00 

Tc99* 

2.5E-01 

Rul03 

4.0E01 

Rul05 

1.8E-01 

RulOe 

3.68E02 

Rhl05 

1.5E00 

Tel29 

4.8E-02 

Tel29* 

3.41E01 

Tel31* 

1.25E00 

Tel32 

3.25 

1131 

8.05C00 

MI  131 

8.05E00 

1132 

l.OE-01 

1133 

8.75-01 

1134 

3.6E-02 

1135 

2.8E-01 

MI  135 

2.8E-01 

XE133 

5.3E00 

XE135 

3.8E-01 

CS134 

7.52E02 

CS136 

1.29E01 

CS137 

1.1E04 

BA140 

1.28E01 

LA140 

1.66E00 

CE141 

3.28E01 

CE143 

1.37E00 

CE144 

2.85E02 

PR143 

1.36E01 

ND147 

l.lEOl 

PM147 

9.6E02 

PM149 

2.2E00 

PU238 

3.2E04 

PU239 

8.7E06 

Power  Operation 


Constant(sec~^) 


3.9E-03 

4.46E-05 

1.51E-04 

6.91E-05 

1.59E-07 

7.64E-10 

2.01E-05 

2.97E-06 

1.36E-07 

1.22E-07 

1.15E-05 

2.29E-07 

2.86E-06 

3.21E-05 

2.01E-07 

4.46E-05 

2.18E-08 

5.35E-06 

1.67E-04 

2.35E-07 

6.42E-06 

2.47E-06 

9.96E-07 

9.96E-07 

8.02E-05 

9.17E-06 

2.23E-04 

2.87E-05 

2.86E-05 

1.51E-06 

2.11E-05 

1.08E-08 

6.22E07 

7.29E-10 

6.27E-07 

4.83E-06  • 

2.45E-07 

4.86E-06 

2.81E-08 

5.89E-07 

7.29E-07 

8.36E-09 

3.65E-06 

2.51E-10 

9.22E-13 


78 
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Table  6 continued 


Isotope 

Source (Ci) 

Fraction  Released 

Total  Release (Ci) 

KR85 

7.2E05 

.8 

5.76E05 

KR85* 

3.12E07 

.8 

2.50E07 

KR87 

6.24E07 

.8 

4.99E07 

KR88 

9.12E07 

.8 

7.30E07 

Sr89 

1 . 32E08 

.05 

6.60E06 

Sr90 

6.24E06 

.05 

3.12E05 

Sr91 

1.S6E08 

.05 

7.80E06 

Y90 

6.24E06 

.003 

5.04E05 

Y91 

1.68E08 

.003 

5.04E05 

Zr95 

1.92E0S 

.003 

5.76E05 

Zr97 

1.92E08 

.003 

5.76E05 

Nb95 

1.92E08 

.003 

5.76E05 

Mo99 

1.92E08 

.4 

7.68E07 

Tc99* 

1.68E08 

.4 

6.72E07 

RulOa 

1.2E08 

.4 

4.80E07 

Rules 

6.96E07 

.4 

2.78E07 

Ruloe 

2.28E07 

.4 

9.12E06 

RhlOS 

6.96E07 

.4 

2.78E07 

Tel29 

3.36E07 

.4 

1.34E07 

Tel29* 

1.2E07 

.4 

4.80E06 

Teiai* 

1.8E07 

.4 

7.20E06 

Tel32 

1.44E08 

.4 

5.76E07 

1131 

1.02E08 

.6 

6.12E07 

MZ131 

1.02E08 

.006 

6.12E05 

1132 

1.44E08 

.6 

8.64E07 

1133 

2.0HE08 

.6 

1.22E08 

1134 

2.4E08 

.6 

1.44E08 

“1135 

1.8E08 

.6 

1.08E08 

MI135 

1.8E08 

.006 

1.08E08 

XE133 

2.04E08 

.8 

1.63E08 

XE135 

3.12E06 

.8 

2.50E07 

CS134 

a.04E06 

.4 

8.16E05 

CS136 

7.2E06 

.4 

2.d8 "06 

CS137 

6.96E06 

.4 

2.78E06 

BA140 

1.92E08 

.05 

9.60E06 

U140 

1.92E08 

.003 

S.76E05 

CE141 

1.92E08 

.003 

5.76E05 

CE143 

1.8E08 

.003 

5.40E05 

CE144 

1.32E08 

.003 

3.96E05 

PR143 

1.8E08 

.003 

5.40E05 

MD147 

7.2E07 

.003 

2.16E05 

PM147 

2.04E07 

.003 

6.12E04 

PH149 

4.8E07 

.003 

1.44E05 

PU238 

1.2E05 

.003 

3.60E02 

PU239 

1.2E04 

.003 

3.60E01 
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Ttble  7 


ACTDEN 

AP 


DOSFAC 

DR 

N 

NEVAC 


NORGAN 

TAVL 

TD 

nME 

TUG 

TOTPDP 

RUN 

XUMBD 

XMANRll 

POPUL 

FX 

FX2 

POP2 


NEW  DEFINITIONS  INTRODUCED  IN  NAIN  PROGRAM 

2 

Ground  deposition  in  pCl/m 

In  the  evacuation  mode,  the  fraction  of  the  population  not 
evacuated. 

2 

Dose  factors  (Rem/hr/Ci/ra  ). 

Dose  rate  rem/hr. 

Number  of  radll--must  be  an  even  nusd>er 

Evacuation  mode  0 > No  evacuation 

-1  « Evacuation  outside  farthest  sector 
1 • Evacuation  outside  radius  a/2 

(away  from  the  cloud's  direction) 

Organ  of  Interest  for  ground  deposition  dose 
0 > VB,  1 > Skin. 

Tima  between  awareness  of  Impending  core  melt  and  leakage 
in  days. 

Tbtal  Integrated  dose  due  to  ground  deposition  nren. 

Integrated  tlme(days)  spent  standing  on  contasdnated  ground. 

Tims  lag  associated  with  Interpretation  of  data  and  Issuance 
of  warning  to  evacuate. 

Total  population  being  considermd. 

Nsasure  of  evacuation  rate  (days). 

Decay  factor  per  nuclide  (days~^). 

Total  Sian  rem  due  to  1 nuclide. 

e* 

Population  In  each  mesh  point  modified  by  evacuation. 
Population  fraction  remaining  at  the  tlise  of  cloud  arrival. 
Population  fraction  remaining  after  TIME  days  of  evacuation. 
Population  remaining  after  TIME  days  of  evacuation. 
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Table  8 


DEFINITIONS  - SUBROUTINE  RIVCON 


AAVCRS  Avex^gc  cross  sectional  area  of  the  river  in  a reach  concen- 

tration (pCi/sec)  in  water  due  to  the  nuclide  in  the  previoxis 

reach. 

CCI  Concentration  due  to  the  contribution  of  direct  deposition 

and  the  plant  discharge  (pCi/1) 

CONCEN  Concentration  in  reach  before  water  treatment  (Efi). 

CROSS  River  cross  sectional  area 

PINAL  Concentration  (pCi/1)  of  isotope  LL  in  each  reach  for  treated 

water. 

FRACT  Fraction  of  dissolved  nuclide  passing  through  water  treatment 

facilities  class  3 or  4 (If  unknown,  set  * 1.). 

IIRCCH  Nuaber  of  reaches 


IKACN  Co«atar  for  IIRECH 

NINOUT  Tide  direction  ■ 0 if  tide  is  going  out  at  beginning  of 

accident.  Tide  direction  ■ 1 if  tide  is  going  in  at  beginning 
of  accident. 

HDCF  Counter  for  NDCFTH 

■DEPTH  lAaber  of  depths /profile  (3,  S,  7,  9 or  11). 

■NPROr  Counter  for  NPROF 

■PROF  Nuaber  of  profiles /reach  (20  aaximua). 


i 

1 

\ 

t 


■■1,11  Reach  nxsaber  and  aesh  point  location  of  the  reach.  Numb  should 

start  at  1 end  of  the  river  and  each  reach  met  numbered 
oooaecutively  whenever  a new  aesh  point  is  reached. 

OUTVEL  Velocity  (ai/hr)  of  water  with  the  tide  going  out. 

PLTDEN  pCi/sec  added  to  the  water  due  to  plant  discharge  for  each 

nuclide . 

QQAVE  Average  flow  rate  (l/sec)  in  each  reach. 

e 

SOLUB  Solubility  factor  for  each  nuclide. 

2 

SURFA  Surface  area  of  river  in  the  reach  (ft  ). 

\ 
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Table  8 continued 


TIMIN 


Time  (hr)  the  tide  flows  in. 


TIMOUT 


Time  (hr)  the  tide  flows  out. 


TOTFIN 


Total  concentration  of  nuclides  (pCi/1)  in  processed  water 
in  reach. 


TVLTIM 


Time  for  a packet  of  water  (avg)  to  traverse  the  length  of 
the  reach. 


TWIDTH 


WATFAC 


XLRCH 


XDEPTH 


XIMVEL 


XXTIN 


River  width . 

Water  treatment  facility  in  the  reach:  0 = no  plant, 

1 = class  1 plant,  2 = class  2 plant,  3 = class  3 plant, 

4 s class  4 plant. 

Length  of  reach  (miles). 

River  depth 

Velocity  (mi /hr)  of  water  with  the  tide  going  in 

Time  delay  desired  before  computing  water  concentrations  only 
1 XXTIM/run  allowed,  must  be  in  whole  days,  0 = RIVDOS 
calciilations  not  desired. 


RADIUS 


MM  location  of  reach. 


SECTOR 


II  location  of  reach. 


COOMC 


Total  concentration  (pCi/1)  in  a reach  for  nuclide  LL. 
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Table  9 


DEFINITIONS  - SUBROUTINE  FISHR 


CF(ll)  Concentration  factor  for  fresh  water  fish  (pCi/kg/pCi/1) . 

FISHCN 

(iREACHjLL)  Concentration  (pCi/kg)  of  nuclide  LL  in  reach  IREACH  in  fish. 
FISHT 

(IREACH)  Total  concentration  (pCi/kg)  in  reach  IREACH  in  fish. 


DATA  REQUIRED  - SUBROUTINE  FISHR 


L Cards 


Column 


Title  Format 

Concentration  Factor 
for  fresh  water  fish 


1-10 


FIO.O 


Table  10 


DEFINITIONS  - SUBROUTINE  LAKEDS 


BOATDF 


BOATSK 


DESDOS 


DELAY 


DOSBSK 


DOSBWB 


DOSSSK 


Time  spent  boating  on  the  contaminated  lake  (hrs). 

Whole  body  dose  factors  for  boating  on  the  contaminated 
lake  ( mrem/hr/pCi/1). 

Skin  dose  factor  for  boating  on  the  contaminated  lake 

Dose  desired:  0=  whole  body,  1 = skin,  2 = both  whole  body 
and  skin. 

Time  (days)  after  lake  contamination  the  dose  is  desired. 
Skin  dose  due  to  boating  on  the  contaminated  lake. 

Whole  body  dose  due  to  boating. 

Skin  dose  due  to  swimming  in  a contaminated  lake. 


DOSSWB 


FINAL 


SETUP 


SKINDF 


SWIMDF 


TOTDOS 


VOLUME 


Whole  body  dose  due  to  swimming. 

Concentration  (pCi/1)  of  treated  water. 

Indicator  - I =run  subprogram  I times  for  I lakes, 

0 or  blank  card  = subprogram  not  required. 

Skin  dose  factor  for  swimming  in  a contaminated  lake. 

Total  concentration  (pCi/1)  of  nuclides  for  treated  water. 

2 

Lake  surface  area  (ft  ). 

Time  spent  swimming  in  the  contaminated  lake  (hrs). 

Whole  bo(fy  dose  factors  for  swimming  in  the  contaminated  lake. 
Total  concentration  (pCi/1)  of  nuclides  for  untreated  water. 
Lake  volume  in  gallons. 

Lake  concentration  (pCi/1)  of  nuclide  LL. 


84 


I 

I 


I 

I 

1 

I 

I 

t 

I 

I 

4 

I 

i 

i 

I 


BIO 

BF 

CD 

CBIO 

DAY 

DRINK 

INC 

LAMEFF 

LEAVE 

XLD 

NUKBER 

P 

QF 

R 

SD 

TAREA 

TOTAL 

TOTDAY 

VOL 

H 

YF 


1 

i 

■'  I 

Table  11  j 

I J 

DEFINITIONS  - SUBROUTINE  MILK  j; 

! 1 
I ] 

t ] 

Biological  half  life  in  forage  crop  (days).  \\ 

! 

Forage  crop  concentration  factor  (pCi/kg  per  pCiAgsoil).  ■ 

I 

Milk  concentration  (pCi/1).  ;j 

Biological  half  life  in  the  cow. 

Time  (days)  after  accident  release. 

Drinking  water  type  - trough,  river,  lake 
0 = trough  (undergroimd  water),  1 = river,  2 = lake. 

Increments  (days)  desired  in  dose  calculations.  ; 

.693 

Effective  decay  constant  = X , . , . , + ■= — : — rr— ; • i ■ 

•'  radiological  Th  Biological  = 14  days 

j , 

Dair;'  cow  evacuation  mode  after  exposure. 

(0  = no  evacuation,  1 = evacuation  after  1 day). 

Ingestion  of  drinking  water  by  the  dairy  cow  (L/day). 

Number  of  mesh  points  not  serviced  by  undergroimd  water  supply.  I 

If  = -1,  go  to  return.  j 

2 

Constant:  224kg  soil/m  plowlayer. 

Fresh  forage  ingestion  by  dairy  cow  (kg/day). 

Deposition  retention  factor.  , 

Coefficient  of  transfer  from  diet  to  milk  (pCi/1  per  pCi/day).  | 

I 

2 

Typical  trough  area  (ft  ).  \ 

Total  milk  concentration  in  each  mesh  point  (pCi/1).  i 

Maximum  number  of  post  accident  days  to  calciilate  doses. 

Typical  trough  capacity  (gal). 

Concentration  (pCi/1)  in  drinking  water.  i 

2 I 

Forage  crop  yield  (kg/m  ). 


I 
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REM/PERSON  TO  THE 

Table  12 

THYROID  FOR  SELECTED  DIRECTIONS 

Average  Radixis 

(mi) 

SSW 

H 

SSE 

N 

.5 

3.0E05 

7.1E05 

— 

2.3E05 

1.5 

1.1E05 

2.5E05 

8.4E04 

8.7E04 

2.5 

5.7E04 

9.5 

4.3E04 

4.3E04 

3.5 

3.5E04 

7.2E04 

2.7E04 

2.7E04 

4.5 

2.4E04 

4.9E04 

1.9E04 

1.9E04 

7.5 

1.1E04 

2.3E04 

8800 

8900 

15 

4100 

7600 

3200 

3200 

22.5 

2100 

4200 

1700 

1700 

27.5 

1600 

2600 

1300 

1300 

35 

1100 

1800 

900 

900 

42.5 

790 

1200 

660 

670 

47.5 

660 

1000 

560 

570 

Table  13 


I 


THYROID  MAN-REM 


i 

I 

i 


I 

i 

i 


Major  Contributing  Nuclides 


M 

2.71E08 

1131  = 60% 

NNE 

4.27E08 

1133  = 30% 

NE 

3.01E08 

1135  = 7% 

ENE 

2.02E08 

E 

8.72E07 

ESE. 

9.10E07 

SE 

1.26E08 

SSE 

1.28E08 

S 

4.50E08 

SSW 

4.98E08 

SH 

2.59E08 

WSW 

4.36E08 

H 

1.33E08 

HNW 

1.33E08 

NW 

9.86E07 

NMH 

9.80E07 

I 

i 
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Table  14 


REM/PERSON  TO  THE  GI  TRACT  FOR  SELECTED  DIRECTIONS 


Average  Radius 


(mi) 

ssw 

W 

SSE 

N 

.5 

2.06E04 

4.86E04 

1.46E04 

1.59E04 

in 

• 

7700 

1.64E04 

5700 

5880 

2.5 

3810 

7700 

2880 

2910 

3.5 

2310 

4500 

1770 

1780 

4.5 

1470 

2980 

1220 

1220 

7.5 

730 

1340 

579 

576 

15 

245 

400 

201 

201 

22.5 

120 

160 

101 

103 

27.5 

84 

110 

72 

75 

35 

56 

71 

50 

52 

42.5 

39 

45 

35 

37 

47.5 

32 

36 

29 

31 

Major  Contributing  Nuclides 
TE132  30% 

MO  99  22% 

RU103  10% 

RU106  9% 

BA140  5% 
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Table  15 

GI  TRACT  MAN- REM 


N 

1.55E07 

NNE 

2.40E07 

NE 

1.79E07 

ENE 

1.15E07 

E 

5.27E06 

ESE 

5.44E06 

SE 

7.58E06 

SSE 

7.66E06 

S 

2.46E07 

SSW 

2.90E07 

SW 

1.55E07 

WSW 

2.56E07 

W 

7.72E06 

WNW 

7.72E06 

NW 

5.32E06 

NNW 

5 . 35E06 

r 


Table  16 


REM/PERSON  TO  THE  LUNGS  FOR  SELECTED  DIRECTIONS 


Average  Radius 


(mi) 

_ ssw 

W 

SSE 

N 

.5 

2.22E04 

2.52E04 

1.57E04 

1.71E04 

1.5 

8360 

1.78E04 

6150 

6340 

2.5 

4140 

8480 

3120 

3150 

3.5 

2520 

5000 

1920 

1930 

4.5 

1724 

3340 

1330 

1330 

7.5 

812 

1530 

637 

633 

15 

278 

474 

225 

225 

22.5 

137 

203 

115 

118 

27.5 

98 

139 

84 

86 

35 

67 

91 

58 

60 

42.5 

47 

59 

42 

43 

47.5 

39 

47 

34 

36 

Table  17 


LUNG  INHALATION  MAN- REM 


N 

1.82E07 

NNE 

2.87E07 

NE 

2.02E07 

ENE 

1.35E07 

E 

5.92E06 

ESE 

6.56E06 

SE 

8.63E06 

SSE 

8.92E06 

S 

2.84E07 

SSW 

3.21E07 

SW 

1.75E07 

WSW 

2.88E07 

W 

8.78E06 

WNW 

7.54E06 

NW 

1.09E07 

NNW 

5.28E06 

! 

I 

i 

i 


I 
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Table  18 


WHOLE  BODY  REM/ PERSON  DUE  TO  CLOUD  SUBMERSION  FOR  SELECTED  DIRECTIONS 


Average  Radiiis(ini) 

ssw 

W 

SSE 

N 

.5 

760 

1730 

530 

580 

1.5 

279 

550 

201 

213 

2.5 

133 

256 

102 

103 

3.5 

79 

13.8 

62 

62 

4.5 

53 

90 

42 

42 

7.5 

23 

35.6 

18.8 

19 

15 

6.4 

8.5 

5.58 

5.6 

, 22.5 

2.87 

3.39 

2.60 

2.73 

27.5 

1.91 

2.07 

1.78 

1.80 

35 

1.17 

1.22 

1.1 

1.15 

42.5 

.77 

.75 

.76 

.78 

47.5 

.61 

.57 

.60 

• 62 
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Table  19 

WHOLE  BODY  CLOUD  SUMBERSION  MAN-REM 


N 

3.08E05 

NNE 

5.93E05 

ME 

4.99E05 

ENE 

E.08E05 

E 

1.50E05 

ESE 

1.60E05 

SE 

2.18E05 

SSE 

2.17E05 

S 

5.70E05 

ssw 

7.60E05 

sw 

4.60E05 

wsw 

6.90E05 

w 

2.16E05 

WNW 

1.73E05 

NW 

1.35E05 

NNW 

1.36E05 

Major  Contributing  Nuclides 
1135  30% 

1132  19% 

1133  12% 

1134  9% 

1131  5% 


Table  20 

WHOLE  BODY  REM/PERSON  DUE  TO  GROUND  DEPOSITION  FOR  SELECTED  DIRECTIONS 


Average  Radius (mi) 

ssw 

W 

N 

.5 

46.4 

93.2 

80.3 

35.6 

1.5 

17.4 

37.2 

31.9 

13.2 

2.5 

8.61 

17.6 

15.2 

6.60 

3.5 

5.24 

10.4 

9.03 

4.02 

4.5 

3.58 

6.89 

6.05 

2.77 

7.5 

1.69 

3.17 

2.81 

1.32 

15 

.573 

1.44 

.901 

.467 

22.5 

00 

CM 

.647 

.395 

.244 

27.5 

.199 

.453 

.275 

.179 

35 

.136 

.179 

.178 

.123 

42.5 

.095 

.115 

.114 

.090 

47.5 

.078 

.091 

.090 

.074 

r 


W ' 


k 


1 

Table  20A 

MAN- REM  DUE  TO  GROUND  DEPOSITION  FOR  ONE  DAY  ’ 


N 

48,950 

NNE 

57,180 

NE 

41,850 

ENE 

27,140 

E 

12,270 

ESE 

12,680 

SE 

17,590 

SSE 

17,850 

j 

S 

58,870 

i 

SSW 

67,980 

j 

SW 

35,990 

WSW 

60,090 

1 

w 

18,100 

1 

WNW 

16,740 

NW 

12,860 

1 

i 

NNW 

13,160  • 

t 

1 

i 

95 
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Table  21 


DOSE  FACTORS  FOR  EXPOSURE  ON  CONTAMINATED  GROUND 


Isotope 

Average 

WB  Ground 

Dose  Factors 
(rem/hr/Ci/m^)^ 

Skin  Dose 
Factors 

KR  85 

0 

0 

0 

KR85* 

.16 

3.0 

6.6* 

KR87 

.82 

12.0 

15* 

KR88 

2.21 

30.0 

39* 

SR89 

0 

0 

0 

SR90 

0 

0 

0 

SR91 

.75 

14.2 

21* 

Y90 

1.739^ 

.01 

.013* 

Y91 

1.21^ 

.01 

.014* 

ZR95 

.76 

10.0 

12.4 

ZR97 

.24 

11.9 

23* 

Nb95 

.77 

10.2 

13.16 

Mo99 

.19 

3.6 

8* 

Tc99* 

.14 

2,0 

4.5* 

Rul03 

.49 

7.2 

8.4 

Rul05 

.79 

9.0 

13.4* 

Rul06 

.20 

3.0 

3.6 

Rhl05 

.02 

1.4 

7 

Tel29 

.07 

1.5 

2.3* 

Tel29* 

.10 

1.5 

3.3* 

Tel31* 

1.49 

16.8 

22.7* 

Tel32 

.28 

3.4 

4.0 

1131 

.39 

5.6 

6.8 

* Calculated  from  average  energy  at  skin  depth  of  7x10 


cm. 
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Table  21  continued 


[sotope 

Modified 

WB  Dos  e 
Factors 

Modified 

Skin 

Dose  Factors 

Soliability  in 
Cold  Water 

KR85 

0 

0 

S 

KR85* 

1.0 

2.2 

S 

KR87 

3.3 

5.0 

S 

KR88 

10 

13 

S 

SR8. 

0 

0 

D 

SR90 

0 

0 

D 

SR91 

4.7 

7.0 

D 

Y90 

.003 

.004 

SS 

D 

Y91 

.003 

.005 

SS 

D 

Zr95 

3.3 

4.1 

I 

Zr97 

3.6 

7.6 

I 

Nb95 

3.4 

4.5 

I 

Mo99 

1.2 

2.6 

I 

Tc99* 

.66 

1.5 

I 

Rul03 

2.4 

GO 

CM 

D 

RulOS 

F.O 

4.5 

D 

RulOe 

1.0 

1.3 

D 

RhlOS 

.46 

2.3 

r 

Tel29 

.5 

.76 

I 

Tel29* 

.5 

1.1 

I 

Tel31* 

5.6 

7.6 

I 

Tel32 

1.1 

1.3 

r 

1131 

1.9 

2.3 

s 

S = Soluble 

I = Insoluble 

SS  = Slightly  soluble 

D = Decomposes 
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Table  22 


Table  22 


Isotope  Fraction 


I 132 


Mil35 

Xel33 

Xel35 

Csl34 

Csl36 

Csl37 

Bal40 

Lal40 

Cel41 

Cel43 

Cel44 

Prl43 

Ndl47 

Pml47 

Pral49 

Pu238 

Pu239 


Blank  s imknown 
and  assumed  = 1 


p ’ ■ 1 

Table  23 

HUDSON  RIVER  DATA 

f * - Water  Treatr.ent  i 


imber 

Radius 

Sector 

Length (mi) 

2 

Surface  Area (ft  ) 

Facility  Type 

1 

7 

2 

11 

1.63E08 

o 

] 

; 

i 

2 

6 

2 

5.3 

8.46E07 

4^^ 

1 

3 

5 

3 

1.1 

1.17E07 

420 

J 

4 

4 

3 

1.0 

1.62E07 

0 

i 

5 

3 

3 

1.0 

2.27E07 

0 

1 

6 

2 

3 

.7 

1.87E07 

0 

1 

7 

2 

2 

1.0 

1.13E07 

0 

j 

8 

1 

3 

.3 

4.32E06 

0 

1 

9 

1 

4 

.3 

6.48E06 

0 

j 

10 

1 

5 

.3 

4.55E06 

0 

11 

1 

6 

.3 

5.40E06 

0 

12 

1 

7 

.3 

5.85E06 

0 

13 

1 

8 

.4 

3.60E06 

0 

14 

2 

8 

.7 

9.72E06 

0 

15 

2 

8 

1.0 

2.11E07 

0 

16 

3 

9 

1.1 

1.84E07 

0 

17 

4 

9 

1.0 

2.17E07 

0 

18 

5 

10 

1.0 

2.27E07 

0 

19 

6 

9 

5.1 

1.36E08 

0 

100 
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Table  24  continued 


Table  25 


24 

WATER  SUBMERSION  DOSE  FACTORS 


Isotope 

Skin 

mrem/hr 

pCi/1 

Whole  Body 
mrem/hr 
pCi/1 

KR  85 

1.8E-07 

4.7E-09 

KR  85* 

5.1E-07 

2.8E-07 

KR  87 

,4.6E-06 

2.7E-06 

KR  88 

4.1E-06 

3.3E-06 

Sr  89 

5.4E-07 

4.6E-09 

Sr  90 

1.5E-07 

5.4E-10 

Sr  91 

2.9E-06 

1.9E-06 

Y 90 

9.6E-07 

1.3E-08 

Y 91 

5.7E-07 

6.7E-09 

Zr  95 

1.8E-06 

1.5E-06 

Zr  97 

2.4E-06 

1.5E-06 

Nb  95 

1.6E-06 

1.4E-06 

Mo  99 

2.4E-07 

2.1E-07 

To  99* 

2.7E-07 

2.4E-07 

Ru  103 

l.lE-06 

8.9E-07 

Ru  105 

1.8E-06 

1.2E-06 

Ru  106. 

1.9E-06 

3.8E-07 

Rh  105 

3.0E-07 

1.7E-07 

Te  129 

7.0E-07 

1.9E-07 

Te  129* 

7.4E-07 

2.1E-07 

Te  131* 

2.7E-06 

2.2E-06 

Te  132 

4.8E-07 

4.0E-07 

I 131 

9.3E-07 

5.8E-07 

MI  131 

9.3E-07 

6.8E-07 

I 132 

5.5E-06 

4.4E-06 

I 133 

1.5E-06 

9.6E-07 

I 134 

5.5E-06 

4.2E-06 

I 135 

4.0E-06 

3.3E-06 

MI  135* 

4.0E-06 

3.3E-06 

Xe  133 

l.lE-07 

5.7E-08 

Xe  135 

7.9E-07 

4.5E-07 

Cs  134 

3.5E-06 

2.9E-06 

Cs  136 

4.8E-06 

4.1E-06 

Cs  137 

1.4E-06 

l.OE-06 

Ba  140 

7.6E-07 

4.9E-07 

La  140 

5.3E-06 

4.1E-06 

Ce  141 

2.4E-07 

1.3E-07 

Ce  143 

l.OE-06 

5.7E-07 

Ce  144 

1.4E-06 

8.6E-08 

Pr  143 

2.8E-07 

1.6E-09 

Nd  147 

5.0E-07 

2.8E-07 

Pm  147 

1.3E-08 

7.5E-11 

Pm  149 

3.5E-07 

1.5E-08 

Pu  238 

4.0E-09 

1.5E-10 

Pu  239 

1.7E-09 

1.2E-10 

103 


Table  26 


r 


WHOLE  BODY  CLOUD  SUBMERSION  DOSE  FACTORS 

WASH  1400  WASH  1400  WASH  1258  WASH  1258 

converted  to  Modified  by  Converted  to  Modified  by 
mrem/sec  1/3  Shielding  mrem/sec  1/3  Shielding 

Isotope  Ci/m^  Factor  Ci/m^  Factor 

KR 

85 

0 

0 

.6 

.2 

KR 

85* 

36 

13 

36 

12 

KR 

87 

360 

130 

361 

120 

KR 

88 

420 

140 

417 

139  ■ 

SR 

89 

0 

0 

.6 

.2 

SR 

90 

0 

0 

.1 

.03 

SR 

91 

160 

53 

247 

82 

Y 

90 

2 

.6 

1.7 

.6 

Y 

91 

2 

.6 

.9 

.3 

Zr 

95 

190 

63 

190 

63 

Zr 

97 

60 

20 

192 

64 

Nb 

95 

180 

60 

178 

59 

Mo 

99 

60 

20 

31 

10 

Tc 

99 

35 

12 

31 

10 

Ru 

103 

110 

37 

114 

38 

Ru 

105 

200 

67 

150 

50 

Ru 

106 

50 

16 

47 

16 

Rh 

105 

5 

1.6 

22 

7 

Te 

129 

18 

6 

27 

9 

Te 

129* 

25 

12.5 

24 

8 

Te 

131* 

375 

125 

278 

93 

Te 

132 

50 

17 

50 

17 

I 

131 

90 

30 

86 

29 

I 

132 

550 

183 

556 

185 

I 

133 

120 

■ 40 

122 

41 

I 

134 

600 

200 

556 

185 

I 

135 

420 

140 

417 

139 

MI 

135 

420 

140 

417 

139 

Xe 

133 

7 

2.3 

7 

2 

Xe 

135 

60 

20 

58 

19 

Cs 

134 

360 

120 

361 

120 

Cs 

136 

460 

153 

528 

176 

Cs 

137 

130 

43 

131 

44 

Ba 

140 

60 

20 

61 

20 

La 

140 

520 

173 

528 

176 

Ce 

141 

16 

5.3 

1-6 

5 

Ce 

143 

85 

28 

72 

24 

Ce 

144 

4 

1.3 

11 

4 

Pr 

143 

0 

0 

.2 

.1 

Nd 

147 

45 

15 

36 

13 

Pm 

147 

0 

0 

0 

Pm 

149 

12 

4 

2 

.6 

Pu 

238 

0 

0 

.02 

.007 

Pu 

239 

0 

0 

.02 

.007 

104 


Table  27 


CRITICAL  ORGAN  DOSES 
DUE  TO  CONTAMINATED  WATER  INGESTION 
IN  ALCOVE  RESERVOIR 


Isotope 

Concentration 

(pCi/1) 

Critical  Organ 

Dose  Factor  to 
Critical  Organ 
(Rem/Ci) 

I 131 

2.5E07 

Thyroid 

1.94E06 

Cs  134 

9.6E06 

Whole  Body 

5.83E03 

Cs  137 

3.8E07 

Whole  Body 

8 . 6E04 

BA  140 

6.8E06 

GI 

4.43E05 

SR  89 

1.8E07 

Bone 

2.19E05 

SR  90 

4.3E06 

Bone 

4.03E06 

Ru  103 

1.1E08 

GI 

3.21E04 

Ru  106 

9.1E07 

GI 

3.47E05 

Isotope 

Dose  Factor  to 
Whole  Body 
( Rem/Ci ) 

Dose  to  Critical 
Organ  ( rem ) 

1 Liter  Ingested 

Dose  to  Whole 
Body  (rem) 

1 Liter  Ingested 

I 131 

3.53E03 

49 

.1 

CS  134 

5.83E03 

.1 

.1 

CS  137 

8.6E04 

3.3 

3.3 

BA  140 

1 . 3E03 

3.0 

.01 

SR  89 

8.78E03 

3.94 

.16 

SR  90 

1.77E06 

14.3 

7.6 

Ru  103 

8.64E01 

3.5 

.01 

Ru  106 

4.43E01 

31.6 

.004 

105 


Table  28 


CRITICAL  ORGAN  DOSES  DUE  TO  INGESTION  OF  1 LITER  OF  MILK  AT  VARIOUS  TIMES 

AFTER  THE  ACCIDENT 

Cows  are  assumed  to  be  fed  contaminated  feed  for  the  duration  of  time  of 
interest.  Distance  from  the  reactor  is  2-3  miles  north. 

Dose  Factor  Critical 


Isotope  (Rem/Ci)  Organ  Day  1 Day  11 


®^Sr 

2.19x10^ 

Bone 

19.5 

10.1 

^°Sr 

4.03x10® 

Bone 

84.6 

48.4 

99 

Mo 

1.93x10*^ 

GI 

6.6 

.33 

129^^ 

2.44x10** 

GI 

0 

0 

131^ 

1.94x10® 

Thyroid 

980 

254 

^^Vs 

5.83x10® 

WB 

1.34 

.82 

136cs 

7.55x10® 

WB 

.36 

.13 

^3^Cs 

8.6x10** 

WB 

80 

47 

^■^Oea 

4.43x10® 

GI 

8.4 

2.9 

Isotope 

Dose 
Day  21 

(rem) 

Day  31 

Day  41 

Day  51 

89c 

Sr 

5.3 

2.6 

1.4 

.72 

90- 

Sr 

29.4 

17.3 

10.1 

6.05 

^^Mo 

.017 

.0009 

. 04rarem 

0 

129Te 

0 

0 

0 

0 

131l 

63 

16 

4.12 

1.03 

'3"CS 

.47 

.27 

.16 

.09 

.044 

.015 

.0054 

.0019 

27 

16 

9.5 

5.76 

1.02 

.354 

.12 

.043 

106 


Table  29 


CRITICAL  ORGAN  DOSES  DUE  TO  INGESTION 
OF  1 LITER  OF  MILK  AT  VARIOUS  TIMES  AFTER  THE  ACCIDENT 

Cows  are  assumed  to  be  fed  uncontaminated  feed  after  1 
day's  ingestion  of  contaminated  forage.  Distance  from 
the  reactor  is  2-3  miles  north. 


Isotope 

Dose  Factor 
(Rem/Ci) 

Critical 

Organ 

Day  1 Day  11 

89_ 

Sr 

2.19x10^ 

Bone 

19.5 

.28 

Sr 

4.03x10® 

Bene 

84.6 

1.61 

^^Mo 

1. 93x10"* 

GI 

6.6 

.001 

129^e 

2.44xlo"* 

GI 

0 

0 

131j. 

1.94x10® 

Thyroid 

980 

3.3 

^^‘^Cs 

5.83x10® 

WB 

1.34 

.026 

7.55x10® 

WB 

.36 

.002 

^3"cs 

8.6xlo"* 

WB 

80 

1.55 

140„ 

6a 

4.43x10® 

GI 

8.42 

.053 

Isotope 

Day  21 

Ct‘ 

Day  31 

• •in') 

Day  41 

Day  51 

«^Sr 

.15 

.079 

.042 

.022 

^°Sr 

.967 

.605 

.359 

.218 

^^Mo 

0 

0 

0 

0 

129Te 

0 

0 

0 

0 

131l 

.83 

.21 

.05 

.014 

^3"CS 

.015 

.0093 

.0055 

.0034 

136c3 

.8mrem 

. 3mrem 

.Imrem 

0 

.95 

.56 

.34 

.21 

.019 

.0066 

.0023 

. 8mrem 

f 


Table  30 


i ■■'■J  IL"  ” 


PROBABILITY  OF  WIND  BLOWING  IN  DIRECTION 
LISTED  (%)  FOR  ALL  STABILITY  CLASSES^^ 


N 

23.58 

NNE 

10.85 

NE 

2.99 

ENE 

1.62 

E 

1.89 

ESE 

4.05 

SE 

6.17 

SSE 

8.28 

S 

11.18 

SSW 

12.6 

SW 

5.58 

wsw 

2.06 

w 

1.31 

WNW 

1.1 

NW 

1.59 

NNW 

5.15 

I 


Table  31 


FRESH  WATER  FISH  CONCENTRATION  FACTORS 

24 

Isotope  pCi/Kg/pCi/1 


KR 

85 

1 

KR 

85* 

1 

KR 

87 

1 

KR 

88 

1 

SR 

89 

30 

SR 

90 

30 

SR 

91 

30 

Y 

90 

25 

Y 

91 

25 

Zr 

95 

330 

Zr 

97 

330 

Nb 

95 

30000 

Mo 

99 

10 

Tc 

99* 

15 

Ru 

103 

10 

Ru 

105 

10 

Ru 

106 

10 

Rh 

105 

10 

Te 

129 

400 

Te 

129* 

400 

Te 

131* 

400 

Te 

132 

100031 

I 

131 

15 

MI 

131 

1524 

Values  vary  considerably  in  the  litera- 
ture. Most  conservative  values  used. 

24 

1 pp.  50-52  unless  otherwise  noted. 

31 

2 Value  from  p 114. 

3 Assumed  from 

24 

4 Assumed  from 


Table  32 


Isotope 

COEFFICIENTS  FOR  MILK  CONCENTRATION^ 

Coefficient  of  Transfer 

Concentration  Factor  (B)  From  Diet  to  Milk 

pCi/Kg  per  pCi/Kg  Dry  Soil  (pCi/1  per  pCi/day) 

KR 

85 

l.OE-00 

2E-2 

KR 

85* 

l.OE-00 

2E-2 

KR 

87 

l.OE-00 

2E-2 

KR 

88 

l.OE-00 

2E-2 

SR 

89 

2.0E-1 

IE- 3 

SR 

90 

2.0E-1 

IE- 3 

SR 

91 

2.0E-1 

lE-3 

Y 

90 

l.OE-1 

IE- 5 

Y 

91 

l.OE-1 

IE- 5 

Zr 

95 

1.7E-4 

5E-6 

Zr 

97 

1.7E-4 

5E-6 

Nb 

95 

9.4E-3 

2.5E-3 

Mo 

99 

1.3E-1 

7.5E-3 

Te 

99* 

l.OE-2 

2.5E-2 

Ru 

103 

l.OE-2 

l.E-6 

Ru 

105 

l.OE-2 

l.E-6 

Ru 

106 

l.OE-2 

l.E-6 

Rh 

105 

2.0E-10 

l.E-2 

Te 

129 

1.3E-00 

l.E-3 

Te 

129* 

1.3E-00 

l.E-3 

Te 

131* 

1.3E-00 

l.E-3 

Te 

132 

1.3E+0 

l.E-3 

I 

131 

2E-2 

l.E-2 

MI 

131 

2.0E-2 

l.E-2 

I 

132 

2E-2 

l.E-2 

I 

133 

2E-2 

l.E-2 

I 

134 

2.0E-2 

l.E-2 

I 

135 

2E-2 

l.E-2 

MI 

135 

2.0E-2 

l.E-2 

Xe 

133 

l.OE-00 

2.E-2 

Xe 

135 

l.OE-00 

2.E-2 

Cs 

134 

2E-3 

5.E-3 

Cs 

136 

2.0E-3 

5.E-3 

Cs 

137 

2E-3 

5.E-3 

Ba 

140 

5E-3 

6.E-4 

La 

140 

2.5E-3 

5.E-6 

Ce 

141 

5E-4 

2.E-5 

Ce 

143 

5.0E-4 

2.E-5 

Ce 

144 

5E-4 

2.E-5 

Pr 

143 

l.OE-00 

5.E-6 

Nd 

147 

l.OE-00 

5.E-6 

Pm 

147 

l.OE-00 

5.E-6 

Pm 

149 

l.OE-00 

5.E-6 

Pu 

238 

l.OE-11 

1.5E-6 

Pu 

239 

l.OE-11 

1.5E-6 
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PART  VIII 
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% 

APPENDIX 

\ 

J 

Listing  of  the  MODAIREM  Computer  Code 

Approximately  half  of  this  code  is  taken  from  the  AIREM  Computer 
‘ Code  written  by  J.  A.  Martin,  C.  B.  Nelson,  and  P.  A.  Cuny.  There  are 

several  minor  additions  in  the  body  of  the  AIREM  Code  and  the  original 
code  essentially  ends  with  statement  Number  1550  in  the  main  program 
plus  function  SIGZ.  A listing  of  the  code  in  its  entirety  follows. 
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QL(20),IS0I20)  ,0CF(?0)  ,XL>'DA(20)  ,nMIL(12,16,20), 

THYDOSt 12f 16) ,80Nn0S( 12,16) , SKMHOS « 1 2 , 1 6 ) , OPWB ( 1 2 , 16 ) , 
DPTHY{ 12,16) ,DPBON( 12, 16) ,DPSKN( 12 , 16 ) ,OPWBM(12) , 
DPSKNM(  12) ,W8D0S( 1 2 , 1 6 ) , OECCN ( 20 ) , 

OPTHYM{ 12) , *LOW{ 12 ) ,XUP< 12 ) ,Y( 13) ,P0P(12,16) , 

OMUL  (12, 16, 6, 20),  CM  100(12,16, 6, 20)  ,OPBONP  ( 1 2 ) , PV  ( 6 , 4 ) , 
TOM(  13)  ,ARE!V(  12)  ,00”!  (12,16)  ,QGA^<I  (12,16),HV(6,4), 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


QG(12,6),P0EPV(16,6) ,H0EPV(16,6)  ,F{13), 

LATWT( 20) ,CR  TORG( 20 ) ,NDEP ( 20 ) 

N gamma  ( 20)  ,NDTR(  20)  ,AGA^(  2 0,10)  , A81JN ( 20 , 10  ) 

GAMeN( 32) ,SIGZM(32) ,DOSEI (32,32)  ,ETITLE(20) 

DIR( 16),CHI0QS( 12) 

WTQGM( 13), WOG( 12,6), WOGAMK 12,16 ) ,WCGMI (12,16) 
CL0TIM( 12, 16) ,FX( 12,16) ,ACTDEN( 12, 16,20 ), OR (12,16) 
T0{12,16),P0?UL(12,16) 

XLAMBD( 20) 


I 


I 

I 





c 


OIMeNSION  FX2( I2f 16 ) ,POP2{ 12 .16) 
REAL  MONTHS 


) 


>1 


J 


. 


~) 


o 


DATA  OIR/'N  '.'NNE  *,'NE  •.•ENE  '.'E  '.'ESE  ‘.'SE  ‘.'SSE  *, 

* 'S  ‘.'SSW  '.'SW  '.'HSW  '.'W  '.'MNW  * , ' NW  • , ' NNW  •/ 

c 

EQUIVALENCE(CHIG0( 1,1, 1,1),DMIjl(1,1 ,1,1) ) 

C EQUIVALENCECWBDOSd,!)  ,THYDOSC  1,1)  , SKNDnS(  1 ,1 ) ,BONnOS  (1 1 1 > ♦ 

C IQGMni.D) 

C EQUIVALENCE(OPWB( 1 , 1 ) , DPTHY ( 1 , 1 ) ,DPBON ( 1 , 1 ) , OPSKM ( 1 , 1 ) , OGAMI ( 1 , 1 ) ) 

C EQUlVALENCE(OPWBM(  1)  .OPTHYMt  1)  ,OPBONM(l)  ,OPSKN»*(i)  ) 

EQUIVALENCE(Y(1),F(1),TQM{1) ) 

C 

PI*3. 141593 
C 

C READ  VARIABLE  INPUT  DATA 

C 

c 

DO  9001  MMsl,12 
DO  9001  11*1,16 
DO  9001  LL=1,20 
POPUL(MM,  II )=0. 

9001  ACTDEN<mm,II,LL)=0. 

HRITE(3,100) 

100  FORMAT ( 1H1,T32, 'PaOGRAM  AIREM04'// 

* 1X,T12,*FQR  INFORMATION  REGARDING  THIS  PROGRAM,*, 

♦ • CONTACT:*// 

♦ 1X,T20, *ENVIRONMENTAL  PROT-CTICN  AGENCY*/ 

• 1X,T20, *0FFICE  OF  RADIATION  PROGRAMS'/ 

* 1X,T20, *FIELD  OPERATIONS  DIVISICN*/ 

♦ IX,T20, *401  M ST  SW*/ 

♦ 1X,T20, *WASHINGTON,  DC'/ 

* IX, T20, *20460* ) 

MRITEI 3, 380) 

READ! 1, 110) ( FAC  I L, MONTHS, MON TH 1, NYR 1 ,M0NTH2 ,NYR2,ETHERM) 
WRITE(3,390) 

WRITE! 3, 112) ( FAC  I L, MONTHS, MONTH  1 ,NYR1,M0NTH2,NYR2,ETHERM) 

110  FORMAT! 5A4, F5.0, 1X,A4, I 5, IX, A4, 1 5,F10.0) 

112  FORMAT! 1X,5A4,F5.2,IX,A4, I5,1X,A4,I5,F10.0) 

120  REA0!1,300)  I , J , M, L , H, S I GMA X ,HOLDUP  , RA I NF  , WA SHCO 
WRITE(3,400) 

WRITE! 3, 300)  I , J , M, L ,H , S IGMA X ,HOLDUP , R A I NF , W ASHCO 
C 

C REAO/WRITE  WIND  FREQUENCY  (7)  BY  STABILITY  CLASS.  FIRST  LINE 

C IS  NORTH,  SECOND  IS  NNE,  ETC. 

C 

READ! 1,310) ! !FREO! I I,JJ)  ,JJ  = 1,J)  ,11=1,1) 

WRITE(3,410) 

T0TFR0=0. 

DO  122  11=1,1 
Cl-0. 

DO  121  JJ=1,J 

121  C1«C1^FRE0I II, JJ) 

T0TF«Q=T0TFPQ+C1 

122  WRiTE(3,124)  ! 0 IR ! 1 1 ) , ! FREQ ! 1 1 , J J ) . J J*1 . J > tCl ) 

124  FORMAT! 1X,A4,7F10.2) 

C 

C READ/WRITE  WIND  SPEEDS  BY  STABILITY  CLASS.  FIRST  LINE  IS  N, 

C SECOND  IS  NNE,  ETC. 
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) 


1 

;) 

) 


REAO( 1,310) ( (UBARf II, JJ) ,JJ=1,J)  ,11=1,1) 
WRITe(3,420) 

HRITE(3,125) <DIR( II ),(U8AR( I I , J J ) , J J»1 , J ) ,II»1,I 
126  FORHATl 1X,A4,6F10.2) 

C 

C READ/WRITE  POPULATIONS  CLOCKWISE  BY  SECTOR 

C FIRST  LINE  IS  POPULATION  IN  FIRST  8 SECTORS  IN 

C INNER  ANNULUS. 

C 

DO  130  MM=1,M 

130  READ(1,320)  ( POP  I MM, 1 1 ) , 1 1 = 1 , 1 ) 

131  NC0UNT=0 
133  CONTINUE 

WRITE(3,430) 

IF(NCOUNT.GT.O)GO  TO  144 
DO  140  MM=1,M 

WRITE (3, 132)  MM, (POP (MM,  II ) , 11=1,1  ) 

132  FORMAT! IX, 12, ( T5,8F10.0)  ) 

C1=0. 

DO  139  11=1,1 

139  C1=C1  + P0P (MM, H ) 

140  WRITE(3,141)  Cl 

141  F0RMAT(1X,T5,F10.0) 

GO  TO  147 

144  CONTINUE 
IF(NEVAC.EQ.1)G0  TO  147 
DO  146  MM=1,M 

WRITE(3,132)MM, (POPUL(MM,II ) ,11=1,1 ) 

C1»0. 

DO  145  11=1,1 

145  C1=C1+P0PUL (MM,  II ) 

146  WRITE(3,141)C1 

147  CONTINUE 

IF(NEVAC.LT.O)GO  TO  6610 
IFINCOUNT.EQ.O)GO  TO  6610 
MMM=M/2 

DO  6602  MM=1,MMM 

WRITE(3,  132)MM,  (P0PUL(MM,II  ) ,H*1,I  I 
Cl=0. 

00  6603  11=1,1 

6603  C1=C1+P0PUL (MM, II ) 

6602  WRITE(3,141)C1 

MMMM=MMM4-1 
DO  6601  MM=mmmM,M 

WRITE (3, 132) MM, (POP (MM, 1 1), II =1,1) 

C1«0. 

00  6604  11=1,1 

6604  C1*C1*P0P (MM,  II  ) 

6601  WRITEO,  141  )C1 
6610  CONTINUE 

IF(NCOUNT.GT.O)GO  TO  561 

READ! 1,330)((X(MM), XLOW( MM) ,XUP( MM) ) ,MM=1,M) 
WRITE! 3,440) 

WRITEO,  142  )(  mm,  ( X(  MM)  ,XL0W(  MM)  ,XUP(  MM)  ) ,MM=1  ,M) 

142  FORMAT(1X,I2,1X,3F10.0) 

C 

C READ  EXPANSION  COEFFICIENTS  FCR  PARTICULATES  (P) 

C 

C THE  INPUT  DATA  FORMAT  CAN  BE  USED  FCR  A DATA  FIT 


AND  HALOGENS  (H). 
UP  TO  FOUO  TERMS, 
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I 

MANY  FSAR"S  ASSUME  THAT  PARTICULATE  AND  HALCGEN  DEPOSITIONS  ARE 
DIFFERENT.  DATA  IN  MET.  AND  AT,  ENERGY  SHPw  A COMPLETE  OVERLAP. 
WITH  A LARGE  SPREAD  IN  BOTH  SETS  OF  DATA. 

MANY  FSAR"S  ASSUME  DEPOSITION  VELOCITY  IS  PROPORTIONAL  TO  WIND 
SPEED,  DATA  IN  MET,  AND  AT,  E^'E^GY  INOICATES  THE  OPPOSITE  TREND. 
IN  THIS  LIGHT  AND  WHEN  ACTUAL  DATA  IS  LACKING,  A SUGGESTION  IS 
TO  USE  A CONSTANT  O.Ol  meTER/SEC  FOR  BOTH  HALCGFN  AND  PARTICU- 
LATES. (FOLLOWED  BY  THREE  BLANK  INPUT  CARDS  IN  EACH  CASE). 

FORMAT  IS  AS  FOLLOWS 


CLASS/ 

TERM 

' A 

B 

C 

0 

E 

F 

INDICES 

CONSTANT 

Hll 

H21 

H3l 

H41 

H51 

H61- 

(CLASS, KOF) 

LINEAR 

H12 

H22 

H32 

H42 

• 

• 

QUADRATIC 

H13 

H23 

H33 

• 

• 

• 

CUBIC 

H14 

H24 

• 

• 

• 

H64 

H6,KK0F 

DEPOSITION  VELOCITY  IN  SECTOR  1 FOR  STABILITY  CLASS  A IS 
OEPV(l,l)  » Hll  ♦ H12  ♦UBAR(l,l)  ♦ H13  * UBAR(l,l)**2 

♦ H14  ♦ UBARd,!)**  3 (M/SEC) 

KKOF  a A 

REAO( 1, 370) ( (PV( JJ,KQF  ) , JJ=1  ,J ) ,KCFal  ,KKOF ) 

RE A0( 1,370) ( IHV( JJ,KOF ) , J J= 1 , J ) , KDF  = l , KKOF ) 

WRITE(3,150) 

150  F0RMAT(63H0PARTICULATE  AND  HALOGEN  DEPOSITION  VELOCITY  EXPANSION  C 
lOEFFICIENTS.  ) 

WRITE(3,370) { (PV( JJ,KOF)  ,JJ=1,J)  ,K0F=1,KK0F) 

WRITE! 3, 370) ( {HV( JJ,KOF) ,JJ=l, J)  ,KQF*1 ,KKOF) 

READ  NUMBERS  OF  SETS  OF  ORGANS. 

) 

READ! I, 340)NWB,NTHY,N80N,NSKN 
WRITE(3,450) 

WRITE! 3,3 AO) NWB,NTHY,NBON,NSKN 


READ  AND  WRITE  ISOTOPE  DATA 
READ  AND  WRITE  OL * S 

THE  ISOTOPE  CARDS  ARE  FORMATTED  TO  CONTAIN  DATA  TO  BE 
REQUIRED  FOR  LATER  VERSIONS  OF  THIS  PROGRAM.  SUCH  ISOTOPE  CAROS 
NEED  BE  PUNCHED  ONLY  ONCE  TO  COMPILE  A LIBRARY  OF  SAME. 
THEREAFTER,  ONLY  THE  OL'S  NEED  BE  PUNCHED  FOR  EACH  PROBLEM*^ 

WHEN  USING  THE  SEMI-INFINITE  CLOUD  APPXOX I MA TI ON , 

THE  DCF  IS  :(  THE  DOSE  COMMITTMENT  FACTOR  IN  MILLIREM/YR) 
DIVIDED  BY  (THE  ICRP  2 MPCA  TIMES  3.157  Ef07  SEC/YR). 

FOR  EXTERNAL  GAMMA  WHOLE  BODY  AND  EXTERNAL  BETA  SKIN 
DOSE  DCF  a 250  EBAR,  WHERE  EBAR  IS  THE  AVEPAGF  GAmma 
OR  BETA  ENERGY  (mev)  PER  DISINTEGRATION,  RESPECTIVELY.  SEE  AEC 
SAFETY  GUIDES  3 AND  4. 

GAMMA  ENERGIES  AND  ABUNDANCES  ARE  READ  AND  THE  FINITE  CLOUD 
CALCULATION  IS  PERFORMED  IF  NGAMMA  IS  NOT  ZERO. 

REF.  M.E.  MEEK  AND  R.S.  GILBERT  SUMMARY  CF  BETA  AND 
GAMMA  ENERGY  AND  INTENSITY  DATA.  NEnn-12037,  JAN  1R70, 

G.E.  VALLECITOS  NUCLEAR  CENTER,  PLEASANTON,  CALIF. 
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WRITE( 3f460) 

00  200  LL»1,L 

R6AO(l,l60)  IS0(LL ) ,LATWT(LL ) ♦CRTCRG(LL) tOCF ( lL ) , XLMOA ( L L I , 

1 NOEP(LL  LL)  ,NDTR(LL) 

WRITE(3.162)  LL,lSn(LL)  , LATV,  T ( LL  ) .CRTURG  ( LL)  t DCF  ( LL  ) , XL-OA  ( LL  » . 

I N0FP(LL),NGAMMA(LL)  ,NOTR(LLI 

160  F0RMAT(3(1X,A4),F10.1»E10.2,3I5) 

162  FORMAT (1X,I2,IX,3(1X,A4>,F10.1,E10.2,3I5) 

C READ  AND  WRITE  GAMMA  ENERGIES  AND  ABUNDANCES 

IFINGAMMA(LL).EQ.O)  go  to  190 
NG*NGAMMA(LL ) 

READ! 1,170) ( A&AM(LL,IG) ,ABUN(LL,IG) ,IG=1,NG) 

170  FORMAT! 12F6,0) 

WRITE(3,180) ( AGAM(LL,IG) ,ABUN(LL,IG) ,IG*1,NG) 

180  F0RMAT(1X,3X,5( 1X,F9,3,F8.4) ) 

190  READ  ( 1,350)  QL ( LL ) ,OECON ( LL  ) 

200  WRITE(3,290)  ( QL ! LL) ,DECON( LL  ) ) 

REAO(1,202)  ETITLE 
202  FORMAT (20A4) 

WRlTe(3,204)  ETITLE 
204  FORMAT! IHO, 20A4) 

READ  !l,210)  NGAM,nSIG 
210  F0RMAT!20I4) 

WRITE!3,220)  NGAM,NSIG 
220  F0RMAT!lX,*NGAM,NSrG»/lX,I4,lX,I4) 

READ  !1,250)  ( GAMEN ! IN ) , IN*1 ,NGAM ) 

WRITE (3, 2 30)  (GAMEN! IN) , IN*1,NGAM) 

230  FORWATIIX, ’GAMEN  IN  ME V / ( 1 X , 8E 1 2 . 4 ) ) 

READ  11,250)  (SIGZM! IN),IN=1,NSIG) 

WRITE (3, 240)  ( SIGZM! IN) , IN=1,NSIG) 

240  FORMATIIX, 'SIGZM  IN  METERS' /! IX, 8E12. 4) ) 

250  F0RMAT!8E10,4) 

00  260  11=1, NGAM 

260  READ  !1,250)  ! DOSE  I ! 1 1 , J J ) , J J=1 , NSIG ) 

DO  270  11=1, NGAM 

270  WRITE(3,280)  GAMEN ! 1 1 ) , ! DOSE  I ! 1 1 , J J ) , J J=1 , NSI G) 

280  FORMAT! IX, * OOSEI  IN  RAO  METERS**2  PER  CURIE  SEC  AT',F6.3,'  MEV'/ 

1 (lX,eE12.4)) 

290  FORMAT!  IX, 4X, 

1 F13.3,3X,7H  CUR IES,5X,7HDeC0N=  ,E10.3/) 

3.00  F0RMAT!4I5,2F10.2,3E10.2) 

310  F0RMAT!6F10.2) 

320  FORMAT!8F10.0) 

330  FORMAT!3F10.0) 

340  F0RMAT!4I5) 

350  FORMAT!F10.6,E10.3I 
360  FORMAT !E10. 2) 

370  FORMAT!6E10.2) 

380  FORMAT  ! IHO, 30X, ' INPUT  DATA:') 

390  FORMAT! 76H0FAC1LITY,  NO. MONTHS  OF  DATA,  PERIOD,  THERMAL  ENERGY  GEN 
lERATEO  DURING  PERIOD) 

400  FORMATdHO, 'NO.  SECTORS,  STABILITY  CLASSES,  RADII  AND  ISOTOPES,', 

1 • STACK  HEIGHT  IN  METERS,  INVERSION  LID'Z 

2 IX, 'HOLDUP  IN  DAYS, RAINFALL  FRAC TI ON , WASHOUT  COEFF  I C I ^ 

3 • (1/SEC)') 

410  FORMAT! IHO, 'WIND  FREQUENCY  IN  PERCENT  BY  STABILITY  CLASS  FOR', 

1 • EACH  SECTOR  AND  TOTAL  FREQUENCY  FOR  SECTOR’/ 

2 lX,'0IR*,ax,’A',9X,'B',9X,’C',9X,’0',RX,.’E’,RX,’F’, 

3 7X, 'TOTAL') 

420  FORMATdHO, 'WIND  SPEED  IN  METERS  PER  SECOND  BY  STABILITY  CLASS 


119 


no  o o o o o o o r»  n o o r>  o oo 


1 'FOR  FACH  COMPASS  DIRECTION'/ 

2 IX, 'DIR* .RX, 'A'tgX, 'R' ,9X,  'C'.qx, 'D* ,9V,'P' ,9X, 'F* ) 

<^30  FORMAT!  IHI, 'POPULATIONS  IN  SECTOR  SEGMENTS,  TWO  LINES  PER  RADIUS', 

* ' (READ  CLOCKWISE!,  AND  TOTAL  POPULATION  IN  ANNULI'/ 

* IX, 'MM', IX, 'N/ S/TOTAL ',3X, 'NNE /SSH ' , 5X, 'NE/SW',3X, 

* 'ENE/WSW',7X, 'E/W',3X,'ESE/WSW',5X,'SE/NW',3X, 

* 'SSE/NNW'l 

440  FORMAT! IHO, 'DISTANCES  IN  METERS:  MIDPOINT,  LOWER  AND  UPPER  OF', 

* ' BINS'/ 

* IX, 'MM' ,8X, 'K' ,7X, 'XLOW' ,7X, 'XUP' ) 

450  F0RMAT(97H0NUMBER  OF  ISOTOPES  PER  SET  OF  CRITICAL  ORGANS.  (FOUR  SE 
ITS  OF  ORGANS,  TWENTY  ISOTOPES,  MAX.  TOTAL) 

460  FORMAT! IHO, • ISOTOPE  DATA:  ISOTOPE,  CRITICAL  ORGAN,  MRAD/SEC  PER  ', 

* 'CI/METER**3,  decay  CONSTANT  (1/SEC),', 

* 'NOEP,  NGAMMA,  NDTR'/ 

* 1X,T16, •AGAM,ABUN:  gamma  ENERGY  (MEV)  AND  ABUNDANCE  ', 

* ' (NGAMMA  PAIRS) '/ 

* IX, T16, 'ACTIVITY  RELEASED  (CURIES)  DURING  PERIOD,', 

* ' NUCLIDE  DECONTAMINATION  FACTOR'/ 

* IX, 'LL') 


N1*NWB+NTHY+NB0N+NSKN“L 
IF!N1)  470,490,470 
470  WR1TE(3,480) 

480  F0RMAT(1X,'THE  NUMBER  OF  ISOTOPE  CARDS  MUST  MATCH  L.  IF  AN  ISOTOPE 
1 IS  ASSOCIATED  WITH  TWO  ORGANS  IT  IS  TO  BE  COUNTED  TWICE.') 

GO  TO  1551 

490  IF! J-6)500, 520,500 
500  WRITE(3,510) 

510  F0RMAT(1H1,12H  J MUST  BE  6) 

GO  TO  1551 

COMPUTE  VALUES  OF  SIGMZ 

520  00  530  JJ»1,J 
DO  530  MM*l,M  ' 

Y(MM)=X(MM)*0.001 
530  SIGMZ!JJ,MM)=SIGZ(JJ,Y(MM)) 

WRITE!3,540) 

540  FORMAT! IHl, 30X, 'OUTPUT  DATA:'/ 

♦ IHO,  2X, ' SIGMA  ZEE"S  BY  STABILITY  CLASS  AND  RADIUS'/ 

♦ 1X,3X,2H  A,8X,2H  B,BX,2H  C,8X,2H  0 ,8X,2H  E,8X,2H  F/) 

550  WRITE(3,370)( (SIGMZIJJ,MM),JJ=1,J) ,MM=1,M) 

CALCULATE  DECAYED  CHI  OVER  0 IN  EACH  SECTOR  SEGMENT, 

FOR  EACH  STABILITY  CLASS  AND  FOR  ='ACH  ISOTOPE  AND  STORE. 

DIFFUSION  EQUATION  IS  EO.  3.144  IN  SLADE,'  MET.  AND  AT.  ENERGY.' 
(1968)  TlO-24190,  NTIS,  US  DEPT  COMMERCE,  SPRINGFIELD,  VA,  22151 
THIS  PROGRAM  SETS  CHI  OVER  Q TO  ZERO  IF  FREQ  CR  UBAR  ARE  LESS  THAN 
0.01  • 

CONVERT  HOLDUP  FROM  DAYS  TO  SECONDS. 

HOL D*H0L0UP*8.64E4 

CONVERT  ORIGINAL  CURIE  INVENTORY  TO  STACK  RELEASE  INVENTORY 
BY  APPLYING  HOLDUP  AND  DECONTAMINATION. 

DO  560  LL«1,L 

C GO  THROUGH  LIST  BACKWARDS  TO  SIMPLIFY  INGROWTH  CALCULATIONS 

IL»L^1-LL 
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C DECAY  ORIGINAL  INVENTORY 

• 0L( IL )=QL(IL»*EXP<-XLMDA(IL)*H0LD) 

C CHECK  FOR  AND  GROW  IN  DAUGHTERS 

IF  ( NDTR(IL).NE.O)  OL ( I L )=QL ( I L ) *01 ( I L-1 I ♦XLHDA ( I L I / ( XLMDA ( IL ) 
f 1 -XLMDA(IL-I) )*(EXP(-XLMDA(IL-l)*HOLO) 

2 -EXP(-XLMDA( IL» ♦HOLD) ) 

C APPLY  DECONTAMINATION  FACTOR 

0 IF  (DECONdD.LT.l.)  OECONI  I L )*1.0 

0L{ IL)=QL(IL)/DEC0N(IL) 

C SET  RELEASES  OF  LESS  THAN  ONE  MICRQCURIE  TO  ZERO  TO 

• C AVOID  underflows. 

IF  (QL( IL).LT.l.E-6)  OLI IL)*0. 

560  CONTINUE 

561  CONTINUE 
XI=I*1. 0/(16.) 

DO  600  LL«lfL 

00  600  II«ltI 
DO  600  JJsltJ 
DO  600  MMal,M 

1 FI FREO( II, JJ)-0. 01 1580,580, 570 
570  1F(UBAR( 1 1, JJ )-0.01 ) 580, 580, 590 
580  CHIOOIMM, II,JJ,LL)«0.0 

GO  TO  600 

590  IFI SIGMAX.lt. 1. ) SIGMAX=1.E*4 
SIGMA=SIGMZ(JJ,MHI 
IFCSIGMA.GT.SIGMAX)  sigma=sigmax 
ARG=H*H/( 2.*SIGMA*SIGMA) 

IFI ARG.GT.IOO. ) ARG*100. 

OENOM=SIGMA»UBAR( II,JJ)*X|MM) 

DECAY=EXP( (-1. )*(XLMDA(LL) )*X(MM)/IUBAR(II ,JJ)) ) 

CHIOQIMM, II, JJ,LL)=(0.0203)*FRE0(I I, JJ)*EXP(-ARGI*Xr*OECAY/OENO« 
600  CONTINUE 

IFINCOUNT.NE.O)GO  TO  611 

*♦♦**♦♦♦*  WRITE  OUT  CHI/0  SUMMED  OVER  STABILITY  CLASSES  ********** 
WRITEI3,2110) 

WRITE! 3,2120)  I X( MM) ,MM=1,M) 

DO  2020  11=1,1 
DO  2010  MM*1,M 
CHI00S(MM)*0. 

DO  2010  JJ«1,J 
LL*1 

IF(FREOIII,JJ».LE.O..OR.U8ARIII,JJ).LE.O.)  GO  TO  2010 
CHIOOS (MM)=CHIOOS(MM)+CHIOQ( MM,I I , JJ,LL) 

* *EXP| XLMDAILL)*XIMM)/UBAR(I1,JJ)) 

2010  continue 

2020  WRITE(3,2130)  D IR I 1 1 ) , I CHIOQ SI  MM ) ,MM=1 ,M) 


1 

. 1 


» 

il 

Ij 


I 


! 


2110  F0RMAT(1H0,4X, 'CHI/O  SUMMED  OVER  STABILITY  CLASSES*)  I 

2120  FOR'-ATI  IHO,  56X, 'DISTANCE  I ME  TERS ) * /I  X , *01  R • , 12F10.0/ ) ! 

2130  F0RMAT(1X,A4,12(1PE10.2) ) ( 

I 

CALCULATE  DEPOSITION  VELOCITIES  AMD  STORE.  | 

DO  610  II  * I, I ^ 

DO  610  JJ  « 1,J  j 

% HDEPVI 1 1 , JJ )=HV( JJ,1)  ♦ HV(JJ,2)«  UBAR(II,JJ)  ♦ HV(JJ,3)  * UBARI  ' 

1II,JJ)  **  2 * HVIJJ,4)*  UBAR(1I,JJ)  ♦*  3 • 


o 
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610  PDEPVf  1 1 , JJ  )»PV(  JJ,1)  ♦ PV(JJt2)*  UB&Rdl.JJ)  ♦ PVCJJt3J  * UBARf 
llItJJ)  **  2 ♦ PV(JJ,4)*  UBAP(II,Jv.i  *♦  3 

calculate  areas  of  sector  segment:  and  store. 

611  continue 

XXI  = fl.O)  * I 

DO  640 

AREA(MM)= (3,1416)*  (XUP(MM)**2  - XI  OW ( PW) ) / ( XXI I 
C»AP.EA(MM)-1,0 
IF(NC0UNT.NE.0)G0  to  641 
IF(C)  620,620,640 
620  WRITE(3,630) 

630  F0RMAT(51H  ONE  AREA  IS  TOO  SMALL.  'HECK  X,  XUP,AN0  XLOH  DATA.) 

640  CONTINUE 

00  DEPOSITION  AND  CLOUD  DEPLETION. 

DRY  deposition: 

CLOUD  DEPOSITION  AND  DEPLETION  ARE  TREATED  USING  A SIMPLE  MODEL. 
DEPOSITION  STARTS  AT  THE  FIRST  (INNER)  RADIUS.  THE  TOTAL  ACTIVITY 
DEPOSITED  ON  THE  GROUND  IS  CALCULATED  ASSUMING  A UNIFORM  RATE  OF 
DEPOSITION  OVER  THE  ENTIRE  SECTOR  SEGMENT.  CHI  OVER  0 AT  THE 
SECOND  RADIUS  IS  DEPLETED  BY  THE  FRACTION  ON  THE  GROUND  AT  THE 
FIRST  RADIUS.  THIS  PROCESS  OF  DEPOSITION  AND  DEPLETION  CONTINUES 
OUT  TO  THE  LAST  RADIUS  IN  THE  SECTOR. 

THIS  ROUTINE  IS  ESSENTIALLY  THE  SA  E AS  CALLED  FOR  BY  EQUATIONS 
5.44  THRU  5.48  IN  METEOROLOGY  AND  ATOMIC  ENERGY  (TID-24190), 

EXCEPT  HERE  THE  DEPOSITION  AND  DEPLETION  DO  NOT  VARY  CONTINUOUSLY 
WITH  DISTANCE  AS  DOES  EO.  5.48. 

CONTINUE 

WET  DEPOSITION: 

WET  DEPOSITION  IS  TREATED  USING  THE  SECTQO  AVERAGED  ANALOG  OF 
EQUATIONS  5.63  TO  5.65  IN  MET.  AND  AT.  ENERGY  1968,  THIS 
implies  that  wet  DEPOSITION  IS  PROPORTIONAL  TO  THE  ACTIVITY 
PER  SQUARE  METER  ABOVE  THE  GROUND,  THE  PROPORTIONALITY 
CONSTANT  BEING  CALLED  WASHCO  HEREIN'.  SINCE  WET  DEPOSITION 
IS  A FIRST  ORDER  PROCESS,  CLOUD  DEPLETION  IS  SIMPLY 
EXPONENTIAL.  WET  DEPOSITION  IS  TREATED  INDEPENDENTLY  OF  ALL  DOSE 
CALCULATIONS.  ALTHOUGH  THE  CLOUD  d DEPLETED  BY  WASHOUT  DURING 
THE  WASHOUT  CALCULATIONS,  STORED  CH)  OVER  C'S  ARE  NOT.  IF  AN 
UNNORMALIZED  WINOROSE  IS  USED,  I.E.  THE  SUM  OF  THE  WIND  FREQUEN- 
CIES IS  MUCH  LESS  THAN  1.0,  RAINF  SHOULD  RE  ENTERED  AS  l.E+00.  IF 
A NORMALIZED  ANNUAL  WIND  ROSE  IS  USED,  I.E.  THE  SUM  OF  THE  WIND 
FREQUENCIES  IS  1.0,  RAINF  SHOULD  BE  THE  DECIMAL  FRACTION  (NOT 
PERCENT)  OF  THE  TIME  it  RAINS.  DEFAULT  VALUES  FOR  RAINF  AND  WASHCO 
ARE  5.E-02  AND  2.E-04  1/SEC  RESPECUVELY  FOR  A MODERATELY  WET 
CLIMATE  AND  PARTICLE  DIAMETERS  OF  bcVERAL  MICRONS.  SEE  MET  AND  AT 
E FOR  DETAILS. 

641  CONTINUE 
DO  900  LL»1,L 

lF(riDEP(LL).LE.O)  GO  TO  900 

IF(NDTR(LL)  .EQ.O.AND.QLILD.EQ.O.)  10  TO  900 

IF(NDTR(LL).E0.1.ANO.OL(LL).EO.O..ANO.OL(LL-1).EO.O.)  go  to  900 
650  C5  * 0.0 
C8*0. 

DO  BOO  II«1,I 
DO  740  JJ>1,J 


122 


oooo  o n n r»  r>  o 


DO  740 
FID  * 1.0 

C CHECK  FOR  CHIOQ  EQUALS  ZERO 

IFI FREOI I I, JJ),GT.O..AND.UBAR( II ,JJ) .GT.O. ) GO  TO  670 

0G(MM«JJ)«0.0 

W0G(MM,JJ|«0.0 

F(MM+1)«0. 

GO  TO  740 
670  K1  « MM 
C*2. 

C SINCE  CHIOQ  INCLUDES  DECAY  THE  DAUGHTER  OS'S  ARE  GROWN  BACK  TO  THE 
C . DAUGHTER  INGROWTH  IS  ADDED  TO  THE  INPUT  DAUGHTER  STRENGTH. 
QS»tLILL  ) 

IF  (NDTRILL ) .NE.O)  OS=OS+OL ( LL-1 ) * XLMDA I LL I 

1 /( XLMDA(LL)-XLMDA(LL-in 

2 ♦(EXP( (XLHDAILL)-XLMDA(LL-1))*X(MM» 

3 /UBARf 1 1 , JJ) )-1.0l 

C PICK  UP  FACTORS  I - F(2)  -F(3)- -F(MM) 

DO  680  MM1>1,K1 
680  C = C-FIMMD 

c the  following  conserves  (Q*FREQ). 

IF(C)  690,690,700 
600  OGIMM, JJI^O.O 

CHIOQIMM, II , JJ,LL)«0.0 
MPLSl»MM+l 
F(MPLSD»2.0 
GO  TO  740 

C THE  FOLLOWING  STEP  REPLACES  THE  OLD  CHI  OVER  Q 

C WITH  THE  DEPLETED  ONE 

700  CHIOQIMM, II, JJ,LL)  = C ♦ CHIOOl MM, 1 1 , JJ,LL) 

IFINDEPILL) .GE.2)  GO  TO  720 
C DO  PARTICULATE  DEPOSITION  AND  DEPLETION. 

710  OG(MM,JJ)=CHIOOIMM,II,JJ,LL)*OS*AREA(MMJ*PDEPVIII,JJI 
GO  TO  730 

C DO  HALOGEN  DEPOSITION  AND  DEPLETION 

720  OGIMM, JJ)=CHI00IMM, II, JJ,LL)*0S*AREA|MM)*HDEPV|II ,JJ) 

730  MPLSl*  MM+1 

0IJ=0S*FREQ( II,JJ»*O.Ol 
FIMPLSD*  QGIMM, JJ)/10IJ) 

CALCULATE  WET  DEPOSITION 

IFIWASHCO.LE.O..OR.RAINF.LE.O. ) GO  TO  740 
DECAY=EXP|-XLMDAILL  »*X(MM)/UBAR| II , JJ) ) 

DENOM  = UBARI 1 1 , J J ) ♦X I MM ) *2*3. 141 6/1 

WQGIMM, JJ )= IRAINF*WASHC0*FRE0I II ,JJ)/DENOM)*OS*AREAIMM)*OECAY*0.01 
**EXPI-WASHCO*XIMM)/UBAR( II, JJ) ) 

ILAST  TERM  INCLUDES  CLOUD  DEPLETION  DURING  WET  DEPOSITION.) 

740  CONTINUE 

SUM  CURIES  ON  THE  GROUND  OVFR  ALL  STABILITY  CLASSES  IN 
EACH  SECTOR  SEGMENT  AND  TOTAL  CURIES  ON  THE  GROUND. 

750  C1»0.0 
C7«0. 

DO  790  MM*l,M 
C2»0.0 
C6«0. 
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DO  760  JJ*  ItJ 

IF(wOG(MM,JJ).LE.(  l.F-12)  ) WQG  ( J J ) *0. 

1F(CG(MM,  JJ).LE.(1,E-12»  ) QG(M*',  JJ)*0.0 
C6=C6+W0G(MM,JJ) 

760  C2  = C2  ♦ 0G(MM,JJI 
WOGMI (MM,  II )*C6 
OGMKMM.in  « C2 

COMPUTE  SATURATED  ACTIVITY  ON  THE  GROUND  FOR  RELEASE  PERIOD  OR 
ONE  YEAR  ( IN  PICOCURIES  PER  M**2  • 

C4=2.6E+06 

C*M0NTHS 

IFtC.LE.O.)  C=12. 

DEC AY=XLMDA(LL) *0*04 

IF( DECAY. LT.I.E-3)  0ECAY=1.E-3 

IF(DECAY-5. ) 780,780,770 

770  OGAMKMM,!!)  = (C2  ♦ 1.0  E + 12 )/(  AREA  ( MMJ  ♦DECAY) 

WQGAMI (MM, 1 1 ) = (C6  ♦ l.O  E+ 12 ) / ( ARE A ( MM) ♦DEC AY) 

781  ACTDEN(MM,  1 1 , LL  ) =OGAM  I ( MM,  I I ) + kiCGAMl  ( MM,  1 1 ) 

GO  TO  785 

780  OGAMKMM, II)  =((C2  ♦ leO  E + 1 2 ) / ( ARE  A ( MM)  *DEC  AY)  ) ♦ ( 1 . 0-EXP  ( -DECAY  ) ) 
WOGAMI (MM,I  I ) *((C6  ♦ 1 . OE-H  2 ) / ( ARE A ( MM) ♦DEC AY) ) ♦ ( 1 . 0-EXP ( -DECAY ) ) 

782  ACTDEN(MM,II,LL)=OGAMI(MM,II )♦ WOGAMI (MM,II) 

785  C7=C7+C6 

790  Cl  = Cl  ♦ C2 
C8*C8+C7 
C5  = C5  ♦ Cl 
WT0G«C8 
800  TOG  « C5 


SUM  CURIES  DEPOSITED  WITHIN  RADII. 

C3«0. 

C4*0. 

DO  810  MM*l,M 
DO  810  11*1,1 

IF( WOGMI (MM, II ) .LE.( l.E-10) ) WOGMI ( MM , 1 1 ) =0. 0 
IF(CGMI (MM, 1 1 ) .LE.( l.c-10) ) OGMI ( MM , 1 1 ) =0. 0 
IF( ACTDEN(MM, I I,LL).LE. ( l.E-10) ) ACTDEN(MM,II ,LL)*0.0 
C4=C440GMI(MM,II) 

C3*C3+WQGMI (MM,I1)  ' 

WTOGMCMM)  «C3 
810  T0M(MM)*C4 

IF(NCOUNT.NE.O)GO  TO  900 

FOR  PARTICULATES  AND  HALOGENS 

WRITE  OUT  CURIES  PER  SOUARE  METER  ON  THE  GROUND, 

SATURATED  ACTIVITY  PER  M^^2.  ^ 

820  WRITE(3,830) { ISO(LL) ,LATWT(LL)) 

9 830  FORMAT! IHl, 20X, 'ACTIVITY  DENSITY  AT  END  OF  RELEASE  PERIOD  FOR  *, 

♦ 'A  UNIFORM  RELEASE  RATE*/ 

♦ 40X, *PIC0CURIES/M^^2  OF  *,44,44) 

C WRITE(3,839) 

839  F0RmaT(44X, *DRY  DEPOSITION*//) 

WRITE(3,840) 

O 840  F0RMAT(47X,7H  SECTOR , //I X,9H  DI STANCE , /I  X ,9H  ( METER S ),/ 16X , 2H  N 

1,7X,4H  NNE,6X,3H  NE,6X,4H  ENE,8X,2H  E,8X,4H  ESE,6X,3H  SE,6X,4H  SSE 


o 
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2,4X,12 

DO  860 
850  FORMAT 
860  WRITE( 
WRITE( 
870  FORMAT 


23H  NW, 
DO  880 
) WRITE( 
QDISP= 
WRITEI 
> FORMAT 


H AREAS(M**2),/95X,19H  OF  SECTOR-SEGMENTS,/) 
(9F10.0,5X,E10.3) 

3,850)  (X(MW),(0GAMI(MM,in,II=l,8),AREA(MMn 
3,870) 

(1H0//17X, 

?H  S,6X,AH  SSH,7X,3H  SW,6X,4H  WSW,7X,2H  W,7X,4H  WNW,8X, 
6X,4H  NNW,6X,17H  CURIES  DE POS I TED, /97X , 13H  WITHIN  RADII,/) 

3,8  50  ) (X(MM) , (OGAMI (MM,I I ) , 1 1 =9 , 1 ) , TCMI MM) ) 
QL(LL)*TOTFRO/lOO. 

3,890) (ISO (LL) ,LATWT<LL) ,XUP(M) ,TQG,QDISP) 

(IHO/ 

IX, lOX, 'TOTAL  CURIES  OF ' , A4, A4 ,4X , 'HE POS I TED  ON  THE*, 

' GROUND  WITHIN* ,F10.0, • METERS  I5',F10.3,*  CURIES*/ 
1X,10X,'0UT  OF  A TOTAL  0F',F10.3,'  CURIES  DISPERSED*, 

• DURING  THIS  PERIOD. *) 


ACCORDING  TO  ONE  REFERENCE  (H.P.  18 ,1 ,P,74, JAN  70) 

0.09  CURIE  PER  LITER  APPEARS  IN  COWS  MILK  PER  CURIE  PER 
SQUARE  METER  OF  I 131  DEPOSITED  ON  THE  GROUND. 


WRITE  WET  DEPOSITION 


1F(RAINF.LE 
WRITEI 3, 830 
WRITE{3,891 

891  F0RMAT(44X, 

* 'RAI 

WRITE(3,840 
DO  892  MM=1 

892  WRITE(3,850 
WRITEI 3, 870 
DO  893  MM=1 

893  WRIT£I3,850 
C=OLILL )*RA 
WRITEI 3, 894 

894  FORMATIIHO/ 

1 IX, 

2 

3 IX, 

4 

900  CONTINUE 


.O..DR.WASHCO.LE.O.)  GO  TO  900 
) I ISOILL)  .LATWTILD) 

) (WASHCO,RAINF) 

•WET  0EP0SITI0N'//20X, 'WASHOUT  CHEF  *' ♦1PE10.2,5X, 
NFALL  FREQUENCY  =',E10.2//) 

) 

,M 

)(XIMM),IW0GAMIIMM,II ) ,11=1,8) , AREA  I MM)) 

) 

,M 

) I XIMM) , I WOGAMI IMM,II ) ,11=9,1 ) ,WTQGM(MM) ) 
INF*TOTFRO/100. 

) I ISOILL) ,LATWTILL) ,XUPIM) ,WTQG,C) 


lOX, 'TOTAL  CURIES  OF • ,A4,A4 ,4X , 'nPPOSITED  ON  THE*, 

* GROUND  WITHIN* ,F10.0, ' METERS  IS',F10.3,*  CURIES'/ 
10X,'OUT  OF  A TOTAL  0F',F10.3,'  CURIES  DISPERSED', 

' DURING  rainfall.*) 


CALCULATE  ALL  IMXJXIXL)  DOSES. 

SUM  DIFFUSED,  DECAYED  AND  DEPLETED  DOSES  OVER  ALL  STABILITY  CLASSE 
PRESERVING  ISOTOPE  DESIGNATION. 


CALCULATE  EXTERNAL  GAMMA  WHOLE  BODY  DOSES  USING  THE  DOSE  INTEGRALS. 


DO  1050  LL«1,L 
DO  1050  11*1,1 
DO  1050  MM*1,M 
WIDTH*2.0*PI*XIMM)/I 
C2*0. 

DO  1040  JJ»1,J 
DMIJ*0. 

IFI FREQI I I, JJ ).LE.0..0R.UBAR 1 1 1 , JJ) .LE.O. ) GO  TO  1040 
IFINDTRILL) .EO.O.AND.OLILL) .EO.O. ) GO  TO  1040 


o 


125 


oooonoo  r>  nnoo 


SINCE  CHIOQ  INCLUDES  DECAY  THE  DAUGHTER  OS'S  ARE  GROWN  BACK  TO  THE 
STACK.  DAUGHTER  INGROWTH  IS  ADDED  TO  THE  INPUT  DAUGHTER  STRENGTH. 


QS*CL(LL) 

IF(NDTR(LL)  .NE.O)  OS=OS+OL  ( LL-1  ) *XL*'DA  ( LL  ) / ( XLMOA  ( LL  ) -XLMDA  t LL-1 ) ) 

1 *(EXP{ ( XLMDA(LL)-XLH0A(LL-1» )*X(MM) 

2 /UBAR( II, JJ) )-1.0) 

IF  (LL.GT.NWB.OR.NGAHMA(LL ) ,EQ.O)  GO  TO  1020 
DECAY=EXP(-XLMDA(LL)*X(HM)/UBAR(  II , JJ) ) 

SI6^^A=SIGMZ  ( JJ,MM) 

IF  (SIGMA.GT.SIGZMI 2) ) GO  TO  910 
IS*1 

IF{SIGHA.LT.SIGZM( 1) ) SIGHA=SIGZH( 1 ) 

GO  TO  950 

910  IF  (SIGMA. LE. SIGZM(NSIG-I) ) GO  TO  920 
IS=NSIG-1 

IF ( SIGMA. GT,SIGZM(NSIG) ) SIGMA=SIGZM ( NSI G) 

GO  TO  950 

920  IS»8.0*ALOG10( SIGMA) 

DO  930  IN*IS,NSIG 

IF  (SIGZM(IN+1).GE. SIGMA)  GO  TO  9A0 
930  CONTINUE 
940  IS* IN 

9 50  PS* (SIGMA-SI GZM( I S > ) /( S IGZM( I S+1 ) -SI GZM(  I S) ) 

DISUM*0. 

NGAML=NGAMMA(LL ) 

DO  1010  NG*1,NGAML 
EGAM=AGAM(LL,NG) 

INTERPOLATE  DOSE  INTEGRAL  TABLE  TO  OBTAIN  01 

IF  (EGAM.GT.GAMEN(2))  GO  TO  960 

IG=1 

GO  TO  1000 

960  IF  (EGAM.LE.GAMEN(NGAM-1 ) ) GO  TO  970 
IG=NGAM-1 
GO  TO  1000 

970  IG*R.0*AL0G10( 100.0*EGAM) 

DO  980  IN=IG,NGAM 

IF  (GAMENdN  + D.GE.EGAM)  GO  TO  990 
980  CONTINUE 
990  IG»IN 

1000  PG* (EGAM-GAMEN( IG) ) / ( GAMEN ( I G+ 1 > -GAMEN ( I G ) ) 

DI=  DDSEI ( IG, IS>*( 1.0-PG)*( 1.0-PS)*DCSEI ( I G, I S+1 ) ♦ ( 1. 0-PG) *PS 
1 +OOSEI  ( IG  + 1,  IS)*PG*(  1.0-PS)  + DCSEI  ( IG+1 , IS-H)*PG*PS 
1010  DISUM=DISUM+OI*ABUN(LL,NG) 

DMI  J*DISUM*OECAY*OS*FREO( 1 1 , J J ) *10. 0/ ( WI DTH^UBAR ( 1 1 , J J ) ) 

GO  TO  1030 

10  20  DM IJ= CHI  00 (MM, 1 1 , J J ,LL ) ♦OS*OCF ( LL ) 

1030  0MIJL(MM, II , JJ,LL)=0MI J 

IFIDMIJ.LE. (l.E-06) ) 0MIJ=0. 

1040  C2=C24-DMIJ 
1050  DMIL(MM,II,LL)*C2 

FOR  ORGANS  FOR  WHICH  THERE  IS  INPUT  DATA: 


SUM  INDIVIDUAL  DOSES  IN  EACH  SECTOR-SFGMENT. 

CALCULATE  ORGAN  MAN-REms  IN  EACH  SECTOR-SEGMENT, 

SUM  AT  RADII  , THEN  TOTAL  FOR  EACH  ORGAN. 


o o o o r»  r* 


PRIM  MILLIREM  BY  SECTOR-SEGMENT* 

PRINT  MAN-REM  BY  SECTOR-SEGMENT* 

PRINT  MAN-REM  BY  RADIUS. 

WRITE  OUT  ISOTOPE  DATA 


IF(NWB)1200, 1200*1060 
1060  JlsNWB 
C6«0.0 

DO  1090  MM«1*M 
DO  1090  II«  I* I 
C1*0.0 

DO  1070  LL*1*J1 

IF(0MIL(MM, II.LL) .LF.( l.E-06)  ) DMI L ( MM* 1 1 *LL I »0. 
C1=C1+0MIL(MM*II*LL) 

1070  WBDOSIMM* II  )«C1 

IF(NCOUNT.NE.O)GO  TO  1081 

1080  DPWBIMM,  I I )=WBDC*S<MM*  1 1 ) ♦POP  ( MM*  1 1 ) ♦ ( 0. 001 1 
GO  TO  1082 

1081  DPWBIMM, I I)=WBOOS(MM*I I )^P0PUL(MM,I1 )*(.001) 

1082  CONTINUE 
C6=C6+DPHB(MM* III 

1090  0PWBM(MM)*C6 
WRITE(3,1100) 

1100  F0RMAT(1H1///25X*82H  DOSE  TO  AN  INDIVIDUAL  IN  THE  INDICATED  CO 

IMPASS  SECTOR  AND  ANNULAR  RING  ( HREM) //4X ,5H  0IST*4X*2H  N* 

24X*4H  NNE* 

24X,3H  NE*3X,4H  ENE*4X,2H  E,4X,4H  ESE,4X,3H  SE,3X,4H  SSE*4X,2H  S*4X 
3»4H  SSW*4X,3H  SW*3X*4H  WSW*4X*2H  W*4X,4H  WNW*4X*3H  NW*3X*4H  NNW/I 
DO  1110  MM*1,M 

1110  WRITE(3»1120){X{MM)*(WBD0S(MM*II )*II»1*II) 

1120  F0RMAT<F9.0,4X,17E7.2) 

WRITE(3, 1130) 

1130  FORMAT( 1H0//15X, 'POPULATION  DOSE  IN  THE  INDICATED  COMPASS*. 

1 ' SECTOR  AND  ANNULAR  RING  ( MAN-REM )' /31 X AVERAGED* * 

2 * TO  CONFORM  TO  THE  POPULATION  WHEEL  WIND  SECTORS  * /41X * 

3 • FOR  THE  ESTIMATED  POPULATION* ,//122X* 8H  PUNNI NG* /4X * 5H  DIST.4X* 

4 2H  N*4X,4H  NNE* 

54X,3H  NE*3X*4H  ENE,4X,2H  E,4X,4H  ESE,4X*3H  SE.3X,4H  SSE,4X,2H  S*4X 
6*4H  SSW*4X*3H  SW.3X.4H  WSW.4X.2H  W*4X*4H  WNWt4X,3H  NW.3X.4H  NNM* 

7 3X.6H  TOTAL) 

DO  1140  MM*1*M 

1140  WRITE (3* 1120) ( XUP ( MM ) * ( OPWB ( MM , 1 1 ) *II»1*I ) *DPWBM(MH) I 
WRITE! 3* 1150) 

1150  FORMAT! IHO/ 

1 30X,35H  THE  CONTRIBUTING  RADIONUCLIDES  ARE*/) 

WRITE! 3*1160) 

1160  FORMAT! lOX, 8H  I SOTOPE, 5X , 15H  CRITICAL  ORGAN, 6X.4H  DCF*6X,16H  CURIE 
IS  RELEASED, 10X.15H  DECAY  CONSTANT.lOX* 'DECON  FACTOR* */ *83X*8H  !1/S 
2EC)  ) 

C2=0. 

DO  1170  LL«1*NWB 
C2*C2+0L!LL) 

1170  WRITE! 3,1180) ! I SOI LL ) *LATWT ! LL ) .CRTORG ! LL ) * DCF ! LL ) * OL ! LL ) * XLMDA! LL 
1 ) ,DECON!LL)  ) 

WRITE!3*1190)  C2 

1180  FORMAT! 10X,A4, IX, A4*13X*A4*4X *F 10.2* 7X*E 11.3, 15X, El 0.2, lOX.GlO. 3) 
1190  F0RMAT!48X* 'TOTAL  Cl.  **011.3) 

C 

1200  IF!*.‘THY)1280*1280*1210 


o 


1210  J2«r>IWB«l 
C7*0.0 

DO  1240  MM«1,M 
DO  1240  II-l.I 
J3=NWB+NTHY 
C2*0,0 

DO  1220  LL«J2tJ3 

IF{DMlLtMH,  ll.LD.LE.f  l.E-06))  DMI L ( , 1 1 ,LL) -O. 

C2=C2+DM1L(MM,II,LLI 
1220  THYD0S(MM,II)«C2 

1230  0PTHY(MM, II )=THYDOSIMM, 1 1 )*POP< PM, 1 1 ) ♦< 0. 001 ) 

C7*C7<-DPTHY(MM,in 
1240  DPTHYM(MM)=C7 
WRITE(3,1100) 

DO  1250  MM*1,H 

1250  WRITE{3,1120)(X(MP),(THYD0S(PP,II) ,11*1,1)1 
WRITE(3, 1130) 

DO  1260  MM»1,M 

1260  WRITE! 3, 1120) ( XUP( MM ) , ( DPTHY ( PM, 1 1 ) ,11  = 1,1) ,OPTHYMIHM) ) 
WRITE(3,1150) 

WRITE!3, 1160) 

C2»0. 

DO  1270  LL=J2,J3 
C2*C2*QL(LL) 

1270  WRITE < 3, 1180) ( ISO! LL),LATWT(LL) ,CRTORG(LL) ,nCF ! LL) , OL ! LL ) ,XLM0A(LL 
1 ) ,0EC0N(LL)  ) 

WRITE!3,1190)  C2 
1280  IF!NB0f4)1350, 1360, 1290 
1290  J4=NWB^NTHY+l 
C8*0.0 

DO  1320  MM»1,M 
DO  1320  II»1,I 
J5=MWB+MTHY+NB0N 
C3»0.0 

DO  1300  LL»J4,J5 

IF!DMIL!MM,  I I.LD.LE.!  l.E-06))  DPI  L!  MM,  1 1 ,LL)-0. 

C3=C3+0MIL!MM, II,LL) 

1300  BONDOS!MM,I I )=C3 

1310  OPBC‘^!MM,  II  ) = B0ND0S!PM,II)*P0P!MM,II)*!0.001) 

C8=C8*DPB0N!MM,II) 

1320  0PB0NM!MM)»C8 
WRITE!3,1100) 

DO  1330  MM«1,M 

1330  WRITE!3,1120)!X!MM),!B0ND0S!PM,1I) ,11=1,1)) 

WRITE!3,1130) 

DO  1340  MM*1,M 

1340  WRITE! 3, 1 120 ) ! XUP ! MM ) , ! DPBON ! MP, 1 1 ) , 1 1 =l , I ) ,DPBONMIPM) ) 
WRITE!3,1150) 

WRITE!3,1160) 

C2-0. 

DO  1350  LL«J4,J5 
C2=C2+0LILL  ) 

1350  WRITE!3, 1180) ! I SO! LL) , LATWT ! LL ) ,CRTORG! LL ) ,nCF ! LL) , OL ! LL) , XLMDA! LL 
1),0EC0N!LL)  ) 

WRITEI3, 1190)  C2 
1360  IF!MSKN)  1440,1440,1370 
1370  J6*NWB>NTHY+NB0N+1 
C9»0.0 

DO  1400  MM«1,M 
DO  1400  11-1,1 
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JTsNWB+NTHY+NBON+NSKN 
C4=0.0 

DO  1380  LL*J6,J7 

IF(DMIL(MM, II,LL).LE.( l.E-Oe) ) DMI L ( KP , 1 1 ,LL) *0. 

C4=C4+0MIL(«M, II.LL) 

1380  SKN00S(MM,  I n*C4 

1390  DPSKN(MNi,  II  ) =SKNOOS(  MM,  1 1 )*POP(MM,II  )*I 0.001) 

C9=C9+DPSKN(MM,II) 
l<f00  DPSKNM(MM)*C9 
HRITE(3,1100) 

DO  1410  MM»1,M 

1410  WRITE (3, 1120) <X( MM), (SKNOOSIMM,! I ) ,11*1,1)) 

WRITE(3,1130) 

DO  1420  MM*1,M 

1420  WRITE(3t 1120) (XUP(MM),(DPSKN(MM,II ) , 1 1 *1 , 1 ) ,DPSKNM( MM) ) 
WRITE(3,1150) 

WRITE(3,1160) 

C2*0. 

DO  1430  LL»J6,J7 
C2=C2+OL(LL) 

1430  WRITE(3,1180) ( ISO(LL)tLATWT(LL) ,CRTORG(LL) ,DCF(LL) , OL ( LL ) , XLMOACLL 
D.DECONtLL)  ) 

WRITE  (3.1190)  C2 
1440  CONTINUE 

C A PATCH  TO  WRITE  INDIVIDUAL  DOSES  MAY  BE  INSERTED  HERE 

Cl  = 0.0 
DO  1450  MM»1,M 
DO  1450  II»1,I 

^ Cl  = Cl  ♦ POP(MM,II) 

) 1450  T0TP0P=C1 

i Cl  = 0.0 

: • HRITE(3.1460)  FAC IL .MONTHl ,NYR1 , 

[ ♦ M0NTH2,NYR2,ETHERM, months, TOTFRO.TOTPOP, HOLDUP, H 

1460  F0RMAT(1H1/46X, 'PROBLEM  SUMMARY • //6X FACl LI TY • ,14X ,' PERI OD* , 

[ % * IIX, 'ENERGY' ,4X, 'MONTHS  OF ' , 2 X , ' TOTAL  * ,5X , 'TOT AL * , 7X , 

I * 'HOLDUP* ,2X, 'HEFF '/24X, 'FROM' ,7X, 'TO* ,6X, ' (MW-D(TH) ) ', 

! * 2X, 'OPERATION* ,2X, 'FREQUENCY' ,2X, 'POPULATION* ,2X, • (DAYS) • , 

• • 2X, ' (METERS) '/ 

* 1X,5A4,2( 1X,A4, 15) ,F11.0,1X,F5.2,6X,F8.1,3X,F11.0,1X,F6.1, 

! * 2X,F7.0) 

• WRITE(3, 1520) 

1520  F0RMAT(1H0,51H  RADIONUCLIDE  CONTRIBUTIONS  TO  THE  TOTAL  POPULATION 
1 DOSE  ARE,/) 

• IF(NCOUNT.NE.O)GO  TO  1541 

DO  1543  LL*=1,L 

I 1 DPL=0. 

' • DO  1530  II»1,I 

DO  1530  MM*1,M 

I 1530  DPL=0PL+0MIL(MM, II,LL)*POP(MM,II )*(0.001) 

( m IF(DPL.LE.l.E-06)nPL=0.0 

WRITE(3, 1540) ISO(LL) ,LATWT(LL),DPL 

1540  F0RMAT(1X,A4,1X,A^^,14X,E10.5,8H  MAN  REM) 

A 1543  CONTINUE 

GO  TO  1550 

1541  DO  1549  LL«1,L 
DPL=0. 

DO  1542  II«1,I 
DO  1542  MM*1,M 

1542  DPL  = DPL  + OMIL(MM, Il,LL)*POPULrMM,II  )*(,001) 
WRITE{3,1540)IS0(LL),LATwT(LL) ,OPL 


o 


129 


1549  CONTINUE 

1550  CONTINUE 

C COMPUTE  TOTAL  DOSE  PER  SECTOR 
WRITE(3,3210) 

3210  FORMAT!  • 1V///40X, 'TOTAL  INTEGRATED  DOSE  DUE  TO  GROUND  't 
I'OEPOSITION'//) 

DO  2001  LL-l.L 

REAr( 1, 2005 )DOSFAC,XLAMBO (LL) f TIME, NORGAN 
2005  F0RMAT(F10.4,E15.3,F5.1,I5) 

DO  3220  MM*1,M 
DO  3230  II«1,I 

0R(MM, 1 1 )=ACTDEN(MM,II ,LL)*( 1.0E-12)*D0SFAC*24, 

TD<MM, 1 1 ) = (DR(MM, II J /XLAMBD ( LL ) ♦ ( 1 .-E XP ( -XL AMBD ( LL ) ♦TI ME ) ) )*1000. 
IF( DR(MM, II ) ,LE.1.E-15)DR(MM,II )*0. 

IFITDIMM,  in.LE.l.E-15)TD!MM,II)*0. 

3230  CONTINUE 
3220  CONTINUE 

WRITE! 3,2030) ! I SO! LL ) , LATWT! LL ) , TI ME ) 

2030  FORMAT!'  ',14X, 'TOTAL  INTEGRATED  DOSE  TO  AN  I NDIV lOUAL ! MREM ) ', 
♦'DUE  TO  GROUND', 

1'  DEPOSITION  OF  RADIONUCLIDE  ',A4,A4/ 

250X, 'TOTAL  TIME  IN  DAYS*  ',F5.1/) 

IF!N0RGAN.E0.1)G0  TO  2033 

2031  WRITE!3,2032) 

2032  F0RMAT!42X, 'THIS  CALCULATION  APPLIES  FOR  WHOLE  BODY  DOSES'//) 

GO  TO  3200 

2033  WRITE!3,2034) 

2034  FORMAT!52X, 'SKIN  DOSE  RECEIVED') 

3200  WR1TE!3,840) 

DO  2000  MM*1,M 

WRITE! 3, 6700) (XI MM), fTDf MM, 1 1 ) ,1 1=1 ,8) ,AREA(MM) ) 

6700  FORMAT!F10.0,2X,8F10.3,5X,E10.3) 

2000  CONTINUE 

3300  FORMAT! 1H0//17X,2H  S,6X,4H  SSW,7X,3H  SW,6X,4H  WSW,7X, 

12H  W,7X,4H  WNW,8X,3H  NW,6X,4H  NNW,6X,/) 

WRITE!3,3300) 

DO  3310  MM»1,M 

3310  WRITE!  3,6700) ! X! MM ) , ! TO! MM, 1 1 ) , I I *9 , 1 ) ) 

IFINQRGAN.EO.DGO  to  2001 

C COMPUTE  THE  MAN-REM  TOTAL  FOR  WHOLE  BODY  CALCULATIONS 

XMNREM«0. 

XMANRM«0, 

DO  2050  MM«l,M 
DO  2050  11=1,1 
IF!NCOUNT.GT.O)GO  TO  4001 
XMNREM=POP! MM, II ) ♦TD ! MM , I I ) ♦«001 
IF! XMNREM.LE.1.E-10)XMNREM=0. 

XMANRM=XMNREM'*-XMANRM 
GO  TO  2050 

4001  DO  4002  MM*1,M 
DO  4002  11=1,1 

4002  FX2!MM, 1 1 ) = AP  + ! 1.-AP)^EXP!-XLAM^TIME-TLAG*TAWL) 

IF|FX2!MM,I I ).LE.1.E-06)FX2!MM,II)*0. 

DO  <*003  MM«1,M 
DO  4003  11=1,1 

4003  P0P2!MM,  I I)=P0P!MM,II)*FX2!Ml«,n  ) 

IF!P0P2!MM, II ),LE..9)P0P2!MM,II )»0, 

DO  4004  MM*1,M 

DO  4004  1 1=  1,  I 

XMNREMxP0P2!MM, II )^T0( MM,II )^.001+! ! PCPUL ! MM , i n _P0P2 ! MM, I I ) ) 


1 

* /2.)*TD(MM, II )*.001 
IF( XMNREM.LE.1.E-10|XMMREM*0. 

XMA*JRM*XMNREM<-XMANR»* 

400A  CONTINUE 
2050  CONTINUE 

HRITE(3.2060)(XMANRM| 

2060  FORMATC  •/ /AOX TOTAL  MANREK  FOR  THIS  NUCLIDE*  •♦F15,6//I 
2001  CONTINUE 

IFINCOUNT.GT.OIGO  TO  3050 
REAC< lt2073) ( AP , XLAM,TLAG , TA WL ,NEVAC ) 

I F(NEVAC) 2071, 3050,2140 

2071  CONTINUE 
2073  FOR''AT(4F10.2,I5I 

WRITE(3,3060) 

3060  FOR‘^AT(  • l'///55X, ‘EVACUATION  MODEL*/// 

141X, 'POPULATION  MOVED  BEYOND  FURTHEST  RADIUS  SHOWN*//) 

WRITE! 3, 2072) ( AP , XL AM , TL AG, TAWL ) 

2072  FORMAT!*  *, 20X, * FRACTION  OF  POPULATION  UNAFFECTED  BY*, 

IIX, 'EVACUATION*  *,F10.2/ 

220X, 'EVACUATION  RATE*  *,F10.2/ 

320X,*TIME  LAG  TO  EVACUATION  NOTICE*  *,F10.2/ 

420X,*TIME  BETWEEN  IMPENDING  CORE  MELT  AND  LEAKAGE*  *,F10.2//) 
WRITE!3,2112) 

2112  FORMAT!*  ',* CLOUD  VELOC I Tl £ S ! M/SEC )* ) 

DO  Z080  11*1,1 
CL0VEL*0. 

DO  2070  JJ*1,J 

CLDV*UBAR!II,JJ)*FREO( II,JJ) 

2070  CLDVEL*CLDV'-CLDVEL 
CL0VEL*CLDVEL/100, 

WRITE!3,2104)CLDVEL 

2104  FOR*'AT('  ',E15.7) 

DO  2077  MM»1,M 

IF!CLDVEL.LE.1.0E-05)G0  TO  2076 

2075  CLDTIM!MM,II)*X!MM)/CLDVEL 
GO  TO  2077 

2076  CLDTIM!MM,II)*1.E  10 

2077  CONTINUE 
2080  CONTINUE 

WRITE! 3,2113) 

2113  FORMAT!*  *, 'CLOUD  TRAVEL  TIME(SEC)*) 

DO  2106  MM»1,M 

2106  WRITE! 3,2105) !CLDTIMIMM,  II ) ,11*1,1 ) 

2105  FORMAT!*  *,16E7,2) 

DO  2090  MM*1,M 
DO  2100  11*1,1 

FX! MM, 1 1 )*AP+!1.-AP)*EXP(-XLAM*! ! CLDTI M ! MM , 1 1 ) / !3600.*24, ) ) 

1-TLAG  + TAWL)  ) 

2100  IF! FX! MM, II ) .LE. 1 • E-06 ) F X ! MM , 1 1 ) *0. 

2090  CONTINUE 

WRITE! 3, 2114) 

2114  FORMAT!*  * , * POPULATI ON  FRACTION  REMAINING*) 

DO  2102  MM*1,M 

2102  WRITE!3,2101) !FX!MM,II ),II*1,I ) 

2101  FORMAT!*  *,16E8.2) 

DO  3330  MM»1,M 
DO  3320  11*1,1 

3320  POPUL !MM, II »xPOP!MM, 11 )*FX! MM,II ) 

IF!POPUL!MM,II ).LE.0.9)POPUL!MM,II )*0« 

3330  CONTINUE 
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21A0 


3070 


3010 


3020 


5000 

3000 


3040 


3030 


3180 


3190 


8000 


3050 

1551 


NC0UNT*NC0UNT+1 
GO  TO  133 
CONTINUE 
WRITE(3,3070) 

FORMAT! ‘O't  55X, 'EVACUATION  MODEL*/ 

12X, 'POPULATION  MOVED  TO  A DISTANCE  MIDWAY  BETWEEN  THE'# 

2*  RELEASE  POINT  AND  THE  OUTERMOST  RADIUS't 
3'  IN  A DIRECTION  AWAY  FROM  THE  CLOUD*/ 

438X,*THE  POPULATION  OUTSIDE  THE  MIDPOINT  IS  NOT  MOVED'//) 

WRITE! 3t 2072) (AP, XL AH, TLAGfTAWL) 

M*M/2 

DO  3000  II«l,I 
CLDVEL»0, 

DO  3010  JJ«1,J 

CLDV=UBAR (I I,JJ)*FREO( II,JJI 

CLDVEL*CLDV4CLDVEL 

CL0VEL=CLDVEL/100, 

DO  3000  MMxl,M 

1F{CL0VEL.LE.1.0E-05)G0  TO  5000 
CLDTIM(MM, I 1)*X(MM)/CLDVEL 
GO  TO  3000 

CLDT1M(MM,I I )*1.0E  05 
CONTINUE 
DO  3030  MM«l,M 
DO  3030  II»1,I 

FX( MM, II ) = AP+(1.-AP)*EXP(-XLAMM{ ( CLDTI M ( MM, 1 1) / (3500.  *24. ) ) 
l-TLAG+TAWL) ) 

IF(  FX(MM,  in  .LE.l.E-06)FX(MM,in«0. 

CONTINUE 
DO  3190  MM*l,M 
DO  3190  11-1,1 

POPUL(MM, II )sPOP(MM,II )*FX(MM,II  I 
IF(  POPULIMM,  II  ) .LE.C.9)P0PUL('^M,I1  )=0. 

CONTINUE 

MMM-M+1 

M*2*M 

DO  POOO  MMsMMM,M 
DO  8000  11-1,1 
POPULIMM, II  )=POP(MM,II) 

NCOUNT-NCOUNT+1 

GO  TO  133 

CONTINUE 

CALL  EXIT 

STOP 

END 

FUNCTION  SIGZtKLASS, X) 

FUNCTIONS  ARE  OAP/EPA  FITS  TO  SIGMA  ZEES 

IN  EPA  OFFICE  OF  AIR  PROGRAMS  DOCUMENT  NO. AP-26. (USGPO ,WASH. ,0C. 
STOCK  NO. 5503-0015.  PRICE  ONE  DOLLAR.) 

FUNCTION  CALCULATES  THE  STANDARD  DEVIATION  OF  PLU^E  CONCENTRATION 
IN  THE  Z DIRECTION  FOR  STABILITY  KLASS  ( 1 - A,  2 = B , 3 » C, 
4=D,  5=E,  5-F)  AND  THE  DISTANCE  ALONG  THE  CENTERLINE  OF 
THE  PLUME  IN  KILOMETERS,  X. 

THE  STANDARD  DEVIATION  IS  IN  METERS. 

AS  WRITTEN,  A CAPPING  LAYER  (SIGMAX)  UPPER  BOUNDS  SIGZ 


THE  RUN  IS  STOPPED  IF  X(MM)  1:5  LESS  THAN  100.  METERS 


100  IF  (X  - O.l)  110,130,130 
no  WRITEI3,120) 


J 


9 


120  F0RMAT( 1X,*0NE  X IS  ToO  SMALL.  CHECK  X,  XUP 
STOP 

130  XX  = ALOG(X) 

SIGMAX=10000. 

GO  TO  ( 140, 170. 180f 190f20C.210) t KLASS 
140  IF(X-1.5)  150,150,160 

150  SIGZ  = EXP(6. 126788  ♦ XX  * (2.214445  ♦ XX  ♦ 
K-0. 379863  ♦ XX*  (-0. 099597  ))»  )) 

IF  (SIGZ-SIGMAX)  220,220,230 
160  SIGZ=500.  * X**2 

IF  (SIGZ-SIGMAX)  220,220,230 
170  SIGZ  = EXP(4. 686302  ♦ XX  ♦ (1.062550  ♦ XX  * 
IF  (SIGZ-SIGMAX)  220,220,230 
180  SIGZ  * 61.141032  * X ♦♦  .914651 
IF  (SIGZ-SIGMAX)  220,220,230 
190  SIGZ  = FXP(3. 416367  ♦ XX  * (0.729577  + XX  ♦ 
IF  (SIGZ-SIGMAX)  220,220,230 
200  SIGZ  = EXP(  3. 057629  ♦ XX  * (0.679089  XX  ♦ 
IF  (SIGZ-SIGMAX)  220,220,230 
210  SIGZ  * eXP(2. 625488  ♦ XX  ♦ (0.658866  ♦ XX  * 
IF  (SIGZ-SIGMAX)  220,220,230 
230  SIGZ  * SIGMAX 
220  RETURN 
END 


AND  XL(^W  DATA.  • ) 

(-0.041129  ♦ XX  • 

(0.018771) ) ) 

(-0.031207) )) 
(-0.044892) )) 
(-0.054137) ) ) 


o 
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SUBROUT  IN  P RIVCONI  ISO,  XL/kRBD  ,/iCTOEN  ,CRTORG , LATWT,  L,  H,  1 1 
INTEGER  RADIUS(20),SECT0R(20) 

DIMENSION  CC0NC(20,20» 

DIMENSION  AVCROS(7C) ,AAVCRS( 201 

DIMENSION  SOLUB(20> ,FRACT(2C ) , TK I DTH ( 20 , 20 ) ,TCTFIN(20I 
DIMENSION  TVLTIM(20> ,OREACH( 20,20) ,CRC SS ( 20 ,20) ,CCAVE(20), 
ISURFAl 20) ,PLT0EN(20) ,XLRCH( 20) ,I S0(2C) ,XLAMB0{20) 
2,XDEPTH(2C,20,20),OAVF(2C),C(20)  ,CCn20)  ,CRT0RG(20),LATWT(20), 
3 FINALI20,20),NUMB(20) 
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DIMFnSION  ACTDEN(  12tl6f 20 

COM^'ON/ON E/ ACT  ( 12,16,20  ,XXT  I f«,CONCEN( 20.20) 

comkcn/two/raoius,  sector,  hr  ech 
INTEGER  WATFAC(20) 

00  gio 

no  sio  II. i, I 

DO  <310  LL«l,L 

910  ACT(MN, ll,LL)=ACTnEN(KM,II,LL) 

REACt 1,1C)XXT1H 

IF( XXTIM.EO.O. )G0  TO  611 

WR1TEI3,660) 

660  FORVAT( • 1*, 32X, ‘THE  FOLLOWING  SECTION  COMPUTES  THE  *, 

1 'CONCENTRATIONS  OF  VARIOUS  NUCL I HE S • /40X , * I N A RIVER  WITHIN  », 

2 'TFE  AREA  OF  INTEREST  OF  THIS  STUDY'/47X, 

3 'DUE  TO  DEPOSITION  ON  THE  RIVER  SURFACE*/) 

10  FOR^-ATIFIO.O) 

WRITE(3,20XXT1h 

20  FORMAT! 'O', 18X, 'THE  DELAY  TIMF  IN  DAYS  BEFORE  COMPUTING  THE  WATER 
ICONCENTRATIONS  IS  •,F10,0//) 

REACH,  30)TIMOUT,T!MIN,OLTVEL,XINVEL,MNOUT,I  IRECH,NDEPTH,NPROF 
30  F0RMAT(4F5. 1,415) 

WRITE! 3,40) TIMOUT,TIHIN,OUTVEL,XINVEL,NINCUT, I IRECH,NCEPTH,NPR0F 
I IRECH.I IRECH+1 

40  FORMAT!'  ',20X,'THE  TIMEIHRS)  THAT  THE  TIDE  FLOWS  OUT=  •,T85,F5.1/ 
120X,'THE  TIME(HRS)  THAT  THE  TIDE  FLOWS  IN=  ',T85,F5.1/ 

220X,'THE  WATER  VELOC ITY! AVG ) WHEN  THE  TIDE  GOES  OUT=  ',T85,F5.1/ 
320X,'THE  WATER  VELOC ITY! AVG ) WHEN  THE  TIDE  COMES  IN=  ',T85,F5.1/ 
42CX,'THE  TICE  INITIALLY  FLOWS  IN  DIRECTION  ' ,T85 , 1 5 , ' ! Q.QUT , • , 

5' l.IN)'/20X, 'THE  NUMBER  OF  REACHES=  ',T85,I5/ 

620X,*THE  NUMBER  OF  DEPTHS/CROSS  SECTIONAL  PROFILE*  ',T85,I5/ 
720X,'THE  NUMBER  OF  PROF  I LES/REACH*  ',T85,I5//) 

DO  700  LL.1,L 

700  READ! 1,60 )PLTDEN ILL) ,SOLUB( LL) ,FR act  ILL) 

DO  710  IREACH.2, IIRECH 

710  REAr!l,90)  XLRCH I TREACH) , SURF A ! I REACH) , W ATFAC ! I RE ACH) 

DO  720  IREACH.2, IIRECH 
DO  720  NNPROF=l,NPRnF 

720  READ!  1, 230) (TWIDTHI IREACH,NNPRCF ) , I XCEPTH ( I RE ACH, NNPROF , NDEP ) , 
1NDEP*1,N0EPTH) ) 

DO  590  LL.1,L 

WRITE!3,9C0)1S0!LL) ,LATWT!LL) 

EXCESS.O. 

WRITE!3,70)PLTDEN!LL), SOLUBILL) ,FRACT!LL) 

DO  590  IREACH=2,IIRECH 
C!  IREACH).0. 

CCI !IREACH).0. 

CONCEN! 1,LL ).0. 

IFILL.GT.DGO  TO  750 

REAC!  1, 340) NUMB  I IR EACH) , RAO  I US! I RE ACH) , SECTOR! I RE ACH) 

750  CONTINUE 

WRITE! 3, 901 )NUMB! IREACH) 

WRITE!  3,  ICO) XLRCH!  IREACH ) ,SURFA! IREACH) ,WATF AC! IREACH) 

DISTAN.O. 

TVLTIM! 1).0. 

DISNEWO. 

DISTNW.O. 

IFININOUT.LE.O. AND. IREACH, EC. 2)G0  TO  110 
C 

C COMPUTE  THE  TVLTIMC,  THE  TIME  IT  TAKES  FCR  THE  AVERAGE  PACKET  OF  WAT 
C ER  TO  TRAVEL  THE  FULL  LENGTH  OF  THE  REACH 
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• C FOR  THE  TICE  COMING  iNt 


C 

distt=xinvel*timin 

» TVLTIMUREACH)*TIMIN*EXCESS 

DISTAN=OUTVEL*TIMnUT-niSTT+Ci;TVEL*EXCESS 

IF(CISTAN.LE.1.E-10)DISTAN.0. 

O IF(XLRCH(  IREACH)-niSTAN) 12C1 .1110,1200 

1200  TVLTIM(  IREACH)*TVLTIMI  IREACH  )*TI  MOljT 
GO  TO  130 

O C FOR  THE  TIDE  INITIALLY  GOING  CUT: 

110  01STAN=0UTVEL*T1M0UT 

IF(CISTAN.LE.1.P-10)DISTAN=0. 

• IF( XLRCH( IR EACH) -01  STAN) 210, 200, 120 

120  TVLTlMnREACH)»TlMOUT 

130  distnw=distan-xinvel*timin 

• IF(  ClSTNW.LE.l.E-lODI  STNW*C« 

TVLTIMU  REACH  )=TVLTl«f  IREACH  ) + TI  MI  N 
DISNEW=01STNW-»-0UTVEL*TIM0UT 

• IF(XLRCH(  IREACH)-niSN'=W)l60, 150,140 

140  TVLT1M( IREACH)sTVLTIM( IREACH )4TI MOOT 

GO  TO  130 

• 150  TVLTIMI lREACH)=TVLTIMt IREACH)+TIM0UT 

GO  TO  170 

160  01STAN=r)ISNEW-XLRCH(  IREACH) 

• IF(riSTAN,LE.l.E-10)OISTAN*0. 

EXC5SS=DISTAN/0UTVEL 

1F( EXCESS. LE,1.E-20)EXCESS*C. 

• TVLTIMI IREACH)=TVLTIM( IREACH )-EXCESS+CUTVEL*TIMOUT 
IF(TVLTlMnREACH).LE,l.E'-15)TVLTIM(  IREACH)=0, 

GO  TO  170 

• 1110  TVLTIM! IREACH)*TVLTIM( IREACHI+TIMOUT 

EXCESS*0. 

GO  TO  170 

• 1201  TVLTIMI  IREACH)=TVLTIM( IREACH)4TIM0UT 

DIS=DISTAN-XLRCH( IREACH) 

EXCESS=CI5/OUTVEL 

• TVLTIMI IREACH )=TVLTIM( IREACH )-EXCESS 
170  CONTINUE 

GO  TO  220 

• 200  TVLTIMnREACH)»TIMOUT 

GO  TO  170 

210  DISTAN=XLRCH( IREACH) 

• TVLTIMC IREACH)»DISTAN/OUTVEL 
IFCTVLTIMHREACH)  .LE.1.E-15)TVLTIM(IREACH)=0. 

GO  TO  170 

• 220  continue 

DO  320  NNPROF»l,NPROF 
C 

• C determine  THE  VERTICAL  CROSS  SECTIONAL  AREAIIN  SC  METERS)  CF  THE 

C REACH(CR0SS( IREACH) )USING  SIMPSCN'S  RULE 

^ C 

• C USE  AN  ODD  NUMBER  OF  DEPTHS  TO  DIVIDE  THE  WIDTH  INTO  ?**  EQUAL 

C SUPINTERVALS 

^ C 

• depth=noepth 

DEL TAX=TWIOTH( IREACH, NNPROF ) /IDEPTH-1. ) 

NDEP»1 

• IF(NDEPTH,E0.3)GO  TO  260 
IF(N0EPTH.E0.5)G0  TO  270 


o 
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^ 

IF(K0EPTH.E0.7)G0  TO  280 
IF(^DEPTH,E0.9)G0  TO  290 
IF{NDEPTH.EC.11)G0  to  300 
IF( NDEPTH.GT.ll  )GO  TO  240 
1F(VOEPTH.LE.2)GO  to  240 
240  WRITE(3,250) 

260  CRO<;S(  IREACH,NNPROF)  = ( C DELTA  X/3.  ) ♦ I XDEPTH 11  RE ACH,  KNPROF  ,NOEP) 

1+4.*XDEPTH( IREACH.NNPROF.NDEP^l )*XDEPTH( 1REACH,NNPR0F,NDEP+2I )»♦ 
2.0929 
GO  TO  310 

270  CROSS ( I REACH, NNPRnF)=( ( DELT A X/3. ) * ( XDEPTH ( I RE ACH , NNPROF ,NDEP) 
1-»4.*XDEPTH(  IREACH,NNPRnF,N0EP-»l)  + 2.*XCEPTH(  IREACH»N'NPR0P,NDEP+2I 
2 ♦4. ♦XDEPTH ( I REACH,NNPR0F,NDEP^3) ♦XDEPTH ( I REACH,NNPROF ,NDEP+4) I »♦ 

3  .0929 
GO  TO  310 

230  CROSS ( IREACH,NNPROF)  = ( (DELT A X/3. )♦( XDEPTH ( I RE ACH, NNPROF, NOEP) 
144.*XDEPTH( IRE ACH, NNPROF, NDEP  + 1 ) + 2. •XDEPTH (I  RE ACH, NNPROF, NDEP^2I 

2 +4. * XDEPTH ( IPEACH,NNPR0F,NnEP^3 )42.*xnEPTH( IREACH,NNPR0F,NDEP+4) 

3 ♦4.*XDFPTH( IRE ACH,NNPR0F,NDEP+5)+ XDEPTH! IRE ACH,NKPR0F,NDEP+6) ))* 

4 .0929 

GO  TO  310 

290  CROSS! I RE ACH, NNPROF ) = ! ! DELTA X/3. )*( XDEPTH ( I RE ACH , NNPROF , NOEP J 

1  +4. ♦XDEPTH! I RE ACH, NNPROF, NDEP+l )♦ 2. •XDEPTH! IRFACH,NNPRCF,NDEP+2» 

2+4. ♦XDEPTH! IREACH,NNPR0F,NDEP+3)+2.*XDEPTH(1REACH,NNPR0F,NDEP+4) 

3+4. ♦XDEPTH 

3 ! I RE ACH, NNPROF, NDEP  + S ) +2. ♦XDEPTH! I REACH ,NNPRCF ,N0EP^6 ) ♦4.* 

4 XDEPTH! IREACH,NNPR0F,NDEP+7 >♦ XDEPTH ! I REACH,NNPR0F,NCEP+8) ) )^.0929 
GO  TO  310 

300  CROSS! IREACH,NNPRnF»=! ! DELTA X/3. )♦! XDEPTH! I REACH, NNPROF, NDEP) 

1 ♦4. ♦XDEPTH! IREACH,NNPR0F,NDEP+1 ) + 2. ♦XDEPTH! I RE ACH, NNPR0F,NDEP^2» 

2 +4. ♦XDEPTH ( IRFACH, NNPROF, NDEP+3  » + 2. •XDEPTH! I REACH , NNPROF, NDEP+4) 

3 ♦4.+X0FPTH! IPEACH,NNPR0F,NDFP45)42.^xnEPTH{ IPFACH,NNPRCF,NDEP+6» 

4 +4. ♦XDEPTH I IRE ACH, NNPROF, ND EP  + 7 ) +2. * XDEPTH! I RE ACH , NNPRCF , NDEP^8 ) 

5 44. ♦XDEPTH! IRE ACH,NNPROF ,NDEP  + 9 )♦ XDEPTH ! I RE ACH .NNPROF, NDEP+10 ) 1 1 

6 ♦.0929 

310  OREACHI IREACH,NNPROF)=! 9.99997E  02 ) ♦ ! ABS ! CUTVEL-XINVEL ) * 

14. 47019 3E-01)* I CROSS! I REACH, NNPRCF) ) 

320  CONTINUE 

OAVE! IREACH)»0. 

DO  330  NNPROF=1,NPROF 

3 30  GAVE! IRE ACH ) = OREACH! IREACH, NNPROF ) ♦CAVE ! I RE ACHI 
Z=NPROF 

OQAVE! IREACH)=OAVF! IREACHI/Z 

COMPUTE  NUCLIDE  CONCENTRATION 

AVCROS! IREACH)»0. 

AAVCRS!1)*1. 

SURFA<1)*0. 

XLRCH(1)»1. 

TVLTIM!1).0. 

CONCEN! 1,LL )-0. 

COAVE! 1)*0. 

DO  350  NNPROF*!, NPROF 

350  AVCROS! I REACH )*CROSSf IRE ACH, NNPROF )♦ AVCROS! IREACH) 

CC»NPROF 

AAVCRS! 1REACH)*AVCR0SI IREACH)/CC 
MM=PADIUS ! IREACH) 

II=SECTnR!IREACH) 

IF  ! NUMB!  IREACH)  .EO.DACTOEN!  PM,  II  ,LL)*0. 
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CONCEN(  IREACH.LL  )*  ( ACTDEN ( , 1 1 , LL  ) • SLRP A ( I PEACH-l ) ♦.09?9 » 

• 1/  (AAVCRS(  IRE ACH-H*XLRCH(  IRE  ACH-H  *5280.  •?,8317E-02*9.9997E  021 
451  IFINUKBI IREACH).GT.2C»G0  TO  452 

GO  TO  470 

• 452  WR1TE<3»460) 

GO  TO  610 

THE  FIRST  TERM  C IS  THE  CONCENTRATION  OLE  TO  THE  CONTRIBUTION  OF 
THE  PREVIOUS  REACH 


470  C(IREACH)  = (QOAVE(  IREACH-H*C0NCEN(IREACH-1,LL»*E>CP(  (-KLAMBOILLI 
1/24. )*TVLTIM( IRE ACH-l) ) » ♦E XP ( -XL AMBD ( LL ) ♦XXTI M) 
XLL=XLRCH(IREACH-1) 

AAX=AAVCRS( IREACH-1) 

SXX=SURFA( IREACH-ll 
TVX=TVLTIM( IREACH-ll 
OX=CQAVE( IREACH-ll 

THE  SECOND  TERM  CCI  IS  THE  CONCENTRATION  DUE  TO  THE  CONTRIBUTION 
OF  DIRECT  DEPOSITION  AND  THE  PLANT  DISCHARGE 


o 


CCniREACH)  = ( ( ACTDENfMM,  II,LL)*EXP(-XLAmbD(LL)*XXTIM)*SURFA(IREACH 

1 )*.0929) /(  AAVCRSt IRE ACH ) ♦XLRCH ( IPEACH ) *52 80 3048*9 ,99972E  021  )* 

2 PLTOEN(LL) 

CCO^X( IREACH,LL)=(C( IREACH ) +00 K IRE ACH ) ) ♦ ( 1 . /CQAVE 
1( IREACH) I 

WRITE(3,480 J ISO(LLI,LATWT(LL) ,NUMB( IREACH) tCCCNCIIREACH.LLI 
C 

C THE  FOLLOWING  DETERMINES  THE  DRINKING  WATER  CONCENTRATION  FOR 
C THE  VARIOUS  WATER  TREATMENT  FACILITIES  AVAILABLE 
IF(WATFAC(IREACH».E0.1»G0  TO  SCO 
IFUATFAC(IREACH).E0.2)G0  TO  530 
IF(WATFAC(IREACH).GT,2)G0  TO  540 
WRITE(3,490) 

500  FINAL  { IREACH,LL)=CCnNCnREACH,LLI 
IF(wATFAC(IREACHI.EO,l)GO  TO  501 
GO  TO  590 

501  WRITE(3,510JFINAL(IREACH,LLI 
GO  TO  590 

530  FINAL! IRE ACH,LL )*CCONC ( IREACH.LL )*SOLUBI LLI 
WRITE (3, 550) FINAL! IREACH.LLI 
GO  TO  590 

540  FINAL! I REACH, LL ) = FR ACT ! LL ) ♦CCCNC ( I RE ACH, LH ♦SCLUB! LLI 
IF!WATFAC(IREACH).e0.4)G0  TO  570 
WRITE! 3, 560)  FINAL! IREACH, LLI 
GO  TO  590 

570  WRITE! 3, 580)FINAL! IREACH, LLI 
590  CONTINUE 

60  F0RMAT!E15.6,F10.2,F10.2) 

70  FORMAT! 'O', 20X, 'PLANT  DISCHARGE  OF  THIS  NUCLIDE  IN  PCI/SEC*  *, 
1T120,E10.2/10X, 'THE  SOLUBILITY  CF  THIS  NUCLIDE  IN  COLD  WATER*  ‘t 

2 T120, F10.2/20X, 'THE  DISSOLVED  NUCLIDE  FRACTION  PASSING  't 

3 'WATER  TREATMENT  FACILITIES  CF  ALL  TYPES  IS  ' ,T120,F10.2// 1 
90  F0RMAT(F1C.3,E15.6,I5) 

100  FORMAT!'  ',20X,'THE  LENGTH  CF  THE  REACH  IN  MILES*  ',770, 

1 F10.3/20X, 'THE  SURFACE  AREA  CF  THE  RIVER  IN  THIS  REACH*  '.T70f 

2 E15.3/20X, 'THE  WATER  TREATMENT  FACILITY  CLASS  TYPE*  ',T70,I5/  I 
230  F0RMAT(12F6,1) 

250  FORMAT! 'O', 2X, 'ERROR. ..THE  NUMBER  OF  DEPTHS  IS  NOT  ALLOWEO'I 
340  F0RMATI3I5I 
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«60  FORN'AT  (•  0'»  lOX, 'ERROR.  CF  REACHES  EXCEEDS  ?0'///l 

A80  FORMAT! IIX, 'THE  AVERAGE  CCNCENTRATICN  OF  'tAAyAAyZXt 
i'  r;  REACH  ',15,'  IS  •,E16.7,'  PCI/L'/ 

141X, 'BEFORE  WATER  TREATMENT'///) 

A90  FORMAT! ' '.BTX, 'THERE  ARE  NO  WATER  TREATMENT  FACILITIES  IN  THIS*, 
1'  REACH  •///) 

510  FORMAT! ' O', 10X» 'THERE  IS  A CLASS  1 WATER  TREATMENT  FACILITY  IN  *, 
I'THIS  REACH  WHICH  PROV IDES: ' /ICX , • NO  TREATMENT  OTHER  THAN  *, 
2'CHLORINATION' /lOX, 'THERE  IS  NO  REDLCTICN  IN  DISSCLVED  OR  *, 
3'SUSPENDEC  R AD lONUCL IDE S ' /IC X , ' THE  CCNCENTRATICN  AFTER  WATER  *, 
^'TREATMENT  THEREFORE  IS  ' , E I 6. 7,  ' PC  I /L * // ) 

550  FORMAT! 'O', lOX, 'THERE  IS  A CLASS  2 WATER  TREATMENT  FACILITY  IN  *, 
I'THIS  REACH  WHICH  PROV 1 DE S : ' / 1 C X , ' F I LTR AT  I CN  AND  CHLORINATION.  *, 
2'  ALL  SUSPENDED  RADIONUCLIDES  APE  REMOVED  BUT  NO  DISSOLVED 
3'RACIONUCLIDES' /lOX, 'THE  CONCENTRATION  AFTE=  WATER  TREATMENT  *, 

A'  IS',E16.7, 'PCI/L'//) 

560  FORMAT! 'O', lOX, 'THERE  IS  A CLASS  3 WATER  TREATMENT  FACILITY  IN  *, 
I'THIS  REACH  WHICH  PROV IDE S : • /I CX , 'CHEM IC AL  orccESSING.  *, 
2'VARIABLE  AMOUNTS  OF  DISSOLVED  AND  NEARLY  ALL  SUSPENDED 
3'RACIONUCLIDES' /lOX, 'THE  CONCENTRATION  AFTER  WATER  TREATMENT  IS  *, 
AE16.7, 'PCI/L'//) 

580  FORMAT! 'O', lOX, 'THERE  IS  A CLASS  A WATER  TREATMENT  FACILITY  IN  *, 
I'THIS  REACH  WHICH  PROV IDE S : ' /IC X , 'F I LTR ATI CN  PLUS  CHEMICAL  *, 
2*PRCCESSING.  ALL  SUSPENDED  AND  VARIABLE  AMOUNTS  OF  DISSOLVED  *, 
3'RACIONUCLIDES  ARE  REMOVED' /lOX, 'THE  CONCENTRATION  AFTER  *, 
A'WATER  TREATMENT  IS  ' , E16.7, 'PCI /L *//) 

900  FORMAT! 'O', 6AX,AA,AA/) 

901  FORMAT! 'O', 6AX, 'REACH* , 2X, 15 ) 

WR1TEI3,630) 

DO  610  IREACH=2,IIRECH 
TOTFIN! IREACH)«0. 

DO  600  LL=1,L 

600  TOTFIN! I REACH ) = F INAL ! I RE ACH , LL ) ♦TCTF 1 N ! I REACH ) 

WRITE! 3, 620)  NUMB! I REACH)  , TOTFIN! I PEACH) 

620  FORMAT!'  ',20X,'THE  TOTAL  DRINKING  WATER  CONCENTRATION  IN  REACH  *, 
115, 2X,'  IS  ',516.4,'  PCI/L'/) 

630  FORMAT! ' l'//47X, 'SUMMARY  OF  RADIONUCLIDE  CONCENTRATIONS*/) 

610  CONTINUE 

CALL  F ISHR!CCONC,PADIUS, SECTOR,! SO,LATWT,L,I 1RECH,NUMB) 

611  CALL  LAKECS1ACTCEN,L,M,I,LATWT,ISC,SCLUR,FRACT,XLAMB0) 

RETURN 

END 


4 


<» 


139 


SUBROUTINE  LAKEDSIACT.L.f*,  I ,LATHT,I  SCf  SOLUBtFRACT. 

IXLAMBD) 

DIMENSION  SWIMOFI 20) ,BOATDF( 20) ,SKINDFI2C) , ACT ( 12 t 16 t 20 ) t 
lD0SSWB(2C),D0SRWB(?C),DnSSSK(2C),IS0(20) ,LATKT(20) ♦C0N(20), 
2 FINAL(20),S0LUfl(2C),FRACT(20l 
DIMENSION  BOATSKt  20) ,DOSBSK( 20 ,XLAMBDC20) 

INTEGER  SETUP, DESDOS 
INTEGER  WATFAC 
REAC<  1.  lOSETUP 


r>  o o o o 
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10  FORMATdS) 

DO  482  JJJ»1, SETUP 

1F( SETUP. EQ.O)GO  TO  481 

READ « 1.20) SURF, SWIM, BOAT, VOLUME 

WRITEO,  11  )SURF,  SWIM,  BOAT,  VOLUME 

11  FORMAT(*0'//20X,*LAKF  SURFACE  AREAISC  FT)*  •,T50,E15.2/ 

1 20X, 'TIMEIHRS)  SPENT  SWIMMING*  •,T5C,F10.l/ 

2 20X, 'TIMEIHRS)  SPENT  BOATING*  •,T5C,F10.1/ 

3 20X,*LAKE  VOLUME(GALLONS)=  • ,T50,E15.2//) 

20  F0RMAT(E15.7,2F1C.1,E15.7) 

REACH,  30  )DESD0S,LK1*LK2 

30  FORMAT! 15, A4,A4) 

DO  50  LL*l,L 

IF!  JJJ.GT.DGO  TO  31 

READ! 1,40)SWIMDF!LL),SKINDF!LL) 

31  CONTINUE 

BOATDFILL )* . 5*SWIMDFf LL) 

50  BOATSKILL  )».5*SKINOF!LLI 
40  F0RMAT!2E15.7) 

READ!  1,60  )MM-,  II, WATFAC 
60  FORMATI3I5) 

IF!  JJJ.GT.DGO  TO  61 
DO  14  LL*1,L 

WRITE(3,16) ISO(LL),LATHTILL) 

WRITE!3,12) ACT!MM,II,LL) 

14  WRITE!3, 13)SWIMOF!LL), SKINDF!LL),BCATDF!LL)  ,BCATSK!LL) 

13  FOR*'AT( 'OV/lOX, ‘WHOLE  BODY  SWIMMING  DCSE  FACTOR ! MREM/HR  PER  •, 
I'PCI/D*  STTE.EIS.B/IOX.^SKIN  SWIMMING  DOSE  FACTOR*, 

2 '(MREM/HR  PER  PCI/D*  ',  T75  , E15 . 3/lCX  ,' WHOLE  BODY  BOATING  •, 

3 'DCSE  FACTOR(MREM/HR  PER  PCI/D*  • ,T75,E15.3/10X, 

4 'SKIN  BOATING  DOSE  FACTOR  ( MREM/HR  PER  PCI/D*  * , T75 , E15.3// ) 

16  FORMAT ( •0'/2X, 'FOR:  *,44, A4/) 

COMPUTE  WATER  CONCENTRATIONS  AND  DCSES  IN  MREM 


FIX  THE  LOCATION  OF  THE  LAKE 


61  CONTINUE  ] 

DO  70  LL*1,L  1 

CON!LL  )*(  ACTCMM,  II,LL)*.092<3*SL'RF)  / (VCLUMEM3.78541) 

12  FORMAT!  •0'/20X, 'ACTIVITY  DEN  S I T Y ! PC  I /M>**2  ) = ' ,E15.2/) 

DOS 5W BILL )*C0N!LL)MSWIMDF!LL)*SWIM 
DOSEWB(LL  )*CON(LDMBOATDF(LL)*BOAT 
DOSSSK (LL )*CON(LL)MSKINnF(LL)*SWIM 
70  DOSESK(LL ) =CON « LL) *BOATSK ! LL) MBOAT 
WRITE(3,80)LK1,LK? 

80  F0RMAT(*1'///////23X,'THE  FOLLOWING  IS  A LISTING  OF  RADIONUCLIDE', 

1'  CONCENTRATIONS  IN  LAKE  ',2A4,*BY  ISOTOPE'///)  I 

DO  100  LL*1,L 

100  WRITE(3,90) ISO(LL) ,LATWT(LL) ,CON!LL) 

90  FORMAT! *0',39X, 'THE  CONOENTRATION  OF  * ,A4,2X, A4, * IS  ',E15.7, 

1 2X,*  PCI/L  */) 

TOTC0N*0. 

DO  110  LL*1,L 
110  TOTCON*CON(  LD+TOTCON 

WRITE! 3,120)LK1,LK2,TOTCON 


o. 


IFfCESDOS.EC.DGO  TO  170 
DO  UO  LL«ltL 

1F(COSSWB(LL).EO  .0.)G0  TO  1«0 

WRITEO,  130)LKl,LK2,SWIM,rSC(Ll  ) tLATWTCLL)  ♦DOSSWBfLLI 
lAO  CONTINUE 

130  FORKATC  ‘.lAXt'THE  KHOLE  BODY  DCSE  OtE  TC  SWI^'KING  IN  LAKE  •* 

1 2AA,'  FOR  • ,F5.1,2X, 'HOURS  DUE  TC  • ,A4, AA,3X, • IS* ,E15.7,2X, 

2 'MREH'/I 
TOTCOS-0, 

00  150  LL*ltL 

150  TOTCOS=TOTOOS4DOSSWBILL) 

WRITE! 3t IfO )LK1,LK2,T0T00S 

160  FORMAT!  • '.igx.'THE  TOTAL  WHCLE  BODY  DCSE  DUE  TO  SWlf'WING** 
I'FRCM  ALL  NUCLIDES  IN  LAKE  •,2A4,2X,»  IS  * »E15.2,2X, 'MREM* » 
IFICESDOS.EO.OIGO  TO  220 
170  CONTINUE 

DO  180  LL«1,L 

IF!COSSSK!LL  ).E0.0.>  GO  TO  180 

WRITE! 3. 190)LKl,LK2tSWIM,ISC!LL) tLATWT!LL)  .DOSSSKILLl 
180  CONTINUE 

190  FORMAT!'  '.ITX.'THE  SKIN  DOSE  DUE  TO  SWIMMING  IN  LAKE  *, 

1 2AA,'  FOR  ' ,F5,1,2X, 'HOURS  DUE  TO  • ,A4, A4,3X,  • IS  •, 

2 E15.7,2X,'MIREM'/) 

T0TC0S*0. 

DO  200  LL»lfL 

200  TOTCOS=TOTDnS+DOSSSKCLL) 

WRITE! 3,2 10) LK1,LK2.T0T00S 

210  FORMAT!'  ',22X,'THE  TOTAL  SKIN  DOSE  DUE  TC  SWIMMING  FROM  ALLS 

1 ' NUCLIDES  IN  LAKE  ',2A4,2X,'  IS  • ,E1 5. 3 ,2X , 'MREM* ) 

220  CONTINUE 

IFIBOAT.EO.O. )G0  TO  320 
IF!CESD0S.E0.1)G0  TO  270 
DO  230  LL*1,L 

IF!COSBWB!LL).EO.O.)GO  TO  230 

WRITE!3,240)LK1,LK2,B0AT,ISC!LL) ,LATWT!LL) ,D0SBWB(LLI 
230  CONTINUE 

240  FORMAT!'  ',14X,'THE  WHOLE  BODY  DOSE  DUE  TO  BOATING  ON  LAKE 
12A4,'  FOR  ' ,F5.1,2X, 'HOURS  DUE  TC  ',A4,A4,3X,'  IS  •,E15«3,2X, 

2 'MREM'I 
T0TC0S=0. 

DO  250  LL*1,L 

250  TOTCOS=TOTDOS+OOSBWB!LL) 

WRITE! 3,260  )LK1,LK 2, TOTDOS 

260  FORMAT!'  ',20X,'THE  TOTAL  WHOLE  BODY  DCSE  DUE  TC  BOATING  S 
1 'FROM  ALL  NUCLIDES  IN  LAKE  ',2A4,2X,'  IS  ' , E15.2 , * MREM* ) 
1F!CESDOS.EO.O)GO  TO  320 
270  CONTINUE 

DO  280  LL*1,L 

IF!COSBSK!LL ).E0.0.)G0  TO  280 

WRITE! 3, 290)LK1,LK2,SWIM, ISC!LL) ,LATWT!LL) ,OCSBSK!LLI 
280  CONTINUE 

290  FORMAT!'  •,18X,'THF  SKIN  DOSE  DUE  TC  BCATING  IN  LAKE  •,2A4, 

1 • FOR  ' rF5.1,2X, 'HOURS  DUE  TC  ',2A4,3X,'IS  ',E15.3,'  MREM*) 
T0TCOS*0« 

DO  300  LL»1,L 

300  TOTCOS=TOTDOS*DOSBSK!LL) 

WRITE! 3,310 )LK1,LK 2, TOTDOS 

310  FORMAT!*  ♦,22X,'THE  TOTAL  SKIN  DCSE  DUE  TC  BOATING  FROM  ALL  *, 
1 * NUCLIDES  IN  LAKE  *,2A4,2X,*  IS  * , El 5. 2 , 2X , * mrfm* ) 

320  CONTINUE 
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SUM*0. 

WRITE(3«5C0)WATFAC 

500  fqrvatc  • , sox, •treat>«ent  facility  class*  **I5I 
DO  ^00  LL*l,L 
IF(KATFAC.E0.1)G0  to  350 
IF(wATFAC.EC.2)G0  TO  360 
IF(kATFAC.GT.2)G0  TO  370 
WRITE(3,3A0) 

GO  TO  A81 

350  FINAL«LL»=C0N(LLI 
GO  TO  380 

360  FINAL(LL)*CON(LL)*SOLUB(LL) 

GO  TO  380 

370  FINALCLL )=FRACT(LL )*CON(LL)*SCLUB(LL) 

380  CONTINUE 

SUM=SUM^FINAL(LLI 
400  CONTINUE 

340  FORMAT!  • SBTX, 'THERE  IS  NO  V.ATER  TREATMENT*) 

WRITE! 3,410)WATFAC 

410  FORMAT!'  ',28X, 'AFTER  WATER  TREATMENT,  THE  FCLLCWING 

1 'CONCENTRATIONS  ARE  PRESENT  IN  THE  WATER'/ 

2 38X,*THE  WATER  WAS  TREATED  IN  A CLASS  *, 1 5 , 2X , 'TREATMENT  », 

3 'FACILITY'  ) 

WRITE!3,420) 

420  FORMAT! 'O', 37X, * ISOTOPE *,10X, 'ATOMIC  WEIGHT' ,10X, 'CONCENTRATION', 

1 ' !PCI/LI '///) 

DO  430  LL»1,L 

4  30  WRITE! 3, 4 40  ) ISO! LL) ,LATWT!LL) ,F INAL! LL) 

440  FORMAT!'  ' , 40X, A4, 14X , A4 , 1 6 X ,F15. 3 ) 

460  WRITE! 3,470»SUM 

470  FORMAT! '0'///39X, 'THE  TOTAL  RADIONUCLIDE  CONCENTRATION  IS 

1 E15.7,'  PCI/L'/39X, 'DECAY  EN  ROUTE  TC  WATER  USAGE  IS  NCT 

2 'CONSIDERED.'////) 

481  CONTINUE 

CALL  FISHL!CON,MM,II ,1S0,LATWT,L) 

482  CONTINUE 

CALL  MILK!L,M, I, IS0,LATWT,XLAVBD,CCN,FRACT,S0LUB) 

480  RETURN 
END 

SUBROUTINE  F I SHL ! CON , MM , 1 1 , 1 SO ,L ATWT , L ) 

DIMENSION  CON! 20), ISO! 20 ) ,LA TWT ! 2C ) ,CR ! 20 ) ,F I SHCN( 20 ) 
C0MM0N/THREE/CFI20) 

DO  10  LL«1,L 

IF!CF!LL ) .EQ.O. )G0  TO  100 
10  CONTINUE 
20  FORMATIFIO.O) 

FISHT*0. 

DO  30  LL*1,L 

FISHCN!LL )*CON ! LL ) ♦CF!  LL ) 

30  FISHT*FISHT*FISHCN!LL) 

WRIT£!3,40) 

40  FORMAT  ! •-•///////  / •♦♦♦♦♦♦♦♦♦•♦♦♦♦****i******#********************i  , 
X I i^««*******«****«4i**«***«**«i»*******«****4>«*«#*44i**«*«4t4i«*«***//  / ) 

DO  42  LL*1,L 

IFILL.GT.DGO  TO  42 

WRIT£!3,41) ISD!LL),LATWT!LL) ,CFILL) 

42  CONTINUE 

41  FORMAT!*  * , 'CONCENTRAT ION  F ACTCR S * /48X , A4 , A4 ,2X , • « 

1 FIC.O,'  PCI/KG  PER  PCI/L') 

WRIT£!3,50) 
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50  FORVAK '-'////////AIX, 'NUCLIDE  CCNCENTRATI CN  IN  FISH  LOCATED*, 

• 1 • IN  THE  LAKE'/////) 

DO  60  LL*ltL 

IFIFISHCN(LL).EQ.O.)GO  TO  6C 
^ WR1TE(3,7C)  ISOILD.LATWTILL)  ,FISHCN(LL) 

60  CONTINUE 

70  FCRMATC  ',46X,A«,AA,'  = ' , 1 C X , E 15, 2 , ' PC  I /KG  ' I 

• WRITE! 3,80)FISHT 

80  FORMAT! '-'////A2X, 'THE  TOTAL  CONCENTRATION  = ' ,E15.2,2X  ,'PCI/KG* 

_ Iff) 

• WRITEt3,90) 

90  FORMAT ( ' 0'/ /'♦*♦♦♦♦♦*»*♦♦*♦*♦*♦♦*♦♦•♦♦*♦*♦♦♦♦♦*♦*♦♦*♦♦♦*♦*♦**♦*♦•  , 
X • n***.m*******.*^n^i^^'ti****n*********^*n*********it***it**'^*i^******  f f f I 

• 100  RETURN 

END 

SUBROUTINE  F I SHR ! CONCEN , R AD  I LS , SECTOR , I SO , L ATWT , L , 1 1 RECH.NUHB) 

• DIMENSION  RADIUS! 20) , SECTOR! 20 ,CCNCEN!20, 20) ,lS0f20) , 

1 LATWT!  20,  FISHCN!  20 . 20  ) , NUMB ! 20  ) , F I SHT ! 20 ) 

COMMON/THREE/CF!  201 

• DO  10  LL*1,L 
REAC(1,20)CF!LL) 

IF!CF!LL ) .EO.O. )G0  TO  130 

• 10  CONTINUE 

20  FORMAT!F10.0) 

DO  AO  IREACH=2, IIRECH 

• FISHT! IREACH)*0. 

DO  AO  LL*1,L 

FlShCN! IREACH,LL )=CONCEN! IREACH,LL)+CF!LL) 

• AO  FIShT! 1REACH)=FISHT!IREACH)4FISHCN!IREACH,LLI 

WRITE!3,60) 

WRITE!3,50) 

• 50  FORMAT!'  ' , AIX, 'NUCL IDE  CONCENTRATION  IN  FISH  LOCATED  ', 

1 'IN  THE  RIVER' ) 

60  FORMAT ! •-•//////•*♦♦♦♦♦♦♦*♦*♦***♦♦•♦♦♦♦♦*•**♦♦♦♦♦♦***♦**♦♦♦♦♦♦♦♦» , 

DO  100  IREACH=2,- IIRECH 
WRITE(3,70NUMB!  IREACH) 

• DO  90  LL*:1,L 
1F1FISHCN!IREACH,LL).E0.C.)G0  TO  9C 

WRITE! 3,80) ISO! LL) ,LATWT!LL) ,F I SHCN! IREACH, LL) 

• 90  CONTINUE 

100  WRITE! 3, llOIFISHT! IREACH) 

70  FOR-^ATl'  ',61X,'RRACH' ,2X,I5) 

O flO  FORMAT!'  ' , A6X, A A, AA , ' = ' ,1CX  ,E15. 2 , ' PCI/KG') 

110  FORMAT!'  ',A2X,'THE  TOTAL  CCNCENTRATICN  * ',E15.2,2X, 

1 'PCI/KG') 

• WRITE!3,120) 

120  FORMAT ! ' 0* // * **************************************************** , 

X • #«*«*«**  **4t****«4i*  *******  «*««*«  *************  **4*«*i,**«i,«****,  / / ) 

• 130  RETURN 

END 


o 
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SUBROUTINE  MlLKCLtH,  1 1 I SC  ,L  A TV.T  , XL  A KBD  ,CCN  ,FR  ACT,  SOLUB) 
C0MN'0N/0NE/ACT(  12,16,20)  , XXT I R .CCNCEN  ( 20,20) 

INTEGER  RR (20), SS{20) 

INTEGER  DRINK(  192)  ,RADIUS(  1'32)  , SECTOR  ( 192  ), PAY, TOTOAY , WATFAC  (192) 
REAL  LAMEFF(20),  ID 
DIMENSION  C0WL(20) 

DIMENSION  ISO)  20),LATVIT(  20  , XLAMBD(2C)  ,FRACT  ( 20)  , S0LU3  ( 20 ) , CON  (20) 
1 ,CC( 12, 16, 20) ,SD( 20 , TOTAL! 12,16) ,BF(20) 

COMMON/TWO/P.R,  SS,  IIRECH 
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REACtl,3C)NUMBER 

WRITE(3,410) 

WRITE(3,qO)NUHBER 

30  FORMATdS) 

IF(NUHBEP..E0,-1  >G0  TO  400 
RFAr(l,lC)OF,R,YF 
WRITE<  3, 1C0)0F,R»YF 

10  F0RMAT(3F10.2» 

REAC( 1.22) XLD,TAREA. VOL 
WRITE(3,120)XL0,TAREA,V0L 
IF(MjM8ER,E0.0)Ga  TO  170 
DO  40  J=l, NUMBER 

40  REAC(1»20)DRINK( JI.RAOIUSi Jl .SECTOR (J) .WATFACIJ) 

WRITE(3.14l) 

DO  130  J*l. NUMBER 

130  WR1TE(3.140)0RINK( J).RADIUSC J) .SECTCR( J)  .WATFACfJ) 

20  F0RMAT(4I5) 

31  FnRMAT( I5.F10.2I 

21  F0R»'AT<2F10.2) 

22  FnRMAT(3F10,2) 

170  REACd.SOTOTOAY.INC 

180  WRITE( 3, 150)TnT0AY,INC 
REAr(1.31)LEAVE.BI0 
WRITE(3.230)LEAVE.BI0 

50  F0RMATI215) 

DO  501  LL*l,L 

501  REAC( 1.502)BF(LL) 

DO  889  LL-l.L 

889  WRrTE(3. llO)BF(LLI.ISn(LL).LAThT«LLI 

502  FORmaT(E10.2) 

DO  270  LL«1.L 

270  REAr(1.280)SD(LLI 
DO  183  LLal.L 

183  WRITE( 3.290)SO(LL) .ISGILL) .LAThTILLI 

280  F0R»'AT(E15.1» 

REAC(1.888)CBIO 

WRITE(3.9C1)CBI0 

901  FORVAK*  »,  lOX. 'BIOLOGICAL  HALF  LIFE  IN  THE  CCW*  •.F10.2/» 

888  FORMATIFIO.2) 

WRITE! 3, 160)XXTIM 

160  F0RMAT(/1CX,  • INITIAL  TIMEIDAYS)  AFTER  ACC IOENT»' .2X.F5. I ) 

90  FORVATI •-•////////5X. I5,2X, •=  THE  NUMBER  OF  MESH  POINTS(MM,in  *. 

1 'THAT  USE  EXPOSED  RIVER  OR  LAKE  WATER  FOP  CCW  DRINKING  WATER*/! 

100  FOR»'AT( 'O'.  5X,F10.2.2X. 'sFRFSH  FORAGE  INGESTICN  BY  *, 

1*  THE  DAIRY  COW  I KG/D AY ) • /5X , F 1C . 2 . 2 X . • =DE PCS ITTCN  *, 

2 'RETENTION  FACTOR  (FRACTION  OF  DEPOSITION  RETAINED  BY  THE  *. 

3 'FORAGE  CROP)'/5X.F10.2,2X, *=FGRAGE  CROP  YI ELD (KG/M**2 ) •/ ! 

110  FOR^'ATC  '.  5X,F10.2.2X. '^FORAGF  CPCP  CCMCENTRATION  FACTOR  *. 

1  'FCR  •,A4,A4,'  (PCI/KG  PER  PCI/KG  SCIL)*/) 

120  FOR*'AT(  *0'.  5X.F10.2,2X. '*INGESTICN  CF  DRINKING  WATER  BY  COW  *, 

1 • (L/DAY) •/5X,F10.2.2Xt •=WATERING  TROUGH  AREA  FROM  AN  UNDERGROU*. 

2 'NC  WELL(FT**2)*/5X, FIO. 2. 2X. 'VOLUME  CF  WATERING  TROUGH*. 

3'(GALL0NS)'//) 

140  FORMAT! 'O'. 5X, 15. '^DRINKING  WATER  TYPE  ( 1=RI VER .2*LAKE ) */ 

1 5X.I5. '-RADIUS  OF  MESH  POINT  USING  THE  A5CVE  DRINKING  WATER*/ 

2 5X.  15, '-SECTOR  OF  MESH  POINT  USI'-G  THE  AVCVE  DRINKING  WATER*/ 

3 5X, 15, '-SERVICING  WATER  TREATMENT  FACILITY  TYPE  *, 

4 '(O-NONE, 1-CLASS  1,2-CLASS  2,3-CLASS  3,4-CLASS  4*/| 

141  FORMAT! *0'////44X, 'LOCATIONS  USING  OTHER  THAN  UNDERGROUND  WELLS* 

1 //! 
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150  FOR^'AK ‘OV/SX,  I5,2)f, ‘-MAXIMCr  Nt'.BER  CF  PAYS  WHEN  COMPUTING  •# 

1 ‘MILK  CONCENTRATIONS  (DAYS  AFTER  ACC  I QFNT ) • /EV , I S, 2X , 

2 •* INCREMENTSCDAyS)  BETWEEN  CCRPLTING  SUCCESSIVE  CONCENTRATIONS*/) 
230  FORMAT (• 0* t 5X, I5t2X, 'sEVACUATIGN  MCOE  (0=NC  EVACUATION,  *, 

1 •1  = EVACUATI0N  AFTER  ONE  PAV'*S  E XPCStRE ) • /5 X ,F10 .2 , 2X , 

2 •=PI0L0GICAL  HALF  LIFE  OF  THE  FORAGE  CROP  IN  THE  CPW(DAYS)*/) 

290  FORMAT! 'O', 5X,E15.1,2X, 'sCOEFFICIENT  CF  TRANSFFP  FROM  DIET  TO  *, 

1*MILK{PCI/L  PER  PCI/DAY>  FOR  *,44,44) 

NXXTIM=XXTIM 

DO  60  0AY=NXXTIM,T0T0AY,INC 
DO  75 

DO  75  n«i,i 
TOTALIMM, II )*0. 

DO  70  LL»1,L 

LAMEFFILL )* XL AMRD( LL )♦ . 693/B 10 

RINC=INC 

XDAYaOAY 

COWL (LL)= XL AMBDCLL ) +.693/06 10 
IFlCOWLtLL) .LE.1.E-20)C0WH  LL)»0. 

0MAX*LAMEFF{LL)*XDAY 
IF(CAY.E0.NXXTIM)G0  TO  79 
IF'LEAVE.EO.DGO  TO  320 
IF! CMAX.GE.40. )DUMMY*0. 

IFICMAX  .GE.AO.IGO  TO  81 

DUMMY=CD!MM, II,LL)*( 1.-EXPI-CCWL(LL)*RINC) ) 

01  CONTINUF 

iFtCUMMY.LE. l.E-20)0UMMY*C. 

IF!  CAY.GT.NXXTIM.AND.LEAVE.EO.OGC  TO  312 
79  00  72  J*lfNUMBER 

71  IF!CRINK! J).E0.1.AN0.RAD1US! J».EC.MM.AND.SECT0R!J).FC.II)G0  TO  190 
IF!CRINK! J) .EQ.2.AN0.RAD1US! J).EC.MM.AND.SECTCR(J).EO.n )G0  TO  221 

72  CONTINUE 

DRINKING  WATER  IS  UNDERGROUND , FROM  RIVER, OR  FROM  LAKE 

191  IFICMAX. GT. 40. )W=0. 

IFICMAX. GT. 40. )G0  TO  260 

W^ACTIMM, II ,LL )*TAREA*.C929*XLD/ ! VOL ♦3.785412)  ♦EXP<-OMAX) 
IF!M.LE.1.E-30)W=0. 

GO  TO  260 

190  IF! WATFACIJ  ).E0.2»GO  TO  210 
IFIWATFACIJ  ).GT.2)G0  TO  220 
DO  700  IREACH=1, IIRECH 

IFIRAOIUSIJ ) .E0.RR!IREACH).AND.SECTCR( J) .EC.SSIIREACH) )G0  TO  701 

700  CONTINUE 

CALL  REACHIRAOIUS, SECTOR) 

701  IFICMAX. GT. 40. )W=0. 

IFrCMAX.GT.40. )G0  TO  260 

W=CCNCEN! IREACH,LL)*XLO  ♦EXPI-DMAX) 

IFIW.LE.1.E-30)W*0. 

GO  TO  260 

210  DO  702  IREACH*!, IIRECH 

IF!R4DIUS!J ).EO.RR! IRE ACH). AND. SECTOR !J). EC. SS(IREACH) )G0  TO  703 

702  CONTINUE 

CALL  REACHIRADIUS, SECTOR) 

703  IFICMAX. GT. 40. )W*0. 

IF! CMAX.GT,40. )G0  TO  260 

W*CCNCEN!  IREACH,LL  )4XLD*SCLL'B!LL  ) ♦EXPI-DMAX) 

IF!W.LE.1.E-30)W«0. 

GO  TO  260 
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220  DO  704  IREACH*lt ITRECH 


X 9 IFtPADIUS(J j.EO.RRITREACHj.AND.SECTORt J).EC.SSCIREACH) )G0  TO  705 

704  CONTINUE 

CALL  REACHIRADIUS, SECTOR) 

C 705  IF( DMAX.GT.40. )W»0. 

IF(CMAX.GT.40.)G0  TO  260 

W = CCNCEN{  IREACH,LL)*XLD*FOACT(Ll)*SCLt'B(LL)  ♦EXPI-Of'AX) 
IF(k..LE.l.E-30)W»0. 

GO  TO  260 

221  IF(WATFACIJ ).EC.2)G0  TO  240 
IF(v,ATFAC(J  ) .GT.2)G0  TO  250 
IF(CMAX.GT,40« )W»0. 

IF(CMAX.GT,40. )G0  TO  260 
W=CCN(LL)*XLO  ♦EXPf-OMAX) 

IF(W.LE.l.E-30)W-0, 

GO  TO  260 

240  IF{CMAX.GT.40. )W*0. 

1F(CMAX.GT,40, )G0  TO  260 
W=CCN(LL )*XLD*S0LUB(LL)*EXP|-DPAX) 

1F(W.LE.1.E-30)W*0. 

GO  TO  260 

250  IF(CMAX.GT.40.)W»0, 

IF! CMAX.GT.40. )G0  TO  260 

W=CCN(LL)*XLD*FRACT(LL )*SCLUB( LL>*EXP(-DMAX) 

IF(W.LE.I.E-30)W=0. 

260  CONTINUE 

COMPUTE  MILK  CONCENTRATIONS 

ID=CF*(P*ACT(MM,II,LL)/YF+BF(LL)*ACT(MM,II ,LL)/224.»+H 
IF(  ID.LE. l.E-50) ID=0. 

C0(MM,II,LL)»SD(LL)*ID 

IF(CD(MWt  II  ,LL)  .LF.l.s=-2C)Cr(MM,II,LL)=0. 

IF(CAY.EQ.NXXTTM)G0  to  3CC 
CDIMM, 1 1 ,LL )=CD(MM, I I,LL )»EXPI-DMAX) 

IF(CO(MM, II .LL ).LE.1.E-2C)C0(MM,II ,LL)*0. 

300  IF(LEAVE.EQ.O.mND.OAY.GT,NXXTIK)GC  TC  312 

310  CONTINUE 
GO  TO  311 

312  CD(MM, II,LL )=DUMMY*EXPC-LAmeFF(LL)*RINC) 

IF(CD(MM, 11,LL).LE.1.E-20)C0(KM,II ,LL)»0. 

311  CONTINUE 

IF( PAY.EO.NXXTIMIGO  TO  322 
IF(LEAVE.EO.O)GO  to  321 

320  CONTINUE 

CD  (MM,  I I,LL  )*CD(MM,II,LL  )*EXP(-CCViL(LLI*RINC» 

IF(CD(MM, II,LL).LE.1.E-20)CD(KM,II ,LL)*0. 

GO  TO  321  I 

322  XMAX=LAMFFF(LL)*XXTIM 

IF( XMAX.GE.40. JCD( MM, 1 1 ,LL 1*0. 

IF( XMAX.GE.40. )G0  TO  R3 
CO(MM, II,LL )*CO(MM,II,LL)*EXP(-XMAXI 
83  CONTINUE 

IF(CD(MM, II ,LL) .LF.1.F-20)CD(MM,II,LL)*0. 

321  TOTAL  (MM,  II  )»TnTAL(MM,in4CD(MM,II  ,LL) 

I F( TOTAL (MM, II ) .LE. 1 .E-20 ) TOTAL ( MM , 1 1 ) .Q. 

70  CONTINUE 
75  CONTINUE 

DO  350  LL*1,L 

WRITF(3,340)IS0(LL),LATWT(LL),DAV 
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WRITE(3,330) 

• DO  350 

350  WRITE( 3.360 ) (MM, (COIMM, II.LL) ) ) 

^ WRITE( 3,A20)DAY 

“ WRITE(3,330) 

DO  430 

430  WRITE(3,360)(MM,(T0T&L(HM,II),II*1,I)) 

• SUM=0. 

DO  370  LL-l.L 
DO  370  MM»1,M 

• DO  370 

IF(CD(MM, lI,LL).LE.l.E-10)CD(MM,n .LL)*0. 

370  SUM=SUM*CC(MM, I1,LL) 

• IF( SUM.EQ.O. )G0  TO  380 
60  CONTINUE 

GO  TO  400 

• 330  FORM&K ' 0*///lX, •MM',3X» 'N • , 5X , • NNE • ,4 X, • NE • , 5X, • ENE • ,5X , • E * , 5X , 

1 'ESE* .4X,'SE*,5X, •SSE',5X, •S*,5X, 'SS^'.^X, 

2 'SW , 5X, 'WSW* ,4X, ‘W*. 6X, 'WNW* ,4X, ‘NW* ,5X, ‘NNW*/! 

• 340  FORMAT! •0'/////44X,A4,A4»3x, 'CCNCENTRATIONSfPCI/D  AFTER  *. 

1 14, 2X, 'DAYS'//) 

360  FORMATC  I2,2X,16E7.2> 

• 380  WRITE(3t390» 

390  FORMAT! ‘O'/Z/Z/lOX, 'PROGRAM  TERMINATED.  ALL  MILK  CONCENTRATIONS  »• 

1 • ARE  LESS  THAN  lE-10  PCI/L*) 

• 410  FORMAT ! • 1 •////////••*♦♦♦**♦*♦♦*♦♦**♦♦♦♦♦♦*♦*♦♦♦♦♦♦*♦♦*♦♦♦♦♦♦♦♦♦*•, 

1 t ****4c*]^«*:^* *4141  #*****:4i #««#*«*#«*** 4 :» **«*«**«* ***4>*****:^i»«4i4i*4i**«  , 

2 « ******** ****#*«4i***«  I / 

9 3 .*************************»*******************•, 

X • *************************** *4c:4***4i*4 »*#***********************.  , 

4* ********************* ///// 

• 5 49X, 'RADIONUCLIDE  CONCENTRATICNS  IN  MILK'///I 

420  FORMAT(*0'/////43X, 'TOTAL  CCNC'^NTR  ATT  CN  ( PC  I /L  » AFTER  ',I4,2X, 

I'DAYS'//' ******************** *♦♦*♦** ♦♦♦**♦♦♦*****♦***♦****♦****• , 

O 2 ' *************************************************************** , 

3 •*«***///) 

400  RETURN 

• END 

SUBROUTINE  REACH!RADIUS, SECTOR ) 


COMPUTES  THE  CLOSEST  REACH  PROVIDING  WATER 


DIMENSION  XRADIU(1>,XSECT0! 1).XRR(12) ,XSS!16I 
COM*'ON/TWO/RR,  SS,  I IRECH 

INTEGER  RADIUS! 1), SECTOR! II, RR! 12)  ,SS( 15 1 
XMINalOOCO. 

XRAEIU! J)=RADIUS!JI 

XSECTOIJ)=SECTOR( Jl 

DO  2 IREACH=2, IIRECH 

IIRECH*IIRECH+1 

XRR ! IREACH)*RR! IREACHI 

XSS !IREACH)=SS! IREACHI 

A=ABS!XRAOIU! J )-XRR( IREACHI I 

B=ARS!XSECTO( J )-XSS! IREACHI I 

C=SCRT! A**2*B**2I 

IF!C.GT.XMIN)GO  TO  2 

XMIN»C 

IRMIN»IREACH 
2 CONTINUE 

1REACH*IRMIN 
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WRITE( 3,3)IREACH 

3 FORyAT( •0»//5X,*THE  CLOSEST  RIVER  REACH  TC  PROVIDE  WATER  S 
1 'IS  NUMBER  •t2X,I5//» 

RETURN 

END 
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